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This  thesis  presents  an  application  of  the  Grant-Madsen-Glenn  bottom  boundary 
layer  model  (Grant  and  Madsen,  1979;  Glenn  and  Grant,  1987)  to  predictions  of 
sediment  transport  on  the  continental  shelf.  The  analysis  is  a  two-stage  process. 
Via  numerical  experiment,  we  explore  the  sensitivity  of  sediment  transport  to  vari¬ 
ations  in  model  parameters  and  assumptions.  A  notable  result  is  the  enhancement 
of  suspended  sediment  stratification  due  to  wave  boundary  layer  effects.  When 
sediment  stratification  is  neglected  under  conditions  of  large  wave  bottom  veloc¬ 
ities  (i.e.  U6  >  concentration  predictions  can  be  more  than  an  order  of 

magnitude  higher  than  any  observed  during  storm  conditions  on  the  continental 
shelf. 

A  number  of  limitations  to  application  emerged  from  the  analysis.  Solutions 
to  the  stratified  model  are  not  uniquely  determined  under  a  number  of  cases  of 
interest,  potentially  leading  to  gross  inaccuracies  in  the  prediction  of  sediment  load 
and  transport.  Load  and  sediment  transport  in  the  outer  Ekman  Layer,  beyond  the 
region  of  emphasis  for  the  model,  can  be  as  large  or  larger  than  the  near-bottom 
estimates  in  some  cases;  such  results  suggest  directions  for  improvements  in  the 
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theoretical  model. 

In  the  second  step  of  the  analysis,  we  test  the  ability  of  the  model  to  make 
predictions  of  net  sediment  transport  that  are  consistent  w'ith  observed  sediment 
depositional  patterns.  Data  from  the  Mid-Atlantic  Bight  and  the  Northern  Califor¬ 
nia  coEist  are  used  to  define  reasonable  model  input  to  represent  conditions  on  two 
different  types  of  shelves.  In  these  examples,  the  results  show  how  the  intensifica¬ 
tion  of  wave  bottom  velocities  with  decreasing  depth  can  iiitroduce  net  transport 
over  a  region.  The  patterns  of  erosion/deposition  are  shown  to  be  strongly  influ¬ 
enced  by  sediment  stratification  and  moveable  bed  roughness.  Also  predicted  by 
the  applications  is  a  rapid  winnowing  out  of  fine  grain  size  components  when  there 
is  even  a  small  variation  of  bed  grain  size  texture  in  the  along-flow  direction. 
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20, n  bottom  roughness  due  to  grain  roughness  =  ^ 

zqc  wave-enhanced  bottom  roughness,  as  felt  by  (quasi-)  steady  flow  above  the 
wave  boundary  layer 

a  a  function  of  the  maximum  near-bottom  wave-current  velocity,  used  to  calculate 
the  wave-current  shear  velocity 

0  empirical  constant  derived  from  atmospheric  studies  of  stratified  flows  =4.7 
')  empirical  constant  derived  from  atmospheric  studies  of  stratified  flows  =0.74 
70  empirical  reference  concentration  parameter;  in  this  thesis  .0005  <  7o  ^  -005 
6  boundary  layer  thickness  =  O(^)  (cm) 

Sc  thickness  of  current  boundary  layer 
S^  thickness  of  wave  boundary  layer  = 

f  net  erosion  or  deposition  of  sediment  over  a  period  of  time  (cm) 


stability  parameter  measuring  effects  of  sediment  stratification  on  boundary 
layer  fiow 

i7p  sediment  porosity  =  0.6 

Ti  sediment  ripple  height  on  seafloor  (cm) 

Ohtci  direction  of  bedload  transport  with  respect  to  wave  direction 
Ocur  angle  of  mean  current  with  respect  to  the  alongshelf  direction 
0y,  angle  of  wave  direction  with  respect  to  the  alongshelf  direction 
K  von  Karman’s  constamt  =  0.4 
A  wavelength  of  surface  waves  (m) 

A  sediment  ripple  wavelength  on  seafloor  (cm) 
u  fluid  kinematic  viscosity  =  0.013  for  seawater  (^) 

i/t.ew  enhanced  turbulent  eddy  viscosity  in  the  wave  boundary  layer  due  to  waves 
and  current  (^) 

i/t  turbulent  eddy  viscosity  in  unstratified  flow 

Utm  turbulent  eddy  viscosity  of  momentum  in  a  stratified  flow  =  ri^ 

Uf,  turbulent  diffusivity  of  sediment  mass 
p  fluid  density  «  1*0^ 
p,ti  sediment  grain  density 

<7  radian  frequency  of  fluid  motion  =  27r/T;  (sec“^) 


w 

''  V.V.V  •  *  -  *  s'  •  ■  s'  s'  s'  s'  •J'  s'  s'  s'  •  ' 


r  turbulent  shear  stress  in  vertical  direction  =  p{u'w')  =  (em-7ee* ) 

To  boundary  shear  stress 

Toc  enhanced  time-averaged  boundary  shear  stress 

^ociu  maximum  instantaneous  wave-current  boundary  shezir  stress 

Tq  skin  friction  component  of  bottom  shear  stress,  in  this  case  the  maximum  in¬ 
stantaneous  wave-current  boundary  shear  stress 

<f>c  angle  between  direction  of  reference  current  and  direction  of  wave, 

0°  <  <  90° 

Shields  parameter  for  initiation  of  sediment  motion  = 

xl)'  Shields  parameter  based  on  skin  friction  boundary  shear  stress 

xl)0  breakoff  Shields  parameter,  designating  shear  stress  at  which  ripples  begin  to 
be  washed  out  =  1.85J  ® 

critical  Shields  parameter  for  initiation  of  motion  of  a  given  grain  size 
V’m  maximum  Shields  puameter  generated  by  a  given  flow 
w  frequency  of  surface  wave  motion  =  y  C>(10°  —  10“*sec"^) 


Chapter  1 

INTRODUCTION 

The  transport  of  sediment  is  affected  strongly  by  the  presence  of  surface  waves, 
particularly  by  the  large  waves  associated  with  storms  on  the  continental  shelf. 
This  intuitively  obvious  statement  introduces  a  significant  complication  into  the 
field  of  continental  shelf  sedimentology.  Observations  of  flow  and  suspended  sed¬ 
iment  during  storms  are  difficult  to  make;  the  response  of  sediments  to  waves  is 
difficult  to  quantify;  models,  once  made,  are  hard  to  verify.  Marine  sedimentology 
has  therefore  been  limited,  usually,  to  conceptual  and  generally  qualitative  modeb 
of  sediment  transport  and  sedimentation.  These  modeb  are  based  largely  on  ex¬ 
trapolation  from  local  studies,  limited  in  temporal  and  spatial  extent,  sometimes 
augmented  by  laboratory  studies.  Such  modeb  we  of  limited  use  outside  their 
immediate  study  area,  and  there  has  been  a  recognized  need  for  a  predictive  model 
of  erosion,  deposition,  and  depth  of  reworking  which  could  be  applied  in  a  wide 
variety  of  continental  shelf  environments. 

Interest  in  sediment  transport  is  not  limited  solely  to  sedimentologists  and 
stratigraphers.  Shelf  physical  oceanographers  are  concerned  with  the  effect  of  sus¬ 
pended  sediment  on  the  drag  generated  by  the  seafloor.  Coastal  engineers  must  con¬ 
sider  sediment  transport  in  any  continental  shelf  construction  or  pipeline  project. 
Contaminants  often  attach  themselves  to  sediment  grains,  so  the  range  and  rate  of 
sediment  transport  is  of  concern  for  pollution  control.  The  methods  developed  for 


these  other  areas  of  study  can  be  applied  to  aspects  of  the  geological  question,  and 
the  results  may  prove  useful  for  all. 

This  dissertation  presents  an  application  of  recently-developed  boundary  layer 
theory  toward  the  goal  of  a  predictive,  general  sediment  transport  model  for  the 
continental  shelf.  The  aim  here  is  to  demonstrate  how  the  physical  elements  con¬ 
trolling  transport  are  represented  in  the  theoretical  model,  then  test  the  model  to 
see  which  factors  are  most  important  in  calculating  sediment  transport  and  load. 

The  model  used  here  has  grown  out  of  the  theory  to  predict  boundary  shear 
stress  due  to  the  nonlinear  interaction  of  waves  and  currents  developed  by  Grant 
and  Madsen  (1979).  That  model  was  combined  with  a  moveable  bed  roughness 
model  (Grant  and  Madsen,  1982)  by  Scott  Glenn  (1983;  Grant  and  Glenn,  1983), 
who  also  added  stratification  by  suspended  sediments  to  the  formulation.  With 
minor  modifications,  the  theory  and  computer  programs  generated  by  Glenn  were 
used  in  this  work. 

The  explanation  and  application  of  the  model  proceeds  as  follows.  The  re¬ 
mainder  of  the  Introduction  is  devoted  to  an  overview  of  continental  shelf  physical 
systems  and  the  role  of  sediments  in  them,  and  a  brief  history  of  the  development 
of  geological  models  of  shelf  sediment  transport.  In  Chapter  2,  a  physical  descrip¬ 
tion  of  the  elements  of  the  boundary  layer  model  is  presented.  Then,  in  Chapter  3, 
the  theory  is  developed  quantitatively;  that  is,  the  relationships  of  the  parameters 
are  demonstrated  using  the  equations  of  the  theoretical  model  rather  than  descrip¬ 
tions  of  the  physical  relationships.  The  calculation  scheme  used  by  the  model  is 
demonstrated  in  Chapter  4. 

Presentation  of  results  begins  in  Chapter  5  with  a  sensitivity  analysis  which 
demonstrates  the  response  of  sediment  load  zmd  transport  predictions  to  variation 
in  input  parameters.  In  Chapter  6,  the  boundary  layer  model  is  used  to  calculate 
net  erosion  or  deposition  and  reworking  depth  under  wave  and  current  conditions 
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Figure  1.1:  Schematic  diagram  of  the  major  elements  of  the  continental  shelf  physical  sys¬ 
tem,  emphasizing  the  boundary  layers  {Sy,,6£,6,).  Note  the  vertical  extent  (not  necessarily 
to  scale)  of  the  surface  and  bottom  boundary  layers.  Temperature  (cTf)  is  usually  used  in 
physical  oceanographic  studies  to  define  the  limits  of  the  surface  and  bottom  mixed  layers, 
(from  Grant  and  Madsen,  1986). 


representative  of  storms  on  continental  shelves  of  the  Mid- Atlantic  Bight  and  the 
Northern  California  coast.  Chapter  7  states  the  conclusions  of  the  study  and 
suggests  plans  for  improving  and  applying  the  model. 


1.1  Sediments  Within  Shelf  Physical  Systems 

Modelling  sediment  transport  on  the  continental  shelf  is  closely  bound  to  physi¬ 
cal  oceanographic  modelling  of  surface  waves  and  currents  on  the  continental  shelf 
(Figure  1.1).  The  forcing  mechanisms  that  control  processes  in  the  bottom  bound- 
au-y  layer  where  sediment  resuspension  and  transport  take  place  represent  entire 
fields  of  study  in  and  of  themselves.  Shelf  physical  oceanographers  devote  cau'eers 
to  modelling,  theoretically  and  numerically,  the  response  of  the  surface  and  core 
fiows  to  winds,  pressure  gradients,  amd  topography.  The  development  and  disper- 
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sion  of  wind  waves  are  the  subject  of  numerous  theories.  Transport  of  sediment 
by  these  forces  is  influenced  by  processes  representing  still  more  independent  dis¬ 
ciplines:  The  generation  of  bedforms  by  steady  or  oscillatory  flows  is  an  empirical 
field,  with  predictive  models  based  on  laboratory  flume  studies;  benthic  biological 
studies  define  the  effects  of  biota  in  binding  together  and  reworking  the  surface 
sediments,  and  in  generating  roughness  elements  through  their  own  mass  or  the 
sediment  mounds  they  build.  Advances  in  each  of  these  fields  are  necessary  to  a 
complete  shelf  sediment  transport  model. 

The  connections  between  the  processes  represented  by  these  fields  of  study  are 
complex,  involving  feedback  between  the  forcing  and  response  elements  of  the  sys¬ 
tem  (Figure  1.2).  The  bottom  boundary  layer,  for  example,  exists  in  response  to 
wave  and  current  flows  generated  by  winds  and  pressure  gradients.  The  interac¬ 
tion  of  the  bottom  boundary  layer  flow  with  the  bottom,  however,  can  generate 
roughness  elements  which  feed  back  to  influence  the  bottom  boundary  layer  flow. 
Appendix  A  contains  a  more  detailed  diagram  of  the  specific  elements  of  a  complete 
continental  shelf  physical  model  and  a  point-by-point  description  of  the  elements 
and  their  interactions. 

Each  element  can  be  studied,  modelled,  and  tested  separately.  As  each  becomes 
better  understood  on  its  own,  the  influence  of  other  factors  can  be  added.  This 
generates  the  need  for  reliable  models  of  what  would  have  been  called  ‘extraneous 
effects’  in  the  initial  stages  of  study.  In  this  way,  physical  oceanographers  have 
become  more  interested  in  the  effects  of  waves  and,  to  a  more  limited  degree, 
sediments  on  the  frictional  drag  on  currents.  In  this  dissertation,  available  models 
of  boundary  shear  stress,  sediment  entrainment,  bedform  development,  suspended 
load  reference  concentration,  vertical  diffusion  of  mass  and  momentum,  sediment 
concentration  and  velocity  profiles,  and  stratification  by  suspended  sediments  have 
been  combined  into  a  predictive  sediment  model,  forming  a  framework  into  which 
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Figure  1.2:  Master  diagram  ot  the  interacting  systems'^  the  continental  shelf,  extending 
from  the  atmosphere  to  below  the  water*sediment  interface.  Arrows  indicate  forcing  and 
feedback  between  systems.  (Grant  and  Graber,  unpublished) 

future  modelling  developments  can  be  plau:ed. 

1.2  Development  of  Geological  Models  of  Shelf 
Sediment  Transport 

The  geological  utility  of  quantitative  predictions  of  sediment  tramsport  is  three* 
fold.  First,  qualitative  models  that  have  been  formulated  from  observations  of 
input  and  depositional  patterns  can  be  tested.  Second,  the  rates  and  episodicity  of 
transport  can  be  modelled  accurately.  Finally,  models  that  have  been  formulated 
to  explain  stratigraphic  components  of  the  geological  record  for  which  there  are  no 
contemporauy  counterparts  also  can  be  tested. 

Early  geologic  models  presented  the  shelf,  qualitatively,  as  a  profile  of  equi¬ 
librium  between  sea  level,  sediment  input  and  the  action  of  waves  and  currents. 
Dietz  (1963,  1964)  and  Moore  and  Curray  (1964)  traced  the  development  of  the 
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concepts  of  ‘wave  base’,  ‘wave  cut  terraces’  and  ‘profile  of  equilibrium’.  ‘Wave 
base’  had  been  the  most  basic  classification  for  marine  sedimentology,  and  some 
workers  even  suggested  the  depth  of  the  shelf  edge  was  determined  by  wave  base 
(Figure  1.3).  Dietz  and  Moore  and  Curray,  using  the  increasing  data  base  of  seis¬ 
mic  profiles  and  samples  of  the  continental  shelf,  as  well  as  fundamental  concepts 
of  wave  theory  and  initiation  of  sediment  motion,  emphasized  that  the  morpho¬ 
logical  generation  of  the  shape  of  the  continental  shelf  was  controlled  primarily  by 
tectonics  and  sea  level  variation.  Waves  are  important  for  their  capacity  to  erode 
bedrock  in  the  surf  zone,  they  contended,  but  in  water  deeper  than  a  few  meters, 
the  term  ‘wave  base’,  used  geologically,  can  refer  only  to  the  wave’s  capacity  to 
move  unconsolidated  sediment.  No  absolute  depth  could  be  set  for  such  a  defini¬ 
tion,  since  it  would  depend  on  the  wave  height  and  period  and  the  grain  size  and 
areal  distribution. 

The  complexity  of  continental  shelf  processes  has  generated  a  number  of  com¬ 
plementary  approaches.  By  the  1970s,  geologists  such  as  Swift  (1974)  were  char¬ 
acterizing  continental  shelves  based  on  the  wave  and  current  climate  (i.e.  storm- 
or  tide-dominated)  and  the  origin  and  distribution  pattern  of  the  sediments  (i.e. 
palimpsest,  autochthonous).  Patterns  and  rates  of  sediment  transport  were  being 
related  to  observed  environmental  conditions,  and  specific  causes  of  transport  were 
being  isolated  in  detailed  laboratory  and  field  studies  (see,  e.g.,  papers  in  Swift, 
Duane,  and  Pilkey,  1972). 

In  general,  regional  geological  studies  combine  a  variety  of  observations  in  a  par¬ 
ticular  locale  in  attempting  to  determine  the  modern  and  recent  processes  shaping 
the  distribution.  For  example,  Kulm  et  al.(l975)  used  sediment  cores,  bottom 
photos,  current  measurements,  transmissometer  profiles,  and  CTD  (conductivity, 
temperature  and  depth)  profiles  to  construct  the  model  of  sedimentation  illustrated 
in  Figures  1.4  and  1.5.  The  model  indicates  that  sedimentation  is  sensitive  to  a 
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wide  variety  of  variables  in  this  location:  seasonal  waves  ajid  currents,  sediment 
input,  stratification  in  the  water  column,  biological  activity,  topographic  variation 
and  irregularities,  bedforms  and  grain  size.  However,  the  relative  importance  of 
each  variable  can  only  be  estimated  using,  generally,  physical  intuition.  Controlled 
field  experiments  to  examine  the  effect  of  variation  of  each  variable  on  a  regional 
scale  are  prohibitively  expensive  or  impossible. 

Since  the  70s,  advances  have  been  made  in  combining  field,  laboratory,  and 
theoretical  work  from  disparate  fields  such  2is  boundary  layer  dynamics,  physical 
oceanography,  geology,  and  biology.  Progress  on  a  number  of  these  fronts  is  re¬ 
viewed  in  Nowell  (1983).  Interdisciplinary  field  work  such  as  the  experiments  of 
Cacchione  and  Drake  (1982)  and  Butman  (1987a  and  1987b)  attempt  to  explain 
local  transport  depositional  patterns  by  combining  physical  oceanographic  obser¬ 
vations  with  theoretical  predictions. 

Some  of  the  theoretical  tools  for  constructing  a  numerical  sediment  transport 
model  have  been  devised.  The  elements  and  their  state  of  development  at  that  time 
were  described  by  Smith  (1977).  Rachel  (1980,  Rachel  and  Smith,  1986)  applied 
Smith’s  theoretical  concepts  in  a  model  to  explain  patterns  of  deposition  on  the 
Washington  continental  shelf. 

The  remainder  of  this  dissertation  describes  the  approach  taken  and  the  results 
gained  in  applying  the  boundary  layer  model  of  Grant,  Madsen,  and  Glenn  to 
the  sediment  transport  problem.  This  theoretical  model  has  been  described  in 
the  literature  in  a  series  of  papers  and  reports  (Grant  and  Madsen,  1979;  Grant 
and  Madsen,  1982;  Grant  and  Glenn,  1983;  Glenn,  1983;  Glenn  and  Grant,  1987). 
These  will  be  referred  to  collectively  hereafter  as  GMG. 


Chapter  2 


Model  Background:  Concepts 


In  this  section,  the  key  elements  of  continental  shelf  sediment  transport  are  dis¬ 
cussed  in  terms  of  fluid  dynamics,  but  to  the  extent  possible  without  resorting  to 
equations.  The  presentation  follows  the  formulation  developed  in  GMG.  A  the¬ 
oretical  section  follows  (Chapter  3),  elaborating  and  expanding  on  the  concepts 
introduced  here.  All  of  the  elements  of  the  theoretical  treatment  are  introduced 
in  the  present  section,  however,  so  that  the  reader  who  wants  to  avoid  wading 
through  the  mathematical  treatment  can  skip  from  the  end  of  the  present  section 
to  the  Sensitivity  and  Application  sections  (Chapters  5  and  6)  and  understand  the 
physical  signiflcance  of  the  results  presented  therein.  Many  of  the  symbols  used 
in  the  theoretical  sections  aure  introduced  here,  auid  they  provide  the  continuity 
between  the  physical  descriptions  of  this  discussion  and  the  precise  mathematical 
deflnitions  of  the  theoretical  section. 

2.1  The  Continental  Shelf  Boundary  Layer 

This  discussion  of  sediment  transport  by  waves  and  currents  begins  with  an  exam¬ 
ination  of  the  bottom  boundary  liyer.  The  boundary  layer  is  the  region  of  vertical 
velocity  shear  at  a  boundary.  The  iurbulent  bottom  boundary  layer  refers  to  the  re¬ 
gion  above  the  bottom  where  the  flow  is  turbulent  because  of  flow  interaction  with 


the  ocean  floor.  (The  laminar  boundary  layer  is  not  of  concern  here)  Turbulent 
eddies  are  the  mechanism  by  which  the  flow  sets  sediment  in  motion  and  keeps  it 
suspended,  so  the  boundary  layer  is  the  focus  for  study  of  resuspension  and  sub¬ 
sequent  transport  of  sediment.  There  is  no  mechanism  for  sediment  resuspended 
from  the  bottom  to  go  higher  than  the  boundary  layer  since  the  core  flow  (between 
the  surface  mixed  layer  and  the  bottom  boundary  layer)  is  defined  as  inviscid  and 
therefore  free  of  turbulent  shear  stress  (Figure  1.1). 

The  existence  of  a  boundary  layer  hinges  upon  the  presence  of  a  flow  in  the 
core  layer  driven  by  some  outside  force  such  as  tides,  winds,  waves,  or  a  large  scale 
pressure  gradient  (this  is  often  called  the  outer  layer  forcing).  In  shallow  water 
and/or  under  strong  wind  forcing,  the  surface  and  bottom  bound2uy  layers  can 
meet,  so  that  the  bottom  boundary  layer  is  driven  more  directly  by  surface  forcing. 
This  complication  is  not  treated  in  this  model,  however. 

The  velocity  of  the  fluid  in  contact  with  the  bottom  is  by  definition  zero;  this 
follows  from  the  conclusion  of  fluid  dynamicists,  based  on  numerous  experiments 
with  liquids  and  gases,  that  the  micro-layer  of  fluid  in  direct  contact  with  a  solid 
surface  adheres  to  that  surfaK:e.  This  boundary  condition  is  known  as  the  no¬ 
slip  condition  (see  e.g.  Schlichting,  1968),  and  sets  up  the  boundary  layer  region 
of  vertical  velocity  shear  between  the  bottom  auid  the  core.  Bottom  roughness 
triggers  the  formation  of  turbulent  eddies,  which  transport  low-momentum  fluid 
upward  and  higher-momentum  fluid  downward  (Figure  2.1).  The  drag  exerted  on 
the  flow  due  to  these  turbulent  eddies  is  known  as  Reynolds  stress  or  turbulent 
shear  stress. 

2.1.1  Boundary  shear  stress 

The  shear  stress  exerted  on  the  flow  immediately  above  the  seafloor  is  known  as 
the  boundary  shear  stress,  tq.  The  magnitude  of  the  boundary  shear  stress  sets 


Figure  2.1:  Schematic  representation  of  turbulent  eddies  in  the  bottom  boundary  layer 
(after  Yalin  and  Karahan,1979) 


the  level  of  turbulent  energy  for  the  boundary  layer.  Initiation  of  sediment  motion 
and  resuspension  in  the  bottom  few  centimeters  of  the  flow  depend  directly  on  the 
boundary  shear  stress. 

In  the  study  of  sediment  transport,  the  drag  exerted  on  the  seafloor  by  a  steady 
flow  of  fluid  (or,  conversely,  in  physical  oceanography  the  resistance  imposed  by 
the  seafloor  on  a  steady  flow)  is  often  represented  using  a  drag  eoeffieient,  Cp-  The 
drag  coefficient  can  be  used  in  a  quadratic  drag  law  to  relate  the  boundary  shear 
stress  to  a  measured  current  velocity  (u,*/)  at  some  height  above  the  bottom  (see, 
e.g..  Monin  and  Yaglom,  1965): 

To  =  (2.1) 

where  p  is  fluid  density  and  Cp  is  0(10"’);  its  exact  value  is  a  function  of  boundary 
roughness  and  height  of  the  reference  velocity  above  the  seafloor. 

In  an  analogous  fashion.  Grant  and  Madsen  (1979)  use  a  wave-current  friction 
factor  {few)  to  relate  the  instantaneous  boundary  shear  stress  to  the  combined 
wave  and  current  velocity  close  to  the  bottom,  both  of  which  vary  with  the  phase 
of  the  wave.  To  determine  the  drag  on  the  steady  flow,  the  instantaneous  shear 
stress  values  are  time-averaged  over  the  period  of  the  wave  to  yield  the  enhanced 
current  boundary  shear  stress,  Toe-  In  the  few  centimeters  closest  to  the  bottom, 
the  turbulent  intensity  is  aissumed  to  be  governed  by  the  maocimum  instamtaneous 


boundary  shear  stress,  Toc„,,  based  on  the  assumption  that  the  time  scale  for  de¬ 
cay  of  momentum-carrying  eddies  is  long  compared  with  the  wave  period  (Grant 
and  Madsen,  1979;  Smith,  1977).  These  sissumptions  about  the  mean  and  maxi¬ 
mum  boundary  shear  stresses  are  fundamental  to  this  work;  their  ramifications  are 
presented  below.  For  now  it  is  simply  noted  that  different  shear  stresses  control 
turbulence  at  different  heights  above  the  bottom. 

It  is  useful  to  represent  the  bottom  shear  stress  in  terms  of  a  velocity  scale;  this 
is  accomplished  by  defining  a  shear  velocity  u,: 


2.1.2  Boundary  layer  height 

The  thickness  of  the  continental  shelf  boundary  layer  (fi)  is  dependent  on  two 
factors:  1)  the  time  available  for  the  trsuisport  of  turbulent  energy  away  from 
the  bottom,  represented  by  the  inverse  of  the  temporal  frequency  (a,  signifying 
the  radian  frequency)  of  the  flow,  and  2)  the  velocity  scale  of  the  turbulent  eddies, 
represented  by  the  shear  velocity  u«.  The  height  of  the  boundary  layer  is  calculated 
(Grauit  and  Madsen,  1979): 

«  =  0{~)  (2.3) 

where  k  =  0.4  is  von  Karman’s  constant  (Clauser,  1956). 

To  demonstrate  how  flow  frequency  limits  boundary  layer  height,  oscillatory 
flow  will  be  used  as  an  example.  The  sense  of  vorticity  of  turbulent  eddies  depends 
on  flow  direction.  Changes  in  direction  of  the  flow,  therefore,  change  the  sense  of 
vorticity  of  the  eddies  being  generated.  When  flow  direction  changes  in  oscillatory 
flow,  newly-created  eddies  begin  to  cancel  those  with  opposite  vorticity,  inhibiting 
their  spread  upward.  For  example,  the  tidal  current  boundary  layer  has  a  thickness 
on  the  order  of  Sg  =  0{^^)  «  ^  «  O(lOm)  (a  here  is  the  tidal  frequency). 
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Figure  2.2:  Schematic  diagram  of  the  continental  shelf  boundary  layer  illustrating  the 
nested  wave  and  current  boundary  layers  (after  Glenn,  1983). 

The  boundary  layer  due  to  waves  is  much  thinner.  The  maximum  bottom  shear 
stress  associated  with  waves,  re«,,  is  often  up  to  an  order  of  magnitude  larger  than 
the  mean  boundary  shear  stress  r*.  However,  the  wave  frequency,  w,  is  much  higher 
than  that  of  the  mean  flow.  The  wave  boundary,  therefore,  has  a  thickness  6^  = 
**  **  O(lOcm).  The  height  of  the  boundary  layer,  therefore,  is 

governed  by  the  distance  vorticity  is  tr2m8ported  away  from  the  boundary  in  the 
time  limits  imposed  by  the  flow. 

On  the  continental  shelf,  both  surface  gravity  waves  and  quasi-steady  flow  due 
to  wind-driven,  density-driven,  or  tidal  currents  a^e  usually  present.  As  the  cal¬ 
culations  in  the  last  paragraph  suggest,  the  boundary  layers  associated  with  the 
two  processes  are  orders  of  magnitude  different  in  scale,  with  the  wave  boundary 
layer  nested  inside  the  current  boundary  layer  (Figure  2.2).  For  sediment  trans¬ 
port  prediction,  the  wave  boundary  layer  is  a  critical  element  in  the  model,  since 
it  provides  the  boundary  conditions  for  both  velocity  and  concentration  profiles. 


Wave  boundary  layers 


Maximum  near-bottom  velocities  on  the  mid-to-outer  shelf  are  often  similar  for 
waves  and  currents,  although  during  storms  wave  velocities  can  be  much  higher. 
Even  when  the  velocities  are  similar,  the  smaller  length  scale  associated  with  the 
wave  (^u,)  creates  a  maximum  velocity  shear  much  greater  than  the  current  shear, 
generating  a  much  stronger  instantaneous  shear  stress  in  the  wave  boundary  layer 
than  above  it.  The  bottom  shear  stress  in  the  wave  boundary  layer  is  a  nonlinear 
function  of  the  instantaneous  wave  plus  current  velocity,  which  varies  over  the  pe¬ 
riod  of  the  wave  (Smith,  1977;  Grant  and  Madsen,  1979).  As  noted  in  Section  2.1.1, 
the  intensity  of  turbulence  in  the  wave  boundary  layer  is  assumed  to  be  controlled 
by  the  maximum  wave-current  shear  stress.  Therefore,  the  maximum  boundary 
shear  stress,  To^maz  —  ^Oew  is  used  to  define  the  properties  of  the  wave  boundary 
layer,  including  the  shear  velocity  u,ew  and  the  wave  boundary  layer  height; 

(2.4) 

w 

Based  on  the  approach  of  the  predicted  wave  velocity  to  the  free-stream  veloc¬ 
ity,  Grant  and  Madsen  (1979)  suggested  a  value  for  the  scaling  constemt  of  ^=2; 
Trowbridge  and  Madsen  (1984)  suggest  N=1  based  on  agreement  with  higher-order 
models.  Variation  of  this  parameter  is  of  little  importance  for  load  and  transport 
calculations  because  sediment  is  well-mixed  in  the  wave  boundary  layer,  so  concen¬ 
tration  at  is  not  strongly  dependent  on  N.  Unless  otherwise  noted,  the  value 
N  =  2  is  used  in  this  thesis. 

Current  boundary  layers 

Above  the  wave  boundary  layer,  all  turbulent  transport  is  assumed  to  be  associated 
with  the  current.  The  shear  velocity  u»c  represents,  in  this  case,  the  boundairy  shear 
stress  (roe)  felt  by  the  current.  In  the  no-wave  case,  the  drag  felt  by  a  given  steady 


current  is  governed  by  the  bottom  roughness;  more  energetic  eddies  are  produced 
and  more  low-momentum  fluid  is  transported  upward  if  the  bottom  is  rougher.  The 
effect  of  the  wave  is  also  to  increase  turbulent  energy  in  the  region  very  close  to  the 
bottom.  The  wave-generated  turbulence  therefore  can  be  treated  as  a  roughness 
element  acting  on  the  flow  above  the  wave  boundary  layer;  u.g  is  greater  because 
the  current  “feels”  a  rougher  boundary. 

The  frequency  governing  the  boundary  layer  height  associated  with  the  mean 
current  on  the  shelf  is  the  Coriolis  parameter  /  (Grant  and  Madsen,  1986),  defined 
as 

/  =  2( — - )sin<^  (2 

^24hours  ^ 

where  is  the  latitude.  For  mid-latitudes  /  =  The  scale  height  of  the 

current  boundary  layer  is  therefore  defined,  using  the  wave-enhanced  current  shear 

velocity,  as 

4  =  ^  (2.6) 

Like  the  tidal  current  boundary  layer,  the  characteristic  scale  of  the  mean  flow  is 

O(lOm),  with  40-60m  not  uncommon  during  storms. 

2.2  Initiation  of  Sediment  Motion 


Sediment  motion  is  initiated  when  the  shear  stress  felt  by  the  seafloor  is  greater 
than  the  critical  shear  stress  for  moving  sediment.  The  boundary  shear  stress  tq 
defines  the  effects  of  turbulence  on  the  flow  in  the  boundary  layer.  This  stress  can 
be  viewed  as  resulting  from  two  separate  components:  1)  the  interaction  of  the  fluid 
with  the  solid  boundary  and  2)  the  turbulence  generated  due  to  pressure  gradients 
introduced  by  roughness  elements  on  the  bottom.  These  two  components  of  the 
boundary  shear  stress  are  referred  to  as  skin  friction  and  form  drag.  The  medium 
sand  and  smaller  grains  ofprimary  interest  for  suspended  sediment  transport  are 


not  set  in  motion  by  the  pressure  gradients  that  make  up  the  form  drag  component. 
For  the  purposes  of  this  model,  initial  sediment  motion  is  considered  to  result  from 
the  skin  friction  component,  denoted  by  the  symbol  Tq.  Turbulent  transport  of 
mass  and  momentum  in  the  boundary  layer  is  governed  by  the  total  boundary 
shear  stress. 

A  sediment  grain  responds  nearly  instantaneously  to  turbulent  fluctuations. 
Therefore,  critical  shear  stress  for  initiation  of  motion  might  be  expected  to  be 
related  to  the  skin  friction  component  of  the  maximum  boundary  shear  stress, 
Tqcw,  (Section  2.1.1).  In  controlled  laboratory  settings,  initiation  of  motion  and 
bedload  transport  in  oscillatory  flow  were  found  to  be  predicted  quite  successfully 
using  the  maximum  shear  stress  (Madsen  and  Grant,  1976).  In  those  laboratory 
tests,  the  bed  was  flat,  so  the  maximum  shear  stress  was  equal  to  the  skin  friction 
shear  stress.  In  the  calculations  of  this  study,  the  skin  friction  component  of  the 
maximum  combined  boundary  shear  stress  is  used  for  all  initiation  of  motion  and 
bedload  transport  predictions. 

A  commonly  used  empirical  criterion  for  determining  the  critical  shear  stress 
for  initiation  of  motion  of  non-cohesive  sediments  is  the  Shields  parameter,  which 
is  defined: 

(2.7) 

(s  -  l)pgd 

The  numerator  represents  the  force  trying  to  move  the  particle  and  the  denominator 
represents  the  gravitational  force  (per  unit  area)  on  the  particle,  which  resists 
motion  (s  =  ^,p,ed  =  grain  density,  g  =  gravity,  and  d  =  grain  diameter).  When 
the  critical  value  of  Shield’s  parameter  for  the  non-cohesive  grains  in  the  bed  is 
exceeded  by  the  flow,  sediment  is  moved.  The  critical  value  is  designated  0^. 

Critical  values  for  the  Shields  parameter  have  been  determined  empirically  in  a 
series  of  laboratory  experiments,  beginning  with  those  used  to  generate  the  original 
Shields  Diagram  (Shields,  1936).  That  original  diagram  plotted  V’c  versus  the 


boundary  Reynolds  number  iZ*  =  This  can  be  reformulated  to  make  the 

independent  variable  a  function  of  sediment  and  fluid  properties  only  (e.g.  Yalin, 
1972;  Madsen  and  Grant,  1976).  The  modified  Shields  Diagram  plots  tp  versus  5., 
where 


The  original  flume  experiments  on  which  the  diagram  was  based  were  performed 
using  only  large  grain  size  material,  corresponding  to  grain  diameters  of  medium 
sand  or  larger  (it!,  and  5,  >  l).  Later  investigators  did  experiments  with  grains 
as  small  as  d  —  .0016  cm  (medium  silt)  and  used  machine  oil  to  increase  viscosity, 
thereby  extending  the  Shields  Diagram  from  il.  =  1.0  to  i2,  =  0.05  (White,  1970; 
Mantz,  1977;  Yalin  and  Karahan,  1979).  For  this  study,  this  extended  Shields 
Diagram  has  been  adapted  to  the  5,  formulation(Figure  2.3).  The  initiation  of 
motion  criterion  can  be  increased  by  biological  adhesion  of  sediment  grains  (Nowell, 
Jumars,  and  Eckman,  1981;  Grant,  Boyer  and  Sanford,  1982)  or  electrochemical 
cohesion,  due  to  the  presence  of  clays  in  even  small  quantities.  This  is  especially 
true  for  silt-sized  grains;  interpretations  of  sediment  load  and  transport  results 
must  be  made  with  this  uncertainty  in  mind. 

2.3  Sediment  Suspension 

Once  sediments  are  dislodged  from  the  seabed,  they  are  available  for  transport 
upward  by  turbulent  eddies.  The  strength  of  these  eddies  governs  the  mixing 
of  mass  and  momentum  in  the  boundary  layer  and  is  determined  by  the  bound¬ 
ary  sb-  ar  stress,  tq.  Mixing  occurs  because  vertical  eddies  are  transporting  high- 
concentration  fluid  up  and  low-concentration  fluid  down,  so  there  is  a  net  upward 
flux  of  sediment.  This  flux  is  balanced  by  the  tendency  of  the  sediment  to  fall  out 
of  suspension  due  to  gravity. 


S*  =  d  ^(s-1)  gcl/4v 

Figure  2.3:  Modified  Shields  diagram.  The  extension  of  the  initiation  of  motion  criteria 
to  include  recent  laboratory  results  for  fine  grain  is  shown  by  the  solid  line  to  the  left  of 
the  break  in  slope  at  5.  =  1.  The  dashed  line  shows  the  criterion  generated  by  linearly 
extending  the  trend  established  by  Shields’  coarse  grain  data  (modified  from  Madsen  and 
Grant,  1976). 

The  tendency  of  sediment  to  fall  is  measured  by  the  particle  fall  velocity  xuf 
and  is  determined  to  first  order  by  balancing  the  submerged  particle  weight  with 
the  fluid  drag  on  the  particle,  the  weight  being  determined  by  particle  size  and 
density  (Figure  2.4).  Stokes  drag  law  holds  for  grains  for  which  the  nondimensional 
sediment  parzuneter  5.  is  less  them  one  (representing,  for  quartz  grains,  a  diameter 
of  0.012  cm,  or  fine  sand  size),  amd  the  fall  velocity  is  calculated  (Madsen  and 
Grant,  1976): 

- =  -S.  (2.9) 


In  nature,  the  fall  velocity  of  a  particle  can  be  affected  by  flocculation  or  biolog¬ 
ical  aggregation,  so  that  in  situ  measurements  of  the  fall  velocity  are  desirable. 
However,  the  medium  silt  to  medium  samd  grain  sizes  considered  in  this  work  are 
not  usually  subject  to  aggregation.  The  high  wave  boundary  layer  sheair  stresses 
discussed  here  would,  at  any  rate,  likely  disaggregate  mamy  aggregates  (McCave, 
1984).  The  calculated  fall  velocity  is,  therefore,  probably  reasonable,  but  may  be 
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Figure  2.4:  Nondimensionalized  particle  fall  velocity  as  a  function  of  S',  (nondimension- 
alized  grain  diameter;  after  Madsen  and  Grant,  1976). 

a  little  low  at  the  lower  end  of  the  grain  size  spectrum. 

By  using  shear  velocity  «,  to  characterize  the  energy  of  turbulent  eddies,  the 
concentration  profile  near  the  seafloor  cam  be  chauracterized  by  the  ratio  The 
more  energetic  the  turbulence,  the  greater  the  flux  upward;  the  larger  the  particle, 
the  greater  the  flux  down. 

Stratification  by  suspended  sediments 

The  introduction  of  sediment  into  the  flow  from  the  bottom  causes  a  vertical  density 
gradient,  with  the  highest  density  at  the  bottom.  'This  configuration  is  known 
as  stable  stratification  and  tends  to  damp  turbulence.  The  mechanism  by  which 
stratification  damps  turbulence  can  be  understood  in  terms  of  energy  conservation. 
Turbulent  kinetic  energy  in  a  boundary  layer  is  generated  by  the  interaction  of 
mean  velocity  shear  with  the  Reynolds  stress,  as  the  interaction  of  the  flow  with 
the  seafloor  transforms  the  kinetic  energy  of  the  outer  layer  forcing  into  turbulent 
energy.  If  the  fluid  is  stably  stratified,  then  some  of  the  kinetic  energy  is  expended 
in  displacing  the  denser  fluid  upward,  rather  than  in  the  maintenance  of  turbulence. 


The  effect  of  stratification  on  the  turbulent  energy  depends  on  the  density  gra¬ 
dient  and  the  turbulent  energy  of  the  flow.  In  the  extreme  case  where  the  density 
gradient  is  small  and  the  energy  is  large,  then  the  stratification  has  little  effect 
on  the  turbulence,  and  turbulent  mixing  may  destroy  the  stratification.  This  is 
exemplified  in  the  surface-  and  bottom-  “mixed  layers”  on  the  continental  shelf, 
where  temperature  and  salinity  stratification  have  been  destroyed  by  boundary 
layer  turbulence.  At  the  opposite  extreme,  where  the  density  gradient  is  large  and 
turbulent  energy  is  small,  the  turbulence  may  not  be  sufficiently  energetic  to  dis¬ 
place  the  higher  density  fluid  upward,  and  ail  turbulent  mixing  between  the  layers 
may  be  stifled.  This  situation  is  not  unusual  in  estuaries  where  a  turbulent  fresh 
water  layer  may  flow  over  a  turbulent  saline  layer  without  substantial  mixing.  In 
intermediate  cases,  stratification  lessens,  but  does  not  entirely  wipe  out,  turbulent 
mixing. 

The  degree  of  stratification  is  generally  denoted  using  a  Richardson  Number, 
which  is  derived  from  the  turbulent  kinetic  energy  equation.  It  is  the  ratio  of  the 
production  of  turbulent  kinetic  energy  due  to  interaction  of  the  mean  shew  with 
Reynolds  stress  to  the  production  (or  absorption)  of  turbulent  kinetic  energy  due 
to  buoyancy. 

Stable  stratification  of  the  flow  by  sediments  can  affect  both  the  velocity  and 
concentration  profiles.  To  explain  the  eflfects,  the  case  of  a  steady  flow  without 
waves  will  be  examined  first.  Consider  a  flow  in  the  inviscid  core  of  velocity  U, 
over  a  rippled  seafloor  (Figure  2.5a).  The  ripples  interact  with  the  flow  to  generate 
turbulent  eddies,  which  advect  low  velocity  fluid  up  into  the  water  column.  At  some 
height  zi  above  the  bottom,  the  flow  has  been  retarded  to  a  velocity  U  -  Ui.  If 
sediments  stratify  the  flow,  less  turbulent  energy  is  available  to  mix  low-momentum 
fluid  up  from  the  bottom.  This  leads  to  a  larger  vertical  velocity  shear,  so  the 
velocity  at  zi  is  higher  ({/  —  uj,  where  Ui  >  u2)  and  the  drag  due  to  the  boundary 
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Figure  2.5:  Schematic  diagram  showing  the  predicted  effects  of  sediment  stratification 
on  velocity  and  concentration  profiles  in  the  bottom  boundary  layer.  Velocity  profiles  are 
solid  lines;  concentration  profiles  are  dashed,  (a)  Current  profiles  with  no  wave  motion, 
(b)  Wave  and  current  profiles. 

dies  out  closer  to  the  bottom  versus  ^e,n)*  The  latter  effect  means  that  the 
height  of  the  current  boundary  layer  6e  is  smaller  when  there  is  stratification  due 
to  suspended  sediments  in  the  water  column. 

The  effect  on  concentration  profiles  is  analogous  to  that  on  the  velocity.  If 
concentration  profiles  are  calculated  without  taking  stratification  into  account,  the 
concentration  at  zi  reflects  the  balance  between  the  upward  flux  of  sediment  due 
to  turbulence  and  the  downward  flux  due  to  gravity.  Since  stratification  decreases 
the  boundary  layer  turbulent  energy,  a  lairger  concentration  gradient  is  necessary 
to  maintain  the  balance  of  upward  and  downward  fluxes,  and  the  concentration  at 
2i  is  therefore  lower  when  stratification  is  tzd^en  into  account. 

The  presence  of  a  wave  boundary  layer  cam  amplify  the  effect  of  stratifica¬ 
tion.  The  high  energy  associated  with  the  wave-generated  eddies  cause  the  wave 
boundary  layer  to  be  well-mixed.  That  is,  the  sediment  concentration  in  the  wave 
boundary  layer  is  high  throughout.  If  wave  bottom  velocities  are  large  relative  to 
the  near-bottom  current  velocities,  the  turbulent  energy  will  drop  quickly  over  a 
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short  distance  as  the  effects  of  the  waves  die  out. 

The  rapid  decrease  in  turbulent  intensity  translates  to  a  rapid  decrease  in  sed¬ 
iment  concentration,  and  therefore  an  increase  in  stratification  effects.  As  formu¬ 
lated  in  GMG  and  used  here,  with  a  discontinuity  in  eddy  viscosity  at  the  top 
of  the  wave  boundary  layer,  the  fiow  immediately  above  the  wave  boundary  layer 
is  suddenly  able  to  hold  substantially  less  sediment  in  suspension,  causing  strong 
stratification  and  drastically  reducing  the  transfer  of  mass  and  momentum  above 
the  wave  boundary  layer.  Under  strong  waves,  this  wave-enhanced  stratification 
causes  dramatic  increases  in  the  vertical  gradients  of  sediment  concentration  and 
velocity  above  the  wave  boundary  layer  (Figure  2.5b),  and  dramatic  decreases  in 
predicted  load  and  transport.  The  decreases  due  to  stratification  are  discussed  in 
detail  in  the  sensitivity  analysis  (Chapter  5).  The  stratification  effect  due  to  waves 
is  artifically  enhanced  by  telescoping  the  distance  over  which  wave  effects  die  out 
down  to  a  point  at  S^.  The  stratification  effects  discussed  here  are,  therefore,  end 
member  cases  for  the  wave,  current  and  sediment  conditions  presented. 

Determination  of  the  effect  on  the  flow  of  stable  stratification  by  sediments  is 
complicated  by  the  fact  that  the  sediment  concentration  profile  is  itself  a  function 
of  the  turbulent  energy  of  the  flow.  The  feedback  between  the  concentration  profile 
and  the  turbulent  energy  (u,)  means  that  the  calculation  process  must  be 
iterative;  the  method  used  is  described  in  Chapter  4. 

2.4  Sediment  Transport 

The  discussion  so  far  has  been  concerned  with  velocity  and  sediment  concentration 
profiles  in  the  boundary  layer  at  a  point  for  a  single  set  of  conditions.  These 
quantities,  by  themselves,  can  be  used  to  calculate  the  depth  of  reworking  smd  the 
net  flux  past  the  point  over  a  given  period  of  time.  To  predict  the  net  sediment 
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transport,  that  is,  erosion  or  deposition,  for  a  finite  region  on  the  continental 
shelf  over  a  finite  time,  a  model  must  account  for  spatial  variation  of  one  or  more 
environmental  parameters.  The  parameter  variation  could  take  a  number  of  forms, 
such  as  a  variation  in  wave  height  due  to  water  depth  change;  a  variation  in  bottom 
roughness  due  to  biological  activity;  a  variation  in  grain  size;  or  a  vwiation  in 
current  speed  or  direction.  In  cases  where  there  is  no  horizontal  change,  the  input 
of  sediment  to  an  area  would  equal  the  output  from  it,  and  there  would  be  neither 
net  erosion  nor  deposition. 

For  purposes  of  this  exercise  we  let  a  region  of  any  continental  shelf  be  repre¬ 
sented  by  a  square  (Figure  2.6),  defined  by  four  points  at  its  corners,  over  which 
sediment,  wave,  and/or  current  parameters  are  allowed  to  vau’y.  The  box  defined 
by  the  grid  square  and  the  boundary  layer  directly  above  it  comprises  the  field  of 
interest.  The  horizontal  extent  chosen  for  a  grid  square  is  dependent  on  the  spatial 
homogeneity  of  conditions  in  the  region  of  interest  and  the  sensitivity  of  the  pre¬ 
dictions  to  changes  in  the  pMameters.  For  example,  changes  in  water  depth  change 
the  wave  parameters  at  the  bottom  so  that  a  length  scale  might  be  determined  by 
the  bottom  slope. 

The  Grant,  Madsen,  and  Glenn  model  is  a  steady  state  model,  reflecting  the 
assumption  that  flow  and  sediment  profiles  adjust  to  changes  in  forcing  (waves  and 
currents)  quickly  relative  to  the  duration  of  the  flow  conditions.  This  assumption  is 
reeisonable  for  the  near-bottom  layer  under  all  circumstances,  as  the  eddy  diffusivity 
associated  with  the  current  will  mix  mass  and  momentum  to  a  height  of  about 
10  meters  in  a  matter  of  minutes.  Storm  conditions  typically  last  for  one  to  three 
days,  and  stationarity  in  storm  conditions  can  be  assumed  for  a  matter  of  hours. 
The  question  of  time  variance  is  left  to  future  models. 

At  the  corners  of  the  box,  the  sediment  load  and  sediment  flux  must  be  cal¬ 
culated.  To  do  so  requires  the  calculation  of  the  velocity  profile  and  the  sediment 

45 


x,y  -f  Ay 


I  +  Ax,y  -h  Ay 


Ay  9«,,  • 


*  +  Ax.y 


Figure  2.6:  Sample  horizontal  grid  square  for  calculating  net  transport  on  an  area  of  the 
continental  shelf 

concentration  profile  in  the  bottom  boundary  layer.  This  can  be  accomplished 
using  GMG;  the  method  of  calculation  is  described  in  the  next  two  sections.  The 
average  erosion  or  deposition  for  the  grid  square  is  determined  by  vector  averag¬ 
ing  the  transport  predictions  at  each  corner  and  dividing  the  net  trauisport  by  the 
area  of  the  square.  The  average  depth  of  reworking  is  calculated  by  averaging 
the  sediment  load  predictions  at  the  four  corners  and  dividing  by  the  bed  grain 
concentration,  here  taken  to  be  0.6  (Yalin,  1972). 


Chapter  3 


Model  background:  Theory 


There  are  three  quantities  of  interest  in  this  study  of  sediment  transport:  1)  the 
net  erosion  or  deposition  of  sediment  over  a  period  of  time;  2)  the  sediment  load, 
or  total  volume  of  sediment  stored  in  the  water  column,  from  which  the  reworking 
depth  can  be  calculated;  and  3)sediment  transport,  or  rate  at  which  the  sediment 
moves  across  the  seafloor.  These  three  quantities  are  related  through  the  equation 
of  conservation  of  sediment  mass: 

=  (3.1) 

The  net  erosion  or  deposition  is  represented  in  the  left-hand  side  of  this  equation, 
where  is  the  change  in  level  of  the  seafloor  (?)  over  time  (t),  and  qp  =  sediment 
porosity.  The  two  terms  on  the  right-hand  side  represent  1)  changes  in  the  load 
with  time  and  2)  the  horizontal  transport  divergence,  where: 

Load{^)  =  f  C(z)dz  (3.2) 

cm*  J 

q,  =  Transport{cm^ ! cmj sec)  =  j  C{z)u{z)dz  (3.3) 

The  symbol  is  the  mathematical  vector  operator  “del” ,  applied  in  the  horizontal 
plane,  defined  byVH  =  i^+j^.  The  dot  product  of  Vh  and  the  transport  vector 
is  designated  the  ‘horizontal  divergence’  and  represents  spatial  change  in  transport 


volume.  It  is  defined: 


where  and  q,  y  are  the  x  and  y  components  of  the  transport  vector  and  represent 
the  net  rate  of  sediment  flux  into  a  control  square  of  unit  size.  The  quantities 
that  must  be  calculated  are  the  velocity  and  sediment  concentration  profiles.  In 
this  study,  these  are  generated  using  GMG.  The  theoretical  foundations  for  these 
calculations  are  discussed  briefly  here;  they  can  be  studied  in  more  detail  in  GMG. 

3.1  Near-bottom  Boundary  Layer:  No  Suspended 
Sediments 

3.1.1  Governing  Equations 

The  horizontal  flow  in  the  boundary  layer  is  governed  by  the  Reynolds  averaged 
equations  of  momentum: 

du  \dp  d  ,  ,  . 

dv  idp  d  ,  ,  . 

where  (x,y)  and  (u,v)  are  horizontal  directions  smd  velocities,  respectively;  and 
w  is  the  vertical  velocity,  positive  upward.  Primed  variables  represent  turbulent 
velocity  components;  p  is  pressure.  The  terms  {u'w')  =  ^  and  (v'u;')  =  ^  represent 
the  Reynolds  stresses  defined  in  Section  2.1;  they  result  from  Reynolds  averaging 
the  turbulent  velocity  components  over  a  time  that  is  long  compared  to  the  time 
scale  of  the  turbulence.  The  terms  in  these  equations  show  that  the  velocity  profile 
is  determined  by  the  Coriolis  force,  horizontal  pressure  gradients,  and  Reynolds 
stresses.  The  assumptions  underlying  this  form  of  the  equations  are  that  the  flow, 
the  boundary,  and  the  roughness  elements  are  all  horizontally  homogeneous,  so 
that  variations  in  x  and  y  directions  can  be  neglected. 
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The  velocity  and  pressure  terms  can  be  divided  into  wave  and  current  contrr 
butions,  assuming  wave  motion  is  in  the  positive  x  direction: 


—  ti-g  “f"  ttti; 

(3.7) 

p  =  Pc  +  Pw 

(3.8) 

Wave  Governing  Equation 


After  substituting  these  into  the  momentum  equation  (Equation  3.5),  the  terms 
governing  the  mean  current  motion  can  be  separated  from  those  governing  the  wave 
motion.  The  equation  governing  the  wave  motion  becomes: 

dutu  1  dpw  9 


{u'w') 


at 


p  dx  dz 


(3.9) 


Current  Governing  Equations 


The  current  is  assumed  to  be  steady  over  the  time  scales  of  interest,  so  that  varia¬ 
tions  in  time  are  neglected.  Assuming  also  that  the  flow  is  driven  by  Coriolis  forcing 
in  the  inviscid  core,  the  pressure  terms  can  be  replaced  using  the  geostrophic  equa¬ 
tions: 


^  ~  pdx 


_  1 5Pc 

/**oo  —  ^ 

P  oy 


(3.10) 

(3.11) 


The  equations  governing  the  mean  flow  become: 


-/(v  -  Voo)  = 


(3.12) 


/(Uc  -  Uoo)  = 


(3.13) 


Coriolis  effects  can  be  neglected  close  to  the  bottom  where  the  assumption  of  a 
constant  stress  layer  is  made;  for  simplicity  in  the  discussion,  waves  and  currents 


are  assumed  to  be  collinear.  The  momentum  equations  governing  the  mean  flow  in 
the  near-bottom  layer  (Equations  3.12  and  3.13)  then  simplify  to  a  single  equation: 

^(uV)=0  (3.14) 

oz 

The  assumption  of  wave-current  collinearity  is  not  necessary  for  the  calculation, 
and  the  model  is  equipped  to  handle  differences  in  wave  and  current  direction. 

The  constant  stress  region  is  referred  to  throughout  this  work  as  the  near¬ 
bottom  region,  distinguishing  it  from  what  is  referred  to  as  the  outer  Ekman  region, 
where  shear  stress  drops  and  the  current  direction  veers  due  to  Coriolis  accelera¬ 
tions.  The  largest  volume  of  shelf  sediment  transport  is  concentrated  in  the  bottom 
few  meters  of  the  water  column  under  most  circumstances.  The  near-bottom  re¬ 
gion  is  the  focus  of  the  GMG  model,  and  the  determination  of  its  vertical  limits 
are  based  on  the  current  shear  stress  as  discussed  in  Section  3.2.4. 

In  the  outer  Ekman  layer,  the  effects  of  changes  in  direction  must  be  taken  into 
account,  and  the  full  Ekman  layer  equations  (Equations  3.12  and  3.13)  must  be 
solved  using  a  numerical  solution.  This  is  not  done  in  the  present  model;  instead, 
the  model  calculates  a  rough  estimate  of  the  total  load  and  transport  in  the  outer 
Ekman  layer  based  on  the  values  of  sediment  concentration  and  eddy  viscosity  at 
the  outer  edge  of  the  constant  stress  layer.  The  method  used  to  make  the  estimates 
is  described  in  Section  3.2.4.  In  cases  where  the  load  and  transport  in  the  outer 
Ekman  layer  are  the  same  as  or  greater  than  the  near-bottom  transport,  the  volume 
and  direction  of  predictions  of  total  transport  must  be  treated  with  caution. 

3.1.2  Turbulent  Closure 

To  solve  the  simplifled  wave  and  current  momentum  equations  (Equations  3.9 
and  3.14)  for  the  near-bottom  velocity  profiles,  the  Reynolds  stresses  must  be 
made  tractable.  By  analogy  with  molecular  frictional  shear  stress,  the  Reynolds 


stress  can  be  modelled  using  an  eddy  viscosity,  Uti 


-  {u'w')  -  —  =  (3.15) 

p  oz 

where  eddy  viscosity  is  a  parameterization  of  the  intensity  of  the  turbulence  based 
on  the  flow  velocity  and  the  length  scale  of  the  eddies  (cf.  the  fluid  viscosity  v 
which  is  a  property  of  the  fluid).  The  turbulent  eddy  viscosity  model  for  the  near¬ 
bottom  layer  chosen  by  GMG  and  used  in  this  work  is  time  invariant  and  linear 
with  distance  from  the  bottom: 

Ut  —  Ku,z  (3.16) 

Distance  from  the  bottom,  z,  is  included  because,  close  to  the  bottom,  the  vertical 
length  scale  of  the  turbulent  fluctuations  can  be  no  larger  than  the  distance  from 
the  bottom.  The  shear  velocity,  u,,  is  a  measure  of  the  shear  stress  on  the  bottom, 
as  defined  in  Equation  2.2.  Close  to  the  boundary,  the  turbulent  shear  stress 
has  been  observed  to  be  approximately  constant  and  equal  to  the  boundary  shear 
stress.  To-  The  region  where  this  approximation  holds  is  referred  to  as  the  constant 
stress  region.  Substituting  the  eddy  viscosity  representation  of  the  Reynolds  stress 
into  the  near-bottom  momentum  equation  governing  the  mean  flow  (Equation  3.14) 
and  solving  it  with  the  boundary  condition  that  the  bottom  shear  stress  is  equal 
to  tq  leads  to  the  constant  stress  layer  definition: 


I't 


du 

dz 


(3.17) 


Recall  from  the  discussion  of  boundary  layers  (Section  2.1.1)  that  the  char¬ 
acteristic  shear  velocity  is  different  for  the  wave  boundary  layer  and  the  current 
boundary  layer.  In  the  wave  boundary  layer,  turbulent  energy  is  assumed  to  reflect 
the  maximum  combined  wave  and  current  boundary  shear  stress: 


where  7  represents  a  functional  relationship.  Above  the  wave  boundary  layer,  tur¬ 
bulent  eddies  are  associated  with  the  time-averaged  current  shear  stress  to^  alone. 
This  means  that  eddy  viscosity  is  defined  differently  inside  the  wave  boundary  layer 
and  above  it: 

Ut  —  KU,cZ  z  >  6^  (3.20) 

The  limits  of  the  constant  stress  layer  and  the  eddy  viscosity  model  for  the  outer 
region  of  the  current  boundary  layer  are  discussed  in  more  detail  in  Section  3.2.4. 


3.1.3  Boundary  Shear  Stress  Calculation 


The  characteristic  boundary  shear  stresses  and  shear  velocities  are  calculated  from 
the  instantaneous  boundary  shear  stress.  The  instantaneous  boundary  shear  stress 
is  defined  in  GMG  using  a  quadratic  drag  law: 


-*  2\r  ^  ^ 

’’o  =  rp/cu.(u  +  V )[-  2V^’7~nr~2V^ 


(3.21) 


where  u,  v  are  the  x,y  components  of  a  combined  wave  and  current  reference  veloc¬ 
ity  vector  close  to  the  bottom.  In  this  discussion,  however,  the  wave  and  current 
are  tissumed  to  be  collinear  in  the  i-direction.  fcv  is  the  combined  wave  and  cur¬ 
rent  friction  factor,  analogous  to  Jonsson’s  (1966)  wave  friction  factor  for  pure 
oscillatory  motion.  The  characteristic  shear  stress  in  the  wave  boundary  layer 
{‘’’ocw  =  P  discussed  in  Section  2.1,  is  defined  as  the  maximum  value  of 

Equation  3.21.  For  the  current  boundary  layer,  tqc  (=  puJJ  is  calculated  by  time 
averaging  Equaion  3.21.  The  solutions  for  the  shear  velocities  are: 


where  a  and  V2  are  functions  of  the  mjiximum  and  time-averaged  velocities,  re¬ 
spectively,  in  the  wave  boundary  layer.  «(,  is  the  maximum  wave  bottom  velocity, 
calculated  using  linear  wave  theory: 


H  LJ 
2  sinhA:h 


(3.24) 


where  H  =  trough-to-crest  wave  height,  k  —  is  wave  number,  and  A  is  wave¬ 
length. 

Ua  is  a  representation  of  the  velocity  of  the  mean  flow  in  the  wave  boundary 
layer,  so  ^  is  a  representation  of  the  relative  strength  of  the  mean  versus  the 
maximum  oscillatory  flow  in  the  wave  boundary  layer.  The  value  of  fc^,  is  calculated 
using  these  definitions  and  the  wave  velocity  profile.  The  solutions  for  the  friction 
factor  and  the  wave  velocity  profile  are  found  in  Appendix  B.  The  functional 
depencence  of  fcui  is: 

feiu  ~  ^ (y^bi^ref  f  ■^b'}  (3.25) 


where  uj  is  the  maximum  wave  bottom  velocity,  4>c  is  the  angle  between  wave  and 
current  directions,  kb  is  the  bottom  roughness,  and  Ab  is  the  bottom  excursion 
amplitude  for  the  wave,  defined: 


u 


(3.26) 


3.1.4  Velocity  Profile  Solutions 


With  the  definitions  of  shear  velocity  and  eddy  viscosity  given  in  Sections  3.1.2  and  3.1.3 
and  the  bottom  roughness  discussed  in  Section  3.3,  solutions  for  the  steady  and 
oscillatory  velocity  profiles  can  be  determined  from  Equations  3.14  and  3.9. 

As  discussed  in  Section  2.1.2,  the  effects  of  the  boundary  on  wave  motion  are 
limited  to  a  thin  layer  close  to  the  bottom,  the  wave  boundary  layer.  The  shear 
velocity  in  that  region  is  a  function  of  the  combined  wave-current  velocity  near  the 
bottom.  Above  the  wave  boundary  layer,  the  turbulent  momentum  flux  associated 


with  the  wave  is  small.  The  solution  for  fluid  motion  due  to  the  wave  is  simply 
the  linear  wave  solution  for  inviscid  flow,  predicting  sinusoidal  motion  decreasing 
in  amplitude  with  depth  (e.g.,  LeBlond  and  Mysak,  1978,  sec.  11).  For  waves  that 
generate  significant  velocities  at  the  seafloor,  (see  Section  5.1),  the  wave  motion  at 
the  top  of  the  wave  boundary  layer  is: 

u,u  =  ujsincjt  (3.27) 


The  solution  for  the  momentum  equation  for  the  waves  inside  the  wave  bound¬ 
ary  layer  (Equation  3.9)  is  not  explicitly  of  interest  for  the  present  problem,  though 
it  is  necessary  for  the  calculation  of  the  boundary  shear  stress.  The  solution  is  given 
in  Appendix  B. 

Because  of  the  difference  in  turbulent  mixing  in  the  wave  boundary  layer  and 
the  area  above  it,  the  solution  for  the  mean  (current)  velocity  is  different  inside  and 
above  the  wave  boundary  layer.  The  profiles  are  calculated  from  Equation  3.17, 
using  the  different  values  of  eddy  viscosity  inside  and  outside  the  wave  boundzu’y 
layer  as  defined  in  Equations  3.19  and  3.20: 


u(2)  =  -U,e - In —  z  <  0^ 

K  ZQ 


(3.28) 


u(z)  =  — u,cln  — 

K  Zoe 


z>  6^  (3.29) 

The  influence  of  the  enhanced  boundary  shear  stress  on  the  mean  flow  inside  the 
wave  boundary  layer  is  demonstrated  by  the  term  in  Equation  3.28.  Since  the 
combined  wave-current  shear  velocity  u.cw  is  by  definition  larger  than  u,e,  the  mean 
flow  is  reduced  in  the  wave  boundary  layer  relative  to  what  it  would  be  if  turbulent 
intensity  were  governed  entirely  by  the  mean  shezu’  stress.  The  lower  velocities 
reflect  the  increased  drag  on  the  flow  caused  by  wave-generated  turbulence. 

The  value  of  zq,  where  the  calculated  flow  velocity  goes  to  zero,  is  based  on  the 
bottom  roughness.  The  method  used  to  calculate  it  is  discussed  in  Section  3.3.  The 


enhanced  bottom  roughness,  zoc,  represents  the  total  roughness  felt  by  the  mean 
flow,  including  the  effects  of  the  increased  eddy  viscosity  due  to  the  presence  of  the 
wave.  The  value  of  2oc  is  determined  by  matching  the  calculated  velocities  for  each 
region  (Equations  3.28  and  3.29)  at  the  top  of  the  wave  boundary  layer,  6^: 


/  ^10  \1-  -***« 
20c  =  Zo{—) 

ZO 


(3.30) 


Since  >  2o  by  definition,  zqc  is  likewise  greater  than  zq.  This  reflects  the  in¬ 
creased  transport  by  wave-generated  turbulence  of  low-momentum  fluid  up  from 
the  seafloor  to  the  top  of  the  wave  boundary  layer. 

3.2  Near-Bottom  Boundary  Layer:  With  Sus¬ 
pended  Sediments 

The  last  section  introduced  the  concepts  underlying  the  study  of  combined  wave 
boundary  layers.  The  neutral  velocity  profile  solutions  (Equations  3.28  and  3.29) 
are  sufficient  for  the  study  of  flow  in  the  boundary  layer  when  the  bottom  is  bedrock 
or  sediment  grains  so  large  that  they  will  not  go  into  suspension.  The  focus  of 
this  study,  however,  is  the  effect  of  waves  on  sediment  transport,  rather  than  the 
effect  of  waves  on  the  drag  encountered  by  the  current.  Although  the  sediment 
concentration  profile  is  a  function  of  the  turbulent  mixing  capacity  of  the  flow, 
as  represented  by  u,,  the  results  from  the  neutral  boundary  layer  model  cannot 
simply  be  applied  to  the  calculation  of  a  concentration  profile.  This  is  because 
stratification  by  sediments,  as  discussed  in  Section  3.2,  can  damp  the  turbulent 
transfer  of  mass  and  momentum  up  from  the  seafloor,  affecting  the  values  of  u,,,. 
To  calculate  a  velocity  profile  that  takes  stratification  into  account,  the  sediment 
concentration  profile  also  must  be  determined.  Since  the  concentration  profile  is 
itself  dependent  on  the  stratification  correction,  the  procedure  has  to  be  iterative. 
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Governing  Equations 
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The  simplified,  Reynolds-averaged  conservation  of  momentum  equations  for  waves 
and  currents  (Equations  3.9  and  3.14)  still  govern  the  flow.  Again,  for  ease  of 
discussion,  co-directional  waves  and  currents  in  the  x-direction  in  the  near-bottom 
layer  are  assumed.  The  mean  distribution  of  sediment  in  the  water  column  is 
governed  by  the  conservation  of  mass  equation  for  sediment  (also  known  as  the 
sediment  continuity  equation); 


(3.31) 


where  C  is  the  time-averaged  volumetric  sediment  concentration  amd  {C'w')  rep¬ 
resents  the  Reynolds  averaged  turbulent  fluctuation  of  sediment  flux.  Inside  the 
wave  boundary  layer,  there  is  a  fluctuating  instantaneous  concentration  associated 
with  the  oscillatory  wave  velocity,  which  is  taken  into  account  in  calculating  the 
reference  concentration. 

Analogous  to  the  eddy  viscosity  representation  for  turbulent  stress,  turbulent 
mixing  of  sediment  likewise  can  be  modelled  using  an  eddy  diffusivity,  so  that 


{C'w')  = 


(3.32) 


It  is  most  often  assumed  that  the  eddy  diffusivity  and  eddy  viscosity  can  be  repre¬ 
sented  in  boundary  layer  flows  using  similar  forms,  and  t/t,  can  be  written  2is 


Ut  K 

l/ft  =  —  =  — u*z 
1  '1 


(3.33) 


where  'i  is  an  empirical  parameter,  assumed  to  be  0.74,  based  on  Businger  and 
Arya  (1974).  The  GMG  model  2issumes  steady  state  conditions,  with  a  balance  be¬ 
tween  upward  turbulent  transport  and  gravitational  transport  down.  Thus  Equa¬ 


tion  3.31  can  be  written 


rvfC  +  =  0 


(3.34) 


This  is  satisfied  by 


C{z)  =  C{zq]{—)  zo  <  z  <  6^  (3.35) 

C{z)  =  z>6^  (3.36) 

Oy, 

where  C{zo)  is  a  reference  sediment  concentration,  discussed  in  Section  3.2.3,  and 
C{6J)  is  the  concentration  at  the  top  of  the  wave  boundary  layer,  determined  from 
Equation  3.35. 

3.2.1  Stratification  Effects 

The  effect  of  stratification  on  velocity  and  concentration  profiles  is  treated  using 
the  eddy  viscosity  models.  By  analogy  with  atmospheric  models,  the  eddy  viscosity 
for  momentum  in  the  presence  of  stratification  is  defined: 

Utm  =  ^  (3.37) 

<Pm 

and  the  eddy  diffusivity  for  mass  is  defined: 

Ut.  =  ~  (3.38) 

9$ 

where  Ut  is  the  neutral  eddy  viscosity,  and  and  <^,  are  the  nondimensional 
velocity  and  concentration  gradients,  respectively,  defined: 

<l>m  =  l  +  PU  (3.39) 

0,  =  7  +  /??,  (3.40) 

where  /?  is  another  empirical  parameter,  with  value  4.7  (Businger  and  Arya,  1974). 

is  a  stability  parameter,  derived  from  the  Richardson  number,  which  measures 
the  effect  of  stratification  on  the  flow.  The  stability  parameter  in  the  near-bottom 
flow  (where  the  constant  stress  assumption  holds)  is  defined: 


where  L  is  the  Monin-Obukov  Length,  a  parameterization  of  the  ratio  of  the  pro¬ 
duction  of  turbulent  kinetic  energy  by  mean  shear  to  its  dissipation  by  stratification 
(or  production  by  buoyancy;  Long,  1981): 


/c5(p'ti;') 


(3.42) 


Near  the  bed,  all  density  fluctuations  are  assumed  to  be  caused  by  the  vertical 
sediment  concentration  gradient,  and  so  the  turbulent  fluctuation  of  mass  is  mod¬ 
elled  using  an  eddy  diffusivity  and  the  sediment  concentration  profile  (Long, 1981): 

{p'w')  =  p{s  -  l)(C'u;')  =  -p[s  -  l)uu^ 

Before  installing  this  expression  into  the  stability  parameter,  the  differential  is 
disposed  of  by  substituting  in  from  the  steady-state  equation  of  conservation  of 
mass  (Equation  3.34),  so  that 


2  vtg{s  -\)wfC{z) 
L~  u* 


(3.43) 


The  value  of  the  stability  parameter  varies  from  near  zero,  when  stratification 
effects  are  very  small,  and  is  of  order  one  when  stratification  effects  are  large. 

3.2.2  Velocity  and  Concentration  Profile  Solutions 

Using  the  definitions  of  eddy  viscosity  altered  by  stratification  effects  (Equations  3.37 
and  3.38),  the  solutions  for  the  steady  current  can  be  determined.  The  solution  for 
the  current  above  the  wave  boundary  layer  becomes: 


u(2)  = -u,c[ln  —  + /?  /  yds]  z>6^ 

K  ZQc  J(^  Li 


(3.44) 


The  effects  of  stratification  are  shown  in  the  second  term  on  the  right  hsmd  side 
of  Equation  3.44.  The  equation  demonstrates  the  effect  mentioned  in  Section  2.3 
that,  for  a  given  boundary  shear  stress  u,,  a  higher  velocity  is  measured  at  a  height 
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z  if  there  is  stratification  by  sediments.  Conversely,  for  a  measured  velocity  in  the 
field  at  a  height  z  above  the  bottom,  less  shear  stress  is  felt  on  the  bottom  if  there 
is  sediment  stratification. 

Because  the  high  energy  of  the  eddies  in  the  wave  boundary  layer  is  expected 
to  keep  that  region  well  mixed,  the  stratification  correction  is  not  included  in 
the  calculation  of  the  velocity  below  and  the  solution  for  the  current  profile  is 
unchanged  from  Equation  3.19.  The  validity  of  this  assumption  for  the  GMG  model 
is  demonstrated  in  an  analysis  by  Glenn  (1983),  where  the  maximum  expected 
correction  to  velocity  due  to  stratification  in  this  model  is  shown  to  be  an  order  of 
magnitude  too  small  to  affect  turbulent  fluxes. 

The  concentration  profile  inside  the  wave  boundwy  layer  is  likewise  unchanged 
from  the  neutral  case  (Equation  3.35).  The  solution  for  the  concentration  profile 
above  the  wave  boundary  layer  is; 


C(z]  =  exp{-^J-  J  jdz}  z>6, 


(3.45) 


Note  that  the  exponential  decay  of  the  concentration  increases  with  an  increase 
in  the  stability  parameter.  The  integral  in  the  exponential  term  must  be  solved 
numerically. 

The  calculation  of  the  concentration  above  the  wave  boundary  layer  (Equa¬ 
tion  3.45)  depends  directly  on  the  calculation  of  the  concentration  at  the  top  of 
the  boundary  layer,  C(d,u),  as  calculated  using  Equation  3.35.  In  addition  to  the 
decay  in  concentration  due  to  the  balzmce  of  the  fall  velocity  with  the  turbulent 
mixing  as  in  the  wave  boundary  layer,  the  exponential  decay  term  reflects  the 
reduction  in  turbulence  due  to  stratification.  If  the  fall  velocity  is  small  relative 
to  the  shear  velocity  {wj  «  u,),  the  concentration  decays  relatively  slowly  and 
stratification  has  little  effect  on  the  concentration  and  velocity  profiles.  For  rela¬ 
tively  large  fall  velocities  (w/  >  u,),  the  concentration  decays  rapidly  regardless  of 


stratification,  so  stratification  again  has  relatively  little  effect.  The  exception  to 
this  rule  occurs  when  there  is  a  strong  wave  and  relatively  weak  current  so  that 
u.cui  >>  u.c-  In  this  case,  the  concentration  will  drop  off  rapidly  above  the  wave 
boundary  layer  regardless  of  grain  size,  creating  a  strong  stratification  effect.  This 
tendency  is  investigated  in  the  Sensitivity  section. 

The  discussion  has  addressed  only  the  case  of  a  single  grain  size  in  the  seabed. 
This  is  for  ease  of  presentation;  the  GMG  model  can  accomodate  up  to  ten  grain 
size  clcLsses.  The  procedure  is  the  same,  with  grain  diameter  and  bed  concentration 
being  specified  for  each  class.  The  stratification  correction  is  calculated  using  the 
sum  of  the  concentrations  for  all  categories. 

3.2.3  Reference  Sediment  Concentration 

The  solution  for  the  sediment  concentration  profile  requires  that  a  reference  con¬ 
centration  be  specified  (C(«o)  in  Equation  3.35).  The  reference  concentration  is 
calculated  in  the  model  using  the  form  suggested  by  Smith  and  McLean(l977): 


CM  = 


(3.46) 


where  Cjed  is  the  bed  concentration  of  the  grain.  70  is  an  empirical  reference 
concentration  parcimeter,  auid  estimates  of  its  value  range  from  0.0001  to  0.005 
(Glenn,  1983;  Wiberg  and  Smith,  1983).  S  is  the  normalized  excess  skin  friction 


5  =  1^ 

Tc 


-  1 


The  primes  refer  to  the  skin  friction  component  of  the  total  boundary  sheair  stress. 
The  skin  friction  component  is  determined  by  calculating  the  combined  wave- 
current  friction  factor  few  using  the  dominant  grain  size  as  the  bottom  roughness 
scale  rather  than  the  physical  boundary  roughness  (Equation  B.l).  To  accomodate 
the  presence  of  waves  in  this  model,  the  instantaneous  normalized  excess  shear 
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stress  is  used  to  calculate  instantaneous  reference  concentrations  which  then  are 
averaged  over  a  wave  period  to  find  the  mean  reference  concentration. 

The  reference  concentration  is  directly  dependent  on  the  critical  shear  stress,  as 
determined  by  the  Shields  parameter.  However,  the  critical  Shields  parameter  is  an 
empirical  value  based  on  laboratory  flume  experiments  on  single  grain-size  sands. 
Mixed  grain  sizes  and  biological  binding  or  mixing,  as  found  in  field  situations, 
may  affect  the  critical  shear  stress.  The  effects  of  these  uncertainties  in  transport 
and  load  predictions  are  addressed  in  the  Sensitivity  section. 

3.2.4  Sediment  in  Outer  Ekman  Layer 

The  Grant-Madsen-Glenn  near-bottom  model  calculates  velocity  and  concentra¬ 
tion  profiles  only  for  the  constant  stress  region  of  the  boundary  layer  flow.  For 
estimates  of  sediment  load  and  transport,  it  is  necessary  to  specify  a  height  limit 
over  which  the  concentration  profile  is  to  be  integrated,  so  a  vertical  limit  to  the 
constant  stress  layer  must  be  specified.  Turbulent  mixing  extends  beyond  the  con¬ 
stant  stress  region,  however,  so  it  is  possible  that  signific2uat  proportions  of  the  total 
sediment  load  are  located  above  the  constant  stress  region,  in  what  will  be  referred 
to  as  the  outer  Ekman  layer.  To  estimate  the  magnitude  amd  relative  importance  of 
the  outer  Ekman  layer  transport,  the  GMG  model  was  altered  as  discussed  below. 
In  the  near-bottom  layer,  eddy  size  scales  with  distance  from  the  bottom,  so  z  is 
an  appropriate  length  scale  for  the  eddy  viscosity.  This  region  can  be  modelled  ais 
a  constant  stress  layer  with  a  linearly  increasing  eddy  viscosity;  numerous  experi¬ 
ments  have  found  that  this  approximation  provides  good  agreement  with  observed 
velocity  profiles.  Above  the  constant-stress  region,  the  velocity  profile  is  relatively 
slowly-varying  and  the  magnitude  is  insensitive  to  the  form  of  the  eddy  viscosity. 
Ellison  (1956)  and  Businger  and  Arya  (1974)  suggest  a  form  for  the  eddy  viscosity 
which  is  linear  close  to  the  bottom,  but  which  is  modulated  by  an  exponential 
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Ut  =  Ku,ze  s 


(3.47) 


The  eddy  viscosity,  then,  reaches  a  maximum  a.t  z  =  h,  and  decays  as  z  >  oo; 
that  is,  approaching  the  top  of  the  boundary  layer.  Long  (1981)  solved  the  neutral 
momentum  equations  (Equations  3.12  and  3.13)  using  the  exponential  form  for  the 
eddy  viscosity  and  found 

,  lu, 

6  / 

where  y  represents  his  boundary  layer  scale  height.  Above  the  constant-stress 
layer,  Ekman  turning  becomes  significant,  so  that  the  direction  of  flow  is  not  the 
same  as  the  direction  of  the  bottom  stress. 

Numerical  solutions  for  the  velocity  and  concentration  profiles  in  the  outer 
boundary  layer  have  been  approached  by  Long  (1981)  and  Glenn  (1983).  In  the 
present  work,  an  estimate  of  outer  Ekman  layer  load  and  transport  which  can  be 
solved  analytically  is  used  instead.  Following  Smith  and  Long  (1976),  the  assump¬ 
tion  is  made  that  the  eddy  viscosity  above  the  constant  stress  layer  can  be  set 
to  a  constant  value.  This  leaves  open  the  choice  of  what  height  to  choose  as  the 
beginning  of  the  outer  Ekman  layer  and  what  eddy  viscosity  to  apply  there. 

Under  the  assumption  that  transport  of  sediment  is  controlled  by  the  most 
energetic  eddies,  a  logical  choice  for  the  outer  Ekmzm  layer  eddy  viscosity  is  the 
maximum  value  determined  using  the  exponential  form  (Figure  3.1).  This  value 


occurs  at 


z  —  h  —  ^  ^ 


At  this  height,  the  value  of  the  eddy  viscosity,  as  defined  by  the  exponential  model. 
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The  height  where  the  linear  model  produces  this  same  value  will  be  defined  as  the 
top  of  the  near-bottom  region.  Using  the  linear  eddy  viscosity  model,  this  value  of 
the  eddy  viscosity  is  reached  at  a  height 


Figure  3.1  shows  the  linear  (Equation  3.16),  the  modulated  (3.47,  and  the  combined 
eddy  viscosity  profiles.  Using  the  combined  (linear  near-bottom  and  constant  outer 
Ekman  layer)  eddy  viscosity  model,  the  eddy  viscosities  in  the  outer  Ekman  layer 


are  defined.  For  the  neutral  case 


I't.Ek  =KU,—  ~ 
o 


(<cu«)^ 

6/ 


and  for  the  stratified  case 


t^tm.Ek  —  ^ 

¥ 

t^t,Ek 

l^U.Ek  —  ^  • 

Tf 


(3.48) 


(3.49) 


(3.50) 


where  the  Ek  subscript  refers  to  the  constant  value  used  in  the  outer  Ekman  layer 
and  refers  to  the  value  of  the  Monin-Obukov  length  at  the  top  of  the  constant 
stress  layer. 

Use  of  the  constant  eddy  viscosity  in  the  outer  Ekman  layer  means  that  for  the 
outer  boundary  layer,  the  sediment  continuity  equation  (Equation  3.31)  takes  the 
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I'tM.Ek-^  +  =  0 


(3.51) 


This  equation,  solved  with  a  boundary  concentration  given  at  z  =  ^  has  the 


solution; 


=  <^(-^)exp[--^z'] 

O  t'tf.Ek 


(3.52) 
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Figure  3.1:  Dimensionless  eddy  viscosity  as  a  function  of  nondimensional  distance  from 
the  bottom.  Solid  curved  line  is  exponentially  decaying  model;  dashed  line  shows  constant 
outer  layer  value;  straight,  solid  line  is  linear  model  (after  Glenn, 1983) 


where  z'  =  z  —  Rather  than  calculate  the  concentration  profile,  the  fact  that 
the  formula  has  an  exact  solution  when  integrated  from  zero  to  infinity  is  used  to 
estimate  the  entire  load  in  the  outer  boundary  layer: 


z-]dz'  =  C(|) 


Wf 


To  estimate  the  volume  of  transport  in  the  outer  layer,  we  multiply  the  total 
load  by  the  velocity  calculated  at  a  =  since  this  height  represents  the  limit 
of  the  constant  stress  layer  where  velocity  varies  rapidly.  The  velocity  profile  is 
approaching  the  geostrophic  value  at  this  point  and  velocity  growth  is  slow,  so 
that  matching  estimated  values  for  the  logarithmic  profile  and  the  velocity  deficit 
profile  calculated  from  the  geostrophic  boundary  layer  equations  gives  an  estimate 
of  error  of  no  more  than  25%  and  a  20-degree  maximum  turning  angle,  which  is 
acceptable  given  the  fact  that,  under  most  circumstances,  transport  is  expected  to 
be  concentrated  in  the  near  bottom  layer. 

This  estimate  has  two  significant  applications.  First,  the  relative  significance  of 
the  sediment  load  and  transport  in  the  outer  Ekman  layer  can  be  assessed.  While 
in  most  cases  sediment  is  expected  to  remain  concentrated  in  the  lower  part  of  the 
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boundary  layer,  fine  sediment  in  a  strong  current  could  diffuse  upward  in  significant 
concentrations  much  higher  in  the  boundary  layer.  In  these  cases,  it  is  necessary 
to  include  Ekman  turning  to  represent  accurately  the  direction  of  transport,  and 
results  of  the  model  should  be  interpreted  cautiously.  This  estimate  shows  when 
Ekman  turning  can  be  neglected  for  calculation  of  transport;  that  is,  when  the 
constant  stress  layer  model  described  and  used  in  this  thesis  is  sufficient. 

On  the  other  hand,  the  load  and  transport  estimates  are  accurate  enough  for 
estimating  net  transport  as  the  model  is  used  here.  The  depth  of  reworking  pre¬ 
dicted  under  given  wave-current  conditions  may  be  increased  by  the  outer  Ekman 
layer  load;  this  technique  provides  an  estimate  of  that  depth.  When  calculating  net 
transport  (erosion  or  deposition)  for  an  area,  the  conditions  will  be  slowly-varying 
enough  from  one  grid  point  to  the  next  that  the  outer  Ekman  layer  transport  di¬ 
rections  and  volumes  should  be  similar.  Most  of  the  error  in  direction  and  volume 
int  roduced  by  assuming  that  transport  is  in  the  direction  of  and  controlled  by  the 
near- bottom  steady  current  is  therefore  cancelled  in  calculating  the  net  transport 
for  a  grid  square. 

3.3  Bottom  Roughness 

The  calculation  of  the  velocity  profile  inside  the  wave  boundary  layer  (Equa¬ 
tion  3.28)  depends  explicitly  on  the  physical  bottom  roughness  length  zq-  This 
length  is  also  necessary  for  the  calculation  of  the  friction  factor,  feu,  on  which  the 
calculation  of  the  shear  velocities  depends.  The  roughness  changes  in  response  to 
the  flow,  due  to  bedform  development  and  the  formation  of  a  layer  of  sediment  in 
motion. 

The  model  used  in  this  work  for  movable  bed  roughness  under  a  combined  wave 
and  current  flow  was  developed  by  Grant  and  Madsen  (1982);  that  work  will  be 
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described  briefly  here. 

The  physical  roughness  felt  by  the  near-bottom  flow  is  the  sum  of  three  com¬ 
ponents:  1)  the  roughness  due  to  the  individual  grain  diameters  in  the  bed  (skin 
friction);  2)  the  roughness  due  to  ripples  and  mounds  on  the  seafloor  (form  drag); 
and  3)  a  roughness  associated  with  dissipation  due  to  moving  sediment  in  the 
near-bed  layer.  The  effect  of  these  elements  will  be  parameterized  in  terms  of  a 
Nikuradse  equivalent  sand  grain  roughness,  fcj.  As  used  here,  the  roughness  height 


is  expressed  as: 


kb  =  kb^n  +  kb^B  +  kbx 


(3.53) 


where  the  three  terms  on  the  right  hand  side  represent  the  roughness  due  to  grains, 
ripples,  and  sediment  transport,  respectively. 

The  grain  roughness,  /cj  „,  is  represented  by  the  grain  diameter  d.  For  a  flat 
bed,  the  grain  roughness  is  the  only  roughness  element,  and  the  skin  friction  is 
the  total  roughness.  In  most  continental  shelf  situations,  however,  there  are  either 
hydrodynamically  or  biologically  generated  roughness  elements  at  least  an  order  of 
magnitude  greater  than  the  grain  size,  so  that  grain  roughness  can  be  neglected. 
It  should  be  noted,  however,  that  the  sand  grain  size  is  the  appropriate  roughness 
length  for  the  skin  friction  component  of  the  total  boundary  shear  stress,  on  which 
initiation  of  motion  and  bedload  calculations  depend. 

The  form  drag  component  of  shear  stress  is  generated  by  the  formation  of 
eddies  in  the  wake  of  the  roughness  element  and  the  reattachment  of  the  flow 
between  elements.  The  roughness  is  dependent  on  the  shape  and  distribution  of 
the  elements.  Using  the  analysis  of  roughness  elements  of  Wooding  et  al.  (1973), 
Grant  and  Madsen  (1982)  derive  an  expression  for  roughness  associated  with  a 
two-dimensional  wave-generated  ripple: 


kb,B  =  27.7r?(^) 


(3.54) 
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where  t]  and  A  are  ripple  height  and  length.  Grant  and  Glenn  (1983)  include  an 
expression  for  a  generalized  roughness  where  the  height,  shape  and  concentration 
of  the  elements  axe  known.  It  is  not  used  in  the  model  since  that  information  is 
not  generally  available;  as  information  about  the  seafloor  becomes  more  detailed, 
that  element  of  the  model  can  be  determined  more  precisely. 

The  dimensions  of  the  ripple  are  determined  best  from  direct  observation  of  the 
seafloor.  When  this  is  impossible,  or  when  the  roughness  is  in  transition,  empirical 
bedform  formulas  can  be  used.  The  model  used  in  this  work,  since  the  cases  of 
interest  are  wave-dominated,  is  a  model  of  wave-generated  ripples  discussed  in 
Grant  and  Madsen  (1982).  For  boundary  shear  stress  only  slightly  greater  than 
that  needed  to  initiate  motion.  Grant  and  Madsen  found  that  ripples  change  slowly 
in  what  they  refer  to  as  ‘equilibrium  range’.  At  some  higher  shear  stress,  ripples 
height  decreases  rapidly  as  they  aire  washed  out.  The  shear  stress  at  which  the 
washing  out  process  begins  is  designated  by  a  breakoff  Shields  parameter; 


V'B  =  i.aso-Vc 


(3.55) 


where  •?«  is  a  non-dimensionalized  measure  of  grain  diauneter  (Equation  2.8)  and  xpe 
is  the  critical  Shields  parameter  for  initiation  of  motion.  The  empirical  relationships 
for  ripple  geometry  under  waves  given  by  Grant  and  Madsen  are: 
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These  values  are  used  in  Equation  3.54  to  calculate  ripple  roughness  in  the  model. 
Larger  scale  bed  roughness  which  would  not  be  perceived  by  the  wave,  but  which 
would  affect  the  mean  flow,  is  not  treated  here,  as  it  would  be  more  appropriate  to 
a  general  circulation  model. 

The  roughness  associated  with  sediment  transport  is  based  on  arguments  ad¬ 
vanced  by  Owen  (1964)  that  the  waice  structure  around  sediment  grains  in  the 
near-bed  transport  layer  cause  the  flow  to  feel  a  roughness  proportional  to  the 
thickness  of  the  layer.  This  concept  was  applied  by  Smith  amd  McLean  (1977)  to 
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steady  flow  in  the  Columbia  River,  and  by  Grant  and  Madsen  (1982)  to  oscillatory 


Grant  and  Madsen  derive  an  expression  for  the  layer  thickness  by  balancing  the 
initial  kinetic  energy  of  a  particle  put  into  motion  with  the  potential  energy  at  its 
highest  elevation.  The  roughness  length  they  derive  using  data  from  Carstens  et 
al.  (1969)  is  expressed: 


hj  =  160(s  +  Cm)dtl>c[{~)^  -  0.7]' 


(3.56) 


where  =  0.5  is  the  coefficient  of  added  mass  of  a  sphere. 

The  roughness  length  zq  in  fully  turbulent  flows  as  considered  in  the  boundary 
layer  model  is  equal  to  or  approximately 


ivj)  +  5-3(s  +  C„)dxl}c\{~)  '^  -  0.7]*  +  — 


(3.57) 


The  three  terms  on  the  right  hand  side  represent,  respectively,  ripple  roughness 
(^,rip)t  sediment  transport  roughness  {z<j,$.t.),  and  grain  size  roughness. 


3.4  Sediment  Transport 

The  Grant-Madsen-Glenn  model  is  concerned  primarily  with  prediction  of  velocity 
and  concentration  profiles.  The  method  used  to  apply  these  to  the  predictions  of 
suspended  sediment  transport  and  a  method  for  estimating  bedload  transport  are 
described  here;  these  applications  are  not  explicitly  covered  in  GMG. 

Three  quantities  derived  from  the  boundary  layer  predictions  are  examined  in 
the  sensitivity  analysis  and  applications:  1)  the  net  transport  for  a  given  location, 
which  can  be  used  to  determine  whether  a  given  area  is  undergoing  erosion  or 
deposition;  2)  the  total  transport  over  time,  which  can  be  used  to  estimate  the 
residence  time  of  sediment  on  the  shelf;  and  3)  the  reworking  depth,  which  demon¬ 
strates  the  degree  to  which  sediment  is  mobilized  during  a  transport  event,  and 
which  may  prove  useful  in  interpretation  of  the  conditions  necessary  to  generate 
observed  geological  strata.  The  first  can  be  calculated  by  integrating  the  transport 
calculations  over  time: 

q,,tot  =  j Cudz)  +  qi4d\dt  (3-58) 

where  the  first  term  on  the  right  hand  side  is  suspended  transport,  calculated  us- 
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ing  velocity  and  concentration  profiles  as  discussed  above,  and  the  second  term  is 
bedload  transport.  Because  this  quantity  is  integrated  over  time,  any  errors  in  the 
transport  calculations  are  multiplied  as  the  time  interval  is  increased.  Integration 
over  time  also  introduces  an  assumption  of  spatial  homogeneity  in  sediment  concen¬ 
tration  profiles  over  the  area  from  which  sediment  is  being  advected  into  the  region 
being  modelled.  For  these  reasons,  estimates  of  total  transport  are  generated  here 
only  for  fairly  short  time  periods,  such  as  the  duration  of  a  single  storm. 

The  net  transport  is  the  divergence  of  the  transport  vector,  defined  in  Equa¬ 
tion  3.4,  which  is  calculated  numerically  for  a  point  inside  the  grid  square  illustrated 
in  Figure  2.6,  using  transport  values  calculated  at  the  corners: 

^  ^  d  d 

^  l/g»,x(a:  + Ax,y)  +  +  Ax,y  +  Ay)  -  q,,x{x,y)  -  y -h  Ay) 

2^  Ax 

g«,v(^.v)  +g«.v(a^  + Ax,y)  -  g«.v(x,y  +  Ay)  -  q,,y(x -h  Ax,y  +  Ay) 

Ay 

where  the  subscripts  denote  the  direction  of  transport,  the  values  in  parentheses 
designate  the  grid  point,  and  Ax  and  Ay  represent  the  grid  spacing  between  points 
in  the  x—  and  y—  directions,  respectively  (Figure  2.6).  The  resultant  transport 
rate  is  a  small  value  reflecting  the  small  difference  between  relatively  large  numbers. 
Considering  the  uncertainties  in  the  component  models  and  data  inputs,  the  error 
would  have  to  be  expected  to  be  large.  The  flux  should  be  integrated  over  only 
relatively  small  time  periods-a  shelf  storm,  for  example-  and  the  prediction  should 
be  interpreted  as  classifying  the  region  as  erosive,  depositional  or  steady-state  and 
suggesting  order  of  magnitude  of  transport  rate.  A  more  precise  interpretation  of 
the  net  transport  predictions  would  be  spurious. 


3.4.1  Bedload  Calculation 


Bedload  is  not  expected  to  be  a  significant  portion  of  the  total  load  in  the  wave- 
dominated  shelf  conditions  on  which  this  study  focuses  except  in  cases  where  trans¬ 
port  is  relatively  small,  and  GMG  does  not  consider  this  transport  element.  Since 
many  of  the  world’s  continental  shelves  are  sand-covered,  however  (Emery,  1968), 
there  is  a  significant  number  of  shelves  where  bedload  makes  up  a  large  percentage 
of  the  transport  under  most  conditions. 

For  this  reason,  a  method  of  estimating  the  relative  importance  of  bedload  to 
the  quantity  and  direction  of  transport  is  developed  here,  using  a  semi-empirical 
bedload  formula.  While  the  order  of  magnitude  and  direction  of  transport  cal¬ 
culated  using  this  method  are  reasonable,  they  should  be  considered  only  rough 
estimates  for  purposes  of  characterization  of  transport  regimes.  The  focus  of  this 
study  is  near-bottom  suspended  transport. 

Empirical  bedload  formulations  often  base  their  predictions  on  estimates  of  the 
skin  friction  component  of  boundary  shear  stress,  raised  to  some  power.  The  Meyer- 
Peter  and  Muller  (1948)  formulation,  an  empirical  formula  based  on  an  extensive 
set  of  laboratory  experiments  in  steady  flow,  raises  the  excess  shear  stress  to  the  | 


power: 


dJ(s  -  l)gd 


(3.59) 


where  q,,bed  is  bedload  transport  in  cm^ fcmfsec  8md  the  prime  indicates  that  the 
Shields  parameter  is  calculated  using  the  skin  friction  component  of  the  boundary 
shear  stress.  To  apply  this  or  any  other  bedload  formula,  which  is  formulated  for 
steady,  unidirectional  flow,  to  the  combined  wave-current  flow,  we  assume  that 
the  response  time  for  the  sediment  is  small  relative  to  the  unsteady  time  scale  (as 
demonstrated  in  Madsen  and  Grant,  1976)  . 

For  this  study,  a  modified  version  of  the  Meyer-Peter  and  Muller  bedload  for- 
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mula  (Equation  3.59)  is  used.  For  ease  in  calculation,  it  is  assumed  that  if  the 
maximum  skin  friction-based  Shields  Parameter  is  less  than  the  critical  value 
for  initiation  of  motion,  then  there  is  no  bedload  transport.  However,  if  the  maxi¬ 
mum  is  greater  than  the  critical  we  make  the  conservative  assumption  that  trans¬ 
port  occurs  throughout  the  wave  cycle,  in  proportion  to  the  instantaneous  bottom 
boundary  shear  stress.  This  method  is  conservative  in  that  it  will  tend  to  overpre¬ 
dict  the  bedload  transport. 

The  instantaneous  bedload  transport  is  thus  calculated  from  : 

— (3.60) 
d^{s  -  l)gd 

where  tp'  (Equation  2.7)is  proportional  to  the  skin  friction  component  of  the  in¬ 
stantaneous  shear  stress  Tq,  which  is  calculated  eis  in  Equation  3.15: 


where  is  the  skin  friction  component  of  the  wave-current  friction  factor  and  u 
is  the  instantaneous  wave  plus  current  velocity  in  the  wave  boundary  layer,  which 
varies  over  the  wave  period.  To  estimate  the  time-averaged  bedload  transport,  the 
velocity  terms  are  raised  to  the  |  power  and  time-averaged  over  a  wave  period. 
The  X-  and  y-components  of  the  instantaneous  velocity  are  defined: 

u  =  Uj  cos  art  +  Ua  cos  <j>c 
V  =  Ua  sin  <f>e 

where  <i>e  is  the  angle  between  the  wave  (which  is  defined  as  the  x-direction)  and 
the  current,  and  cosu;t  defines  the  phase  of  the  wave. 


Inserting  these  definitions  into  the  shear  stress  equation  raised  to  the  |  power 
gives: 
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Chapter  4 


Model  Background:  Calculation  Method 


The  interconnections  of  these  disparate  elements,  and  the  generation  of  the  results, 
can  be  illustrated  best  by  using  a  step-by-step  examination  of  how  the  boundary 
layer  model  works.  The  programs  used  in  the  application  of  the  model  were  de¬ 
veloped  by  Scott  Glenn  (1983),  and  modified  slightly  for  the  present  work.  The 
near-bottom  model  Fortran  code  can  be  found  in  Gr2mt  and  Glenn  (1983c).  The 
computational  procedure,  as  discussed  in  the  preceding  sections  and  applied  here, 
is  traced  in  Figure  4.1.  Each  line  in  the  fiow  chart  is  labelled,  and  those  labels  are 
referred  to  in  this  discussion. 

There  are  three  inputs  (Line  1)  to  the  model  at  each  point:  (l)  current  ve¬ 
locity  (ur)  at  some  height  within  the  current  boundary  layer  zuid  above  the  wave 
boundary  layer,  (2)  wave  conditions,  consisting  of  maximum  wave  bottom  velocity 
(ui)  and  wave  excursion  amplitude  (At)  (or,  equivalently,  wave  height  [H)  and 
period  (T)  and  water  depth  (/i)),  and  (3)  sediment  size  (d),  density  (p,)  and  tex¬ 
ture.  For  simplicity,  co-directional  wave  and  current  and  a  single  grain  size  bed 
are  assumed  in  this  discussion.  The  example  presented  is  a  moderate  storm  wave, 
with  a  26  cm/sec  current  measured  one  meter  above  a  bottom  composed  of  coarse 
silt  (Table  4.1). 

The  first  step  in  the  model  is  to  estimate  the  contribution  of  the  current  to 
boundary  shear  stress  on  the  bottom  (Line  2a).  This  is  represented  zis  where 


Input  parameters  for  sample  model  run 
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Ah 

96  cm. 

H 

2.6  meter 

T 

15  seconds 

h 

50  meters 

la 

0.002 

Ur 

26^ 

Zr 

1.0  meter 

d 

0.006  cm. 

Table  4.1:  Input  parameters  for  sample  model  run  discussed  in  text 

Ua  represents  the  mean  velocity  at  some  unspecified  height  within  the  wave  bound¬ 
ary  layer.  The  parameters  a  and  V2,  used  to  calculate  the  shear  velocities  in 
Equations  3.22  and  3.23,  are  direct  functions  of  The  model’s  initial  estimate  is 

Ua  _  Uj^ 

Ufc  u&  * 

This  shear  stress  estimate  is  used  with  the  grain  roughness  (d)  to  calculate 
the  skin  friction  component  of  the  friction  factor  Line  2b;  Equation  B.l). 
This  friction  factor  is  necessary  to  test  for  initiation  of  sediment  motion.  It  is 
used  with  Equation  3.22  to  calculate  the  skin  friction  component  of  the  maximum 
bottom  shear  stress  Tq  and,  from  that,  the  maximum  Shields  parameter  for  the 
flow  (Equation  2.7),  If  the  Shields  parameter  is  less  than  the  critical  value  for  the 
sediment  on  the  seafloor,  no  sediment  moves  (Lines  2c-  2e).  In  that  case,  the  skin 
friction  shear  stress  is  the  same  as  the  total  shear  stress,  unless  some  roughness 
element  representing  pre-existing  bedforms  has  been  included  as  input.  For  now, 
we  use  only  the  movable  bed  roughness  estimates. 

If  sediment  is  moving,  then,  the  boundary  roughness  due  to  ripples  zmd  sediment 
transport  is  calculated  according  to  Equation  3.57,  and  that  roughness  is  used 
in  Equation  B.l  to  calculate  the  total  friction  factor  fcv  (Line  2f).  The  total 
fetu  is  used  to  calculate  the  mean  and  maximum  shear  velocities  (Equations  3.22 
and  3.23;  Line  2g).  From  these,  the  first  estimate  of  the  predicted  reference  velocity 
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Figure  4.1:  Flow  chart  tracing  computational  procedure  for  boundary  layer  model.  Labels 
are  referred  to  in  text 
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is  calculated  using  Equation  3.29  (Line  2h).  If  the  predicted  velocity  is  not  within 
one  per  cent  of  the  given  reference  velocity  (as  it  certainly  will  not  be  on  the  initial 
try),  the  model  chooses  another  value  of  and  proceeds  again  through  the  steps 
just  described. 

If  the  predicted  value  was  too  high,  the  value  of  ^  is  halved;  if  too  low,  the 
parameter  is  multiplied  by  a  factor  of  2.05  (this  factor  is  chosen  rather  than  2.0 
to  avoid  a  halving-doubling  flip-flop  in  iterations).  Iterations  continue  until  the 
predicted  and  given  currents  match  within  one  per  cent.  At  that  point,the  neutral 
velocity  profile  is  calculated  using  Equation  3.29  (Line  2j). 

If  sediment  was  put  in  motion,  the  sediment  concentration  profile  is  calculated, 
first  without  including  stratification  corrections  to  either  velocity  or  concentration 
profiles.  The  particle  fall  velocity  is  determined  from  Equation  2.9  and  the  sedi¬ 
ment  reference  concentration  is  determined  from  Equation  3.46.  These  are  used  in 
Equations  3.35  and  3.36  to  calculate  the  sediment  concentration  profile  (Line  2n). 

Finally,  the  velocity  and  concentration  profiles  are  integrated  to  determine  the 
neutral  load  and  transport  predictions  (Equations  3.2  and  3.3).  The  estimated 
bedload  is  calculated  using  Equation  3.60.  The  estimated  outer  Ekman  layer  load 
and  transport  are  calculated  from  Equation  3.51. 

The  neutral  results  for  the  wave  case  described  above  are  shown  in  Table  4.2. 
Note  that  the  value  of  Ug  drops  by  a  factor  of  three  from  the  first  estimate.  Most  of 
the  roughness  is  generated  by  ripples  (compare  2o,r«p  vs.  2d,».t.)j  and  the  additional 
roughness  beyond  the  grain  diameter  increases  the  friction  factor  significantly  (com¬ 
pare  vs.  few)-  The  wave-current  shear  velocity  is  more  than  twice  the  current 
shear  velocity.  The  predicted  bedload  transport  is  insignificant  compared  with  the 
suspended  transport.  The  outer  Ekman  layer  load  and  transport,  however,  are 
estimated  at  more  than  three  and  a  half  times  the  near-bottom  load  and  transport 
for  this  coarse  silt,  suggesting  that  in  some  cases,  at  least,  more  attention  must  be 
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Neutral  results  for  sample  model  run 
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Table  4.2:  Some  results  for  neutral,  near-bottom  model  run  for  a  moderate  storm  wave  on 
the  continental  shelf,  as  described  in  text.  N.B.  refers  to  near-bottom  load  and  transport 
{z  <  ^).  O.E.  refers  to  the  outer  Ekman  layer  load  and  transport  (z  >  ^).  Input 
parameters  are  shown  in  Table  4.1. 

paid  to  the  outer  Ekman  layer  transport.  More  attention  is  devoted  to  this  ques¬ 
tion  in  the  sensitivity  analysis.  The  predicted  neutral  velocity  and  concentration 
profiles  are  shown  in  Figure  4.2  (a)  and  (b). 

The  stratified  calculation  begins,  as  does  the  neutral  one,  with  an  estimate 
of  Initially,  the  value  that  produced  the  neutral  case  solution  is  used.  This 
results  in  a  prediction  of  reference  velocity  which  is  too  high,  since  stratification 
increases  velocities  above  the  wave  boundary  layer.  The  same  steps  as  for  the 
neutral  case  are  followed  through  calculation  of  the  shear  velocities  at  Line  3f  in 
the  flow  diagram  (Figure  4.1).  At  this  stage,  the  calculation  of  the  stratification- 
corrected  concentration  profile  begins. 

First,  the  concentration  profile  without  the  stratification  correction  is  calcu¬ 
lated,  as  for  the  neutral  case.  That  profile  is  integrated  and  used  in  Equation  3.39 
to  procure  an  initial  estimate  for  the  integrated  Monin-Obukov  Length  (/  j;dz\ 
Lines  3g  and  3h).  This  value  is  used  to  calculate  a  stratification-corrected  con¬ 
centration  profile  (Equation  3.45),  which  is  then  integrated  to  generate  a  revised 
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estimate  of  the  integrated  Monin-Obukov  Length  (Lines  3i-3j).  If  the  difference 
between  the  old  and  new  integrated  values  of  L  is  greater  than  the  allowable  error, 
the  new  value  is  used  to  calculate  another  revised  concentration  profile.  These 
iterations  (Lines  3i-3k)  continue  until  the  integrated  Monin-Obukov  Length  values 
converge. 

Once  the  concentration  profile  is  determined  for  this  ^  value,  the  new  ref¬ 
erence  velocity  prediction  including  stratification  effects  can  be  calculated  using 
Equation  3.44.  As  in  the  neutral  csise,  if  the  predicted  and  given  values  are  outside 
acceptable  error  limits  (which  we  take  to  be  one  per  cent  since  the  mean  current 
is  not  measurable  more  precisely  than  this),  iterations  begin  again  with  a  new  ^ 
value,  revised  in  the  same  manner  described  above  for  the  neutral  case  (Figure  4.1, 
Lines  3a-3m).  Once  the  velocity  values  converge,  the  final  stratified  velocity  and 
concentration  profiles  and  the  transport  and  load  predictions  are  calculated. 

The  stratified  results  for  the  wave  case  described  above  are  shown  in  Table  4.3 
and  Figure  4.2  (c)  and  (d).  The  results  differ  markedly  from  the  neutral  case  results 
shown  in  Table  4.2  and  Figure  4.2  (a)  and  (b).  The  bottom  line  of  Table  4.3  notes 
that  the  neutral  case  ^  value  gave  a  reference  velocity  prediction  of  39  cm/sec. 
That  is,  stratification  led  to  an  increase  in  predicted  velocity  at  one  meter  above  the 
seafloor  of  about  50  %.  The  sharp  drop  in  concentration  and  increase  in  velocity 
above  the  wave  boundary  layer  can  be  seen  by  comparing  the  neutral  and  stratified 
profiles  in  Figure  4.2. 

The  largest  changes  from  the  neutral  case  result  from  the  reduced  current  shear 
velocity:  u.c  is  approximately  one-half  the  neutral  value.  For  this  reason,  ^  drops 
to  9.4  m  in  the  stratified  case  from  16.4  m,  and  the  near-bottom  suspended  load  and 
transport  estimates  drop  by  an  order  of  magnitude  or  more.  Even  more  notaole, 
the  estimates  of  outer  Ekman  layer  load  and  transport  drop  by  two  or  more  orders 
of  magnitude,  so  that  the  near-bottom  load  and  transport  volumes  are  about  three 


Chapter  5 


Sensitivity  Analysis 


The  sediment  transport  and  load  predictions  made  by  the  boundary  layer  model 
depend  on  a  broad  array  of  parameters  and  assumptions  (Table  5.1).  As  described 
in  the  Introduction  and  Background  Sections,  the  influence  of  each  parameter  on 
sediment  transport  is  related  in  a  complex  way  to  the  others,  so  that  the  effects 
of  changing  a  parameter  are  not  readily  predictable.  For  example,  a  decrease 
in  grain  size  under  constant  wave  and  current  conditions  might  suggest  that  the 
volume  of  sediment  load  would  go  up.  However,  the  grain  size  is  strongly  related 
to  the  stratification  correction  and  to  roughness  (sq),  not  so  much  because  the 
grain  roughness  kf,_n  is  important,  but  because  grain  size  governs  ripple  dimensions 
and  sediment  transport  layer  roughness  and  fck.r).  If  conditions  are  such  that 
roughness  decreases  and  stratification  increases,  the  effect  of  decreasing  grain  size 
could  be  only  a  very  small  increase,  or  even  a  decrease,  in  total  load  and  transport. 

The  purpose  of  this  chapter  is  to  establish  the  relative  significance  of  the  vari  ¬ 
ous  parameters  in  the  model’s  predictions  of  load  and  transport.  This  serves  three 
useful  purposes.  First,  the  examples  provide  physical  intuition  on  the  roles  of  pa¬ 
rameters  and  their  interactions.  Second,  this  analysis  demonstrates  strengths  and 
weaknesses  of  the  r.heoretical  model:  it  identifies  circumstances  where  the  model’s 
results  are  physically  reasonable,  but  also  defines  conditions  where  field  and  lab- 


Model  parameters  describing  physical  environment 
H  wave  height 

T  wave  period 

h  water  depth 

Ur  reference  current  velocity  at  100  cm  off  bottom 
dn  grain  size(s) 

Atj  physical  boundary  roughness  =  «o  x  30 

part  of  boundary  roughness  due  to  ripples 
Tj,  A  ripple  height  and  wavelength 
6^  wave  boundary  layer  thickness 

<^e  angle  between  waves  and  current 

Empirical  parameters 

')o  sediment  reference  concentration  parameter 
/?,  7  stratification  parameters 

Implicit  assumptions  in  theory  and/or  model  formulation 
Effects  of  stratification  on  velocity  and  concentration  profiles 
can  be  modeled  by  revising  the  eddy  viscosity  to 

The  wave  spectrum  can  be  represented  by  a  single  wave  height 
and  frequency,  such  as  the  frequency  of  the  spectral  peak. 

Grain  distribution  can  be  represented  by  1  or  2  modal  grain 
sizes. 

Effects  of  armoring  can  be  modeled  so  that  there  is  equilibrium 
between  calculated  load  and  bottom  concentrations. 

Parameters  vary  slowly  enough  in  time  and  space  so  that  there 
is  steady  state  in  the  vertical  concentration  and  velocity  pro¬ 
files. 

Table  5.1:  Parameters  and  assumptions  in  the  boundary  layer  model 
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f  oratory  experiments  are  needed  to  constrain  the  model’s  predictions.  Third,  it 

i  provides  initial  estimates  of  the  relative  volumes  of  sediment  transport  in  differ- 

I  ent  circumstances,  thus  discriminating  between  the  importance  of  storms  versus 

I  everyday  regional  currents  in  transporting  sediments. 

I 

The  organization  of  this  section  reflects  the  focus  of  this  study  on  the  influence 
of  waves  on  sediment  dynamics.  As  far  as  sediment  response  is  concerned,  the 
general  climate  is  set  by  the  waves,  and  so  the  sediment  transport  results  are 
presented  by  wave  condition.  The  variation  of  transport  and  load  due  to  other 
factors  is  discussed  within  the  context  of  a  particular  set  of  wave  parameters. 

In  the  first  subsection  the  depth  limit  of  wave  influence,  as  imposed  by  wave 
period,  is  discussed.  Each  of  three  subsequent  subsections  treats  a  single  wave 
condition  (represented  by  bottom  velocity  ut  and  excursion  amplitude  ^5)  rep¬ 
resentative  of  a  general  physical  climate.  The  first  of  these  examines  the  effect 
of  moderate  storm  waves  on  transport.  It  is  used  as  a  case  study  of  the  effects 
of  stratification  on  load  and  transport,  and  contains  detailed  explanations  for  the 
observed  changes  due  to  stratification.  Results  for  the  other  five  wave  conditions 
are  discussed  in  the  final  section,  primarily  in  contrast  with  the  three  that  are 
examined  in  detail.  Within  each  subsection,  the  effects  of  variation  of  several  of 
the  environmental  parameters  listed  in  Table  5.1  are  described. 

5.1  Maximum  Depth  of  Wave  Influence 

“Wave  base”,  or  the  maximum  depth  at  which  waves  affect  the  bottom,  is  often 
described  as  one-half  the  deepwater  wavelength  (A).  This  depth  is  derived  from 
the  dispersion  relationship  of  linear  wave  theory,  which  relates  wave  period  (T), 
water  depth  (A),  and  wavelength  (see,  e.g.,  LeBlond  and  Mysak,  1978): 

w*  =  gkia.Tihkh  (5.1) 
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Wave  height  H  (meters)  for  Ub=10  cm/sec 
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Wave  period  T  (seconds) 

Figure  5.1;  The  approximate  maximum  depth  of  influence  of  waves,  calculated  as 
one-half  A.  Top  axis  shows  wave  height  required  to  generate  a  10  cm/sec  maximum  wave 
velocity  at  that  depth 

In  this  equation,  w  =  y  is  wave  frequency,  g  is  gra-vity,  and  fc  =  ^  is  the  wavenum¬ 
ber.  For  short  waves  or  deep  water  {kh  >•  1),  the  hyperbolic  tangent  approaches 
unity.  The  wavelength  in  deep  water  can,  therefore,  be  calculated  directly  from 
the  wave  period,  as  the  dispersion  relationship  reduces  to: 

A  =  ^T‘  (5.2) 

This  formula  is  accurate  within  one  per  cent  for  hjL  >  j;  shallower  than  this, 
water  depth  becomes  significant  in  calculating  A.  This  ‘wave  base’  (Figure  5.1) 
can  be  used  in  conjunction  with  the  formula  for  calculating  wave  bottom  velocities 
(Equation  3.24)  to  establish  a  rough  criterion  for  maximum  depth  of  influence  for 
waves. 

The  significance  of  the  wave  for  sediment  transport  is  dependent  on  the  maxi¬ 
mum  wave  bottom  velocity,  U6.  The  maximum  wave  bottom  velocity  at  the  depth 
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/i  =  j  (Figure  5.1)  is  calculated,  using  Equation  3.24, 

Tsinh^f  “  r 

To  sustain  a  particular  bottom  velocity  at  a  depth  of  |A,  therefore,  wave  height 
must  increzise  linearly  with  period.  The  wave  height  necessary  to  generate  a  velocity 
of  10  cm/sec  is  plotted  on  the  top  axis  of  Figure  5.1.  (A  5  cm/sec  velocity  would 
require  waves  one-half  as  high;  20  cm/sec  would  require  waves  twice  as  high,  etc.) 
Using  this  graph,  one  can  quickly  determine  whether  a  particular  wave  would 
significantly  affect  the  seafloor  at  a  particular  depth.  For  example,  for  a  wave 
to  generate  a  bottom  velocity  of  10  cm/sec  at  a  depth  of  one  hundred  meters, 
it  must  be  at  least  4.2  meters  high  and  have  a  period  of  eleven  and  one  half 
seconds.  Deeper  water  would  require  higher  and  longer  waves.  This  constraint 
is  of  particular  importance  on  the  central  and  outer  continental  shelf,  where  only 
swell  and  very  large  storm  waves  would  affect  the  bottom. 


5.2  Sensitivity  Test  Conditions 

Wave  conditions 

The  boundary  layer  model  was  run  for  8  different  wave  conditions  (Table  5.2)  to 
represent  a  broad  range  of  wave  conditions  that  might  affect  sediment  transport 
on  the  continental  shelf.  The  primary  consideration  was  to  test  a  range  of  wave 
bottom  velocities.  The  lowest  velocity  chosen  was  20—,  because  lower  velocities 
set  little  or  no  sediment  in  motion.  Wave  bottom  velocities  were  increased  by 
20^  steps,  and  at  least  two  wave  periods  were  used  for  each  bottom  velocity  in 
order  to  demonstrate  the  effect  of  period  on  sediment  transport  predictions.  One 
case  is  presented  for  =  100^,  to  demonstrate  the  results  under  extreme  wave 
conditions.  Wave  periods  from  10  to  20  seconds  were  used  in  order  to  encompass 


fVave  Chmates  Tested 


Waves 


Moderate  windsea 
Low  swell 
Large  storm,  early 
Moderate  storm,  late 
Moderate  swell 
Large  storm,  late 
Large  swell 
Extreme  swell 


•t  cm 

Ai,(cm) 

T{sec) 

H{\n  m,if  h 

20 

32 

10 

2.5 

20 

64 

20 

1.1 

40 

64 

10 

5.0 

40 

96 

15 

2.6 

40 

128 

20 

2.2 

60 

144 

15 

3.9 

60 

191 

20 

3.3 

100 

319 

20 

5.4 

Fine  sand 
Very  fine  sand 
Coarse  silt 
Medium  silt 


0.0005 


0.001 


Grain  sizes  tested 
d{cm) 


iwiin 


<70  values  tested 


0.002  0.003  0.004 


.965 
.950 
0.250 
0.043 


0.005 


Table  5.2:  Variation  of  parameters  in  sensitivity  analysis 


the  range  of  waves  affecting  most  shelf  depths. 

Each  wave  bottom  velocity  and  period  condition  is  given  a  qualitative  label,  to 
reflect  an  example  of  the  surface  wave  conditions  which  could  impose  the  bottom 
parameters  (Table  5.2).  These  labels,  from  moderate  windsea  zmd  low  swell  to 
extreme  swell  conditions,  do  not  signify  the  only  surface  conditions  applicable  to 
the  bottom  wave  conditions  they  describe.  Both  stratified  and  neutral  predictions 
are  made  for  every  input  condition. 


[%T*: 


Sediment  grain  sizes 


For  each  wave  condition,  the  model  was  run  using  4  sediment  grain  size  classes 
ranging  from  fine  sand  to  medium  silt  (Table  5.2).  These  represent  grain  sizes 
typically  found  on  continental  shelves.  They  also  cover  a  size  range  that  is  sus¬ 
ceptible  to  suspended  transport,  but  not  subject  to  cohesive  forces,  where  cohesion 
is  defined  as  the  electrochemical  binding  forces  by  grains  less  than  10~®cm  (Mc- 
Cave,  1984)  which  can  bind  sediments  into  a  plastic  mass.  However,  biologically- 
produced  mucous  can  cause  sediments  to  stick  together,  increasing  the  initiation 
of  motion  criterion  and  fall  velocity.  Both  biological  and  electrochemical  adhesion 
become  increasingly  effective  and  increasingly  common  as  the  mean  grain  size  of 
the  seabed  decreases.  For  these  reasons,  no  sediments  smaller  than  medium  silt 
(d  =  2.5  X  10“®cm)  are  included  in  the  sensitivity  analysis,  since  fEw:tors  not  con¬ 
sidered  in  this  model  are  so  important  in  determining  such  fundamental  properties 
as  fall  velocity  zmd  initiation  of  motion  for  smaller  grains  (McCave,  1984). 

Although  sands  and  medium  silts,  especially  the  latter,  are  not  completely 
immune  to  adhesive  forces,  it  is  assumed  in  this  study  that  binding  effects  are 
limited  to  an  increase  in  the  critical  Shields  parameter  for  the  initiation  of  motion, 
which  directly  affects  only  the  sediment  reference  concentration  (Equation  3.46). 
The  same  effect  on  load  and  trzmsport  predictions  is  thus  achieved  by  varying  the 
sediment  reference  concentration  parameter  '/o- 

Reference  concentration  parameter 

Sediment  load  and  transport  depend  not  only  on  grain  size  but  also  on  how  much 
sediment  is  mobilized  by  the  bottom  shear  stress,  expressed  in  this  work  as  a 
reference  concentration,  C{zo)  (Equation  3.35).  The  reference  concentration  is 
directly  dependent  on  the  empirical  parameter  7o(Equation  3.46).  Estimated  values 
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for  70  range  from  5  x  10~'*to  4  x  10"*  (Glenn,  1983;  Wiberg  and  Smith,  1983).  To 
gauge  the  effect  of  varying  the  reference  concentration  boundary  condition,  the 
model  was  run  using  six  values  of  70 1  representing  an  order  of  magnitude  variation 
for  each  grain  size  (Table  5.2).  This  variation  in  70  also  reflects  the  uncertainty 
in  all  parameters  involved  in  calculating  the  concentration  at  the  top  of  the  wave 
boundary  layer,  including  initiation  of  motion  criteria  (V’c)>  exact  values  of  the  skin 
friction  component  of  boundary  shear  stress  {tq  amd,  therefore,  xjj'),  and  the  height 
of  the  wave  boundary  layer. 

Current  velocities 

Finally,  for  each  wave  condition,  for  each  grain  size,  and  for  each  reference  concen¬ 
tration,  sediment  transport  and  load  were  calculated  for  a  range  of  current  shear 
stress  values.  This  variation  is  expressed  as  a  reference  velocity  calculated  at  one 
meter  above  the  seafloor,  and  ranges  from  a  few  centimeters  per  second  to  30- 
100  cm/sec  for  each  grain  size  and  wave  condition.  The  reason  for  the  variation 
in  reference  velocity  range  from  case  to  case  is  that  the  velocity  is  an  intermediate 
result  of  the  model;  the  results  reflect  variation  in  bottom  shear  stress,  which  is 
translated  into  a  velocity  prediction,  as  discussed  below. 

Presentation  of  results:  format 

In  the  interest  of  brevity,  full  discussions  of  the  results  are  presented  for  only  three 
of  the  wave  conditions  investigated;  the  other  C2ise8  are  discussed  more  generally  in 
relation  to  those  three.  The  three  wave  conditions  detailed  are  listed  in  Table  5.2 
as  low  swell,  moderate  storm,  and  large  swell.  Fot  each  of  those  three  cases,  a  set 
of  16  plots  is  presented  that  demonstrates  the  effects  of  stratification  on  predicted 
load  and  transport,  and  of  the  variation  of  grain  sizes,  reference  concentration,  aud 
current  velocity  on  both  neutral  and  stratified  load  and  transport  . 


The  contour  plots  that  serve  as  the  fundamental  illustrations  for  these  sections 
embrace  a  great  deal  of  information;  unfortunately,  they  are  not  readily  interpreted. 
A  detailed  explanation  of  the  format  is  presented  here,  referring  to  Figures  5.2-5. 5 
as  examples.  Each  figure  shows  results  of  model  runs  for  moderate  storm  wave 
conditions  (uj  =  40—,At,  =  96  cm),  where  the  seabed  is  covered  with  sediment 
of  a  single  grain  size.  Each  figure  consists  of  four  contour  plots  showing  variation 
in  predicted  near-bottom  load  or  transport  for  that  grain  size  by  that  wave,  in¬ 
cluding  or  neglecting  the  effects  of  sediment  stratification  as  labeled.  The  load  and 
transport  volumes  are  calculated  by  integrating  the  predicted  concentration  and 
velocity  profiles  (Equations  3.2  and  3.3)  to  a  height  6c/6  equal  to  one-sixth  the 
Ekman  depth  (Section  3.2.4).  Bedload  transport  and  outer  Ekman  layer  load  and 
transport  are  not  included  in  these  graphs.  Note  that  the  use  of  ‘near-bottom’  here 
conforms  v/ith  its  use  in  GMG  as  the  lowest  few  meters  of  the  water  column. 

Each  ‘x’  represents  a  load  or  transport  value  calculated  by  the  model.  These  are 
the  data  points  on  which  the  contours  are  based.  The  locations  of  the  points  with 
respect  to  the  x-axis  show  the  value  of  the  concentration  at  the  top  of  the  wave 
boundary  layer  (C(5,u))  for  the  given  set  of  wave,  current,  and  sediment  conditions. 

The  variation  in  load  and  transport  predictions  within  each  plot  results  from 
changing  two  of  the  model’s  parameters:  the  sediment  concentration  parameter 
(70)  and  the  current  contribution  to  the  boundary  shear  stress  within  the  wave 
boundary  layer  (ua).  Rather  than  use  these  parameters  as  the  axes  of  the  plots, 
however,  I  have  chosen  to  use  physically  significant  intermediate  results  of  the 
model’s  calculations  that,  under  most  circumstances,  are  uniquely  determined  by 
the  parameters  being  varied  (that  is,  70  and  Uo).  The  x-axis  (C(fi,„))  represents  a 
potentially  measurable  quantity.  The  y-axis  is  the  current  velocity  at  one  meter 
above  the  bottom,  u^,  non-dimensionalized  by  the  maximum  wave  bottom  velocity 
uj.  The  reference  current  is  a  physical  parameter  used  to  specify  the  problem,  and 


is  used  aa  the  matching  condition  for  the  model  iterations.  That  is,  the  velocity 
calculated  by  the  model,  based  on  its  guess  of  the  bottom  shear  stress  and  strat¬ 
ification  effects,  is  matched  against  the  specified  velocity  to  determine  the  end  of 
calculations  (see  Chapter  7 ).  Its  variation  with  respect  to  input  parameters  is 
therefore  important,  and  circumstances  where  it  is  not  uniquely  determined  by  the 
model  require  an  additional  boundary  condition,  as  discussed  below. 

C{6„)  represents  the  sediment  available  at  the  top  of  the  wave  boundary  layer 
for  transport  by  the  current  (Equation  3.35).  It  is  chosen  rather  than  70  as  the 
x-variable  for  a  number  of  reasons.  First,  it  is  linearly  related  to  7o>  so  that  the 
order-of  magnitude  change  in  70  results  in  a  similar  variation  in  C{6tu).  Unlike  70, 
however,  C(5«,)  is  a  measurable  physical  quantity,  so  that  it  provides  insight  into 
the  predicted  conditions  near  the  boundary,  and  can  be  judged  as  realistic  or 
unrealistic  on  physical  grounds.  Second,  C(S,u)  is  also  related  directly  to  the 
seabed  concentration  of  the  sediment,  and  the  excess  bottom  shear  stress 
5  (Equations  3.35  and  3.46).  Although  the  figures  in  this  section  were  produced 
by  varying  70,  they  could  just  as  well  result  from  varying  the  critical  boundary 
shear  stress,  the  skin  friction  component  of  shear  stress,  or  the  bed  concentration. 
Using  C(5*)  as  the  x-axis  variable,  therefore,  shows  the  sensitivity  of  load  and 
transport  predictions  to  four  variables  that  are  subject  to  uncertainty,  rather  than 
one.  Third,  the  magnitude  of  C(Sm)  indicates  the  degree  of  sediment  response 
to  shear  stress  in  the  wave  boundary  layer.  Compare  the  ranges  of  the  x-axes  in 
Figures  5.2  and  5.4,  for  example.  For  the  sarnie  wave  and  current  conditions,  the 
concentration  of  the  finer  sediment  (Figure  5.4)  at  the  top  of  the  wave  boundary 
layer  is  as  much  ats  twice  that  of  the  couser  sediment  (Figure  5.2). 

Because  the  velocities  represented  by  the  y-axis  are  calculated  by  the  model 
(Equation  3.29  or  3.44),  vertical  ramge  of  the  data  points  varies  within  a  wave 
condition  between  neutral  to  stratified  cases  and  between  grain  sizes.  To  make 
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comparison  easier,  the  y-axis  maximum  is  the  same  for  all  plots  for  a  single  wave 
condition,  although  the  contours  only  extend  as  far  zis  the  limits  of  the  ‘data’.  This 
is  the  reason  for  the  frequent  blank  spots  in  the  figures. 

The  predicted  load  values  are  easily  translated  to  estimates  of  erosion  depth,  or 
reworking  depth  once  the  sediment  settles  back  down  after  the  transporting  event, 
by  dividing  by  the  bed  concentration  (0.6  in  this  work).  In  Figure  5.2,  for  example, 
the  load  of  0.16  ^  at  the  location  circled  on  the  plot  would  signify  a  reworking 
depth  of  0.27  cm. 

5.3  Results:  Moderate  Storm  Waves 

On  some  continental  shelves,  moderately  large  storms  generate  waves  up  to  several 
meters  in  height,  with  periods  of  15  seconds  or  more.  The  wave  discussed  in  this 
section,  with  a  maximum  bottom  velocity  of  40  cm/sec  and  an  excursion  amplitude 
of  96  cm  (denoting  a  15  second  period),  would  could  occur  in  such  a  storm.  In  50 
meters  water  depth,  a  2.6  meter  wave  height  would  be  required  to  generate  that 
velocity. 

The  four  figures  accompanying  this  section  (Figures  5.2-5.5)  are  contour  plots  of 
sediment  load  and  transport,  with  and  without  the  effects  of  stratification  by  sed¬ 
iments  in  suspension.  The  neutral  plots  are  primarily  for  reference:  although  one 
might  be  interested  in  velocity  profiles  for  the  neutral  case,  ignoring  stratification 
by  sediments  is  inconsistent  in  the  calculation  of  a  sediment  profile.  In  some  cases, 
however,  sediment  stratification  effects  are  small,  and  the  iteration  needed  to  solve 
for  the  stratification  effects  produce  immeasurably  small  changes  in  the  concen¬ 
tration  profile.  In  those  cases  the  concentration  profile  above  the  wave  boundary 
layer  can  be  predicted  without  the  stratification  correction  (Equation  3.36).  The 
contrausts  between  the  stratified  and  unstratified  load  and  transport  predictions 


suggest  the  amount  of  error  introduced  by  using  the  simpler  procedure  and  can 
provide  guidelines  for  conditions  under  which  the  neutral  model  is  sufficient. 

5.3.1  Neutral  results 

The  results  for  the  unstratified  cases  (Figures  5.2-5. 5)  are  easily  predictable  from 
the  neutral  equations  for  the  velocity  profiles  (Equations  3.28  and  3.29)  and  con¬ 
centration  profiles  (Equations  3.35  and  3.36).  These  equations  and  the  two  for 
calculating  load  and  transport  (Equations  3.2  and  3.3)  should  be  referred  to  in  the 
discussion  that  follows. 

In  parts  (a)  and  (c)  of  each  figure,  both  the  neutral  sediment  load  and  the  neu¬ 
tral  transport  are  seen  to  depend  directly  on  the  sediment  reference  concentration 
(the  x-axis):  doubling  C(S^)  leads  to  doubling  of  the  load  or  transport  for  any 
velocity.  This  follows  from  the  proportionality  of  the  concentration  profile  to  the 
concentration  at  the  top  of  the  wave  boundary  layer. 

In  the  neutral  case,  a  higher  reference  velocity  always  implies  a  higher  current 
shear  velocity  u.^.  The  shear  velocity  affects  the  sediment  concentration  profile 
in  two  ways:  First,  the  concentration  decays  exponentially  with  —  Higher  u.^ 
values  generate  smaller  concentration  gradients.  Second,  higher  u,«  increases  the 
height  of  the  current  boundary  layer  (4;  Equation  2.6).  This  increases  the  vertical 
distance  over  which  load  and  transport  are  integrated  (5c/6)-  The  slower  decay  of 
concentration  with  height  and  longer  integration  distance  combine  to  increase  tlu- 
predicted  values  of  load  and  transport.  Transport  predictions  depend  directly  c)n 
the  velocity  profile  also  (Equation  3.3),  so  transport  values  rise  more  -^apidly  wnr 
increaising  current  velocity  than  load  values  do. 

The  similar  concave  cj.'ves  of  all  of  the  neutral  plots  (Figures  5.2-5  '>  (ai  and 
illustrate  the  predictability  of  the  response  of  neutral  predictions  to  vana’  •  - 
current  and  reference  concentration.  In  fact,  the  neutral  load  and  tr.oiK;. 
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Figure  5.2.  Predicted  near-bottom  load  (^^y)  and  transport  (cm*/cm/sec)  values  for  fine 
sand  in  moderate  storm  waves.  T=15  sec,  and  the  wave  height  in  50  m  of  water  would  be 
2.6m.  Ur  is  reference  current  velocity.  ranges  from  21-25  cm;  ranges  from  12  to  38 
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Figure  5.5:  Predicted  near-bottom  load  (^)  and  transport  (cm’ /cm/sec)  values  for 
medium  silt  in  moderate  storm  waves.  The  irregular  contour  patterns  for  the  stratified 
plots  occur  in  regions  where  the  model  gives  non-unique  solutions,  as  discussed  in  text. 
ranges  from  8-10  cm;  ^  ranges  from  4  to  18  m. 
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for  all  eight  wave  conditions  tested  share  this  pattern  of  response  to  increasing 
velocity  and  C{Su,)-  The  scale  of  the  response  (i.e.  the  values  of  the  contours)  is 
set  by  wave  conditions  and  grain  size.  Because  of  this  predictability,  the  neutral 
results  are  not  examined  in  detail  in  the  other  discussions  of  results. 

Bedload  and  outer  Ekman  Layer  transport 

In  the  neutral  case  for  the  moderate  storm  wave,  bedload  is  generally  insignificant, 
but  transport  in  the  outer  Ekman  layer  is  often  very  important.  The  contours  in 
the  plots  of  Figures  5. 6-5. 9  show  the  relative  magnitudes  of  the  bedload  transport 
and  outer  Ekman  layer  transport,  compared  with  the  near-bottom  transport.  The 
format  is  the  same  as  that  of  the  load  and  transport  contour  plots,  except  that 
the  velocities  on  the  y-axes  are  dimensional.  These  plots  are  included  only  for  the 
moderate  storm  wave  case,  in  order  to  demonstrate  typical  patterns  of  variation. 
For  other  wave  conditions,  the  results  are  discussed  in  the  text,  but  plots  are  not 
included. 

In  the  neutral  case,  bedload  under  this  strong  wave  is  uniformly  insignificant 
compared  with  the  near-bottom  transport  (plot  (a)  in  Figures  5.6-5.9).  The  largest 
bedload  contribution  is  one  per  cent  of  the  near-bottom  transport;  this  occurs  only 
for  the  coarsest  grain  size  at  the  lowest  current  velocity  and  reference  concentration 
(Figure  5.6(a),  lower  left-hand  corner).  The  relative  value  of  bedload  decreases  with 
decreasing  grain  size  and  with  increeising  current  because,  although  the  predicted 
bedload  volume  rises  with  each  of  those  variations,  the  suspended  transport  rises 
much  faster.  The  ratio  of  bedload  transport  to  near-bottom  transport  in  the  neu¬ 
tral  case  also  decreases  with  increasing  reference  concentration  since  near-bottom 
transport  rises  with  C(6*,),  zis  discussed  above.  Bedload  is  uninfluenced  by  that 
parameter  in  this  formulation  because  it  reflects  an  increase  in  nfoj  of  which  bedload 
is  independent. 
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Figure  5.8:  Relative  values  of  bedload  and  outer  Ekman  layer  transport,  compared  with 
near  bottom  transport,  for  coarse  silt  in  moderate  storm  waves. 
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Figure  5.9:  Relative  values  of  bedload  and  outer  Ekman  layer  transport,  compared  with 
near  bottom  transport,  for  medium  silt  in  moderate  storm  waves.  The  contours  for  the 
stratified  results  demonstrate  the  effects  of  the  nonunique  solutions  to  the  model,  and  are 
discussed  in  detail  in  the  text. 


Outer  Ekman  layer  transport  (Section  3.2.4)  is  much  more  likely  than  bedload 
to  be  a  significant  factor  in  predicting  neutral  transport  volume  (plot  (c)  in  Fig¬ 
ures  5.6-5.9).  Since  both  neutral  near-bottom  emd  outer- Ekman  layer  transports 
are  proportional  to  the  reference  concentration,  there  is  no  variation  with  C(5») 
(the  x-axis)  in  their  relative  values.  This  accounts  for  the  uniformly  flat  contours. 
Since  the  concentration  profile  decays  exponentially  with  —  ^  (Equation  3.36), 
however,  the  relative  value  of  outer  Ekman  layer  transport  rises  with  increasing 
current  and  with  decreasing  grain  size.  An  increase  in  current  from  10  to  40  ^ 
causes  a  three-  to  six-fold  increase  in  the  proportion  of  the  transport  that  occurs 
in  the  outer  Ekman  layer. 

The  relative  rate  of  outer  Ekman  layer  transport  is  strongly  dependent  on  grain 
size,  but  in  a  sense  opposite  to  the  relative  bedload  transport  rate.  The  outer 
Ekman  layer  transport  is  ten  to  sixty  per  cent  of  the  near-bottom  transport  for 
fine  sand  (Figure  5.6(c)),  but  is  eight  to  twenty-eight  times  greater  than  the  near¬ 
bottom  transport  for  medium  silt  (Figure  5.9(c)).  Estimates  of  relative  load  for 
the  outer  Ekman  layer  are  slightly  lower  than  the  relative  transport  rates,  ranging 
from  0.013  to  0.44  for  fine  sand,  0.21  to  1.6  for  very  fine  sand,  0.9  to  4.3  for  coarse 
silt  and  5.6  to  24.6  for  medium  silt. 

These  large  proportions  of  transport  in  the  outer  Ekman  layer  might  be  cause 
for  some  concern,  since  the  present  work  concentrates  on  calculating  the  near¬ 
bottom  sediment  profile.  The  large  proportions  are  tempered  by  four  factors,  how¬ 
ever:  1)  stratification,  as  discussed  in  Section  5.3.2;  2)  sediment  availability,  since 
these  predictions  assume  an  infinite  mixing  depth;  3)  temporal  and  vertical  spatial 
limitations,  since  the  water  depth  may  be  less  than  the  predicted  Ekman  layer 
thickness  or  there  might  not  be  adequate  time  to  establish  equilibrium  sediment 
concentration  conditions  for  the  full  water  colunm  (as  opposed  to  the  near-bottom 
layer)  before  environmental  conditions  change;  and  4)  predicted  suspended  sedi- 


ment  concentration  values,  some  of  which  are  up  to  an  order  of  magnitude  or  more 
greater  than  any  ever  observed  in  continental  shelf  environments. 

The  third  factor  implies  that  the  large  outer  Ekman  layer  load  and  transport 
reflect  a  problem  in  the  method  used  to  calculate  them.  This  is  true  under  some 
circumstances:  the  heights  of  the  ‘near-bottom’  layer  {^)  for  these  cases  are  as 
large  as  38  meters  in  the  stronger  currents  (Figures  5. 2-5. 5).  Obviously,  integrat¬ 
ing  sediment  concentration  over  a  water  column  six  times  that  depth  violates  the 
depth  constraints  imposed  by  the  continental  shelf.  This  problem  can  be  addressed 
simply  by  limiting  integration  to  the  water  depth;  however,  that  introduces  the 
complication  of  merging  of  the  surface  and  bottom  boundary  layers  and  the  need 
to  model  a  fully  mixed  shelf.  That  problem  will  therefore  be  deferred,  particularly 
in  light  of  the  other  considerations. 

The  other  three  reservations  with  respect  to  the  neutral  model  results  reflect 
near-bottom  processes.  Silt  concentrations  of  O(10~*)  are  predicted  in  some  cases 
at  heights  up  to  20  m  or  more  above  the  seafloor;  these  are  one  to  two  orders  of 
magnitude  greater  than  concentrations  observed  during  continental  shelf  storms, 
even  much  closer  to  the  bottom  (e.g.  Butman,  1987a).  These  large  overpredictions 
cannot  be  corrected  simply  by  adjusting  the  reference  concentration  parameter  or 
initiation  of  motion  criteria.  At  least  under  the  conditions  of  this  model  run,  the 
neutral  model  does  not  give  physically  reasonable  results. 

Stratification  offers  one  remedy  for  the  excessive  concentrations,  to  wit,  that  the 
sediment  is  responding  to  a  different  turbulent  structure  than  the  one  the  neutral 
model  produces.  Specifically,  the  stratified  model  works  on  the  assumption  that  the 
sediment  itself  modifies  the  turbulent  structure  and  decreases  the  concentration. 

Bed  armoring  of  a  seafloor  with  mixed  grain  sizes  provides  another  possible 
solution.  If  the  bottom  shear  stresses  and  sediment  responses  occur  zis  modelled 
in  the  neutral  caise,  then  the  flow  would  suspend  the  supply  of  fine  sediments  in 


the  surface  layer,  leaving  a  layer  of  coarser  grains  to  armor  the  surface.  The  large 
differential  in  predicted  reworking  depths  for  different  grain  sizes  under  the 

same  wave-current  conditions  suggests  that  this  could  be  at  least  partly  responsible 
for  the  overprediction  of  sediment  concentrations  in  the  neutral  case. 

5.3.2  Stratification  effects 

In  the  neutral  case,  the  concentration  gradient  decays  with  the  ratio  of  the  fall 
velocity  to  the  shear  velocity.  When  stratification  is  included,  the  concentration 
gradient  is  influenced  not  only  by  —  but  also  by  the  gradient  itself.  In  a  strat- 
ified  flow  with  only  one  length  scale  (a  current  boundary  layer  without  waves,  for 
example),  stratification  is  significant  only  when  the  fall  velocity  is  large  enough  so 
that  the  sediment  is  not  mixed  homogeneously  through  the  entire  constant  stress 
layer,  but  is  not  so  large  that  all  the  sediment  falls  out  very  close  to  the  bottom 
(Glenn,  1983).  In  a  stratified  flow  with  two  length  scales  such  as  those  considered 
here,  the  gradient  in  turbulent  intensity  between  the  wave  boundary  layer  and 
the  current  boundatry  layer  introduces  more  stratification.  The  wave  shear  stress 
mobilizes  sediment  from  the  seafloor  and  transports  it  up  into  the  water  column. 
Over  a  short  distance,  however,  the  wave-gene^'ated  shear  stress  dies  out  and  only 
the  current-generated  turbulence  is  available  to  transport  sediment  higher.  The  a 
sediment  concentration  gradient  depends  on  the  ratio  of  Utcv>  to  u,c,  and  reflected 
in  the  change  in  eddy  viscosity  over  that  short  distance.  The  effect  of  stratification 
on  the  sediment  load  and  transport  predictions  in  this  case  hinges  on  three  factors; 
reference  concentration  C(zo),  the  ratio  of  the  sediment  fall  velocity  to  the  shear 
velocity  (^),  and  the  value  of  the  stability  parameter  (£)  immediately  above  the 
wave  boundary  layer. 

The  effects  of  sediment  stratification  on  sediment  load  and  transport  predictions 
are  shown  in  Figures  5.2-5.5(compare  (a)  and  (c)  with  (b)  and  (d)).  These  effects 


are  less  ejisily  predictable  from  inspection  of  the  governing  equations  than  are 
variations  in  the  neutral  cases.  Initial  inspection  shows  three  primary  changes:  1) 
changes  in  contour  shape,  2)  decrease  in  maocimum  transport  and  load  values,  and 
3)  smaller  blank  areas  at  the  top  of  the  plots.  These  will  be  explained  in  reverse 
order. 

The  smaller  blank  areas  result  from  the  effect  of  stratification  on  the  velocity 
profile  (Section  2.3).  The  neutral  and  stratified  cases  in  each  figure  represent  the 
same  range  of  bottom  shear  stress.  Because  the  stratification  dampens  the  transfer 
of  low-momentum  fluid  up  from  the  bottom,  however,  the  predicted  velocities  at 
one  meter  above  the  bottom  are  higher  in  the  stratified  case  than  in  the  neutral 
case.  This  effect  is  demonstrated  in  the  distribution  of  the  ‘data’  points  (xs)  in  the 
plots.  The  magnitude  of  the  effect  of  stratification  on  the  velocity  can  be  gauged 
by  comparing  the  location  of  the  corresponding  neutral  and  stratified  data  points. 
Each  of  the  circled  points  in  Figure  5.2,  plots  (a)  and  (b)  represents  the  same  u.^ 
and  u.ev,  values  for  this  wave  and  reference  concentration  (C(6«,)  ~  5  x  10“^)  for 
fine  sand.  The  predicted  neutral  current  velocity,  however,  is  ~  40^  ~  1.0), 

and  the  predicted  stratified  current  velocity  ~  50  ~  1.25).  The  difference 

of  10  cm/sec  represents  the  value  of  the  stratification  term  in  the  velocity  profile 
calculation  (Equation  3.44). 

A  decrease  in  predicted  sediment  load  and  transport  is  expected,  since  the 
upward  transfer  of  mass,  like  the  upward  transfer  of  momentum,  is  damped  by 
stratification.  The  magnitude  of  the  response  to  stratification  is  strongly  depen¬ 
dent  on  the  relative  wave-current  and  current  shear  stresses.  The  response  of  the 
sediment  load  and  transport  predictions  to  stratification  will  be  discussed  in  three 
categories:  1)  small  concentrations  at  the  top  of  the  wave  boundary  layer;  2)  low 
current  velocity;  and  3)  high  reference  concentrations  with  large  currents.  Within 
these  categories,  the  sediment  response  hinges  on  the  balance  of  terms  in  the  sim- 
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plified  sediment  continuity  equation  (Equation  3.34).  For  the  sake  of  reference  in 
the  discussion,  the  equation  is  repeated  here: 


WfC  +  -  0 

Note  that  the  product  of  the  stratification-corrected  eddy  viscosity  {yt$)  and  the 
concentration  gradient  is  balanced  by  the  product  of  the  fall  velocity  and  the  con¬ 
centration. 

Small  C(S^),  Large  current  velocity 

Under  this  relatively  strong  wave  forcing,  when  the  sediment  concentration  at  the 
top  of  the  wave  boundary  layer  is  small  (0(10“^)  and  there  is  a  strong  current, 
the  response  to  stratification  is  minor.  These  conditions  are  demonstrated  in  Fig¬ 
ures  5. 2-5. 5  by  the  range  of  contours  near  the  left  hand  margin  of  each  plot.  The 
highest  values  vary  only  slightly  between  the  stratified  and  neutral  cases.  For  ex¬ 
ample,  the  upper  left  hand  data  points  in  Figure  5.2(a)  and  (c)  are  located  at  the 
vertical  coordinate  1.15  and  have  load  and  transport  values  of  0.04  cm®/cm*  and 
2.0  cm^ ! cm/ sec  respectively;  the  corresponding  point  for  plots  (b)  and  (d)  are  at 
the  vertical  coordinate  1.2  and  have  values  of  about  0.028  and  1.7.  Stratification 
does  not  substantially  alter  the  load,  transport,  or  velocity  value  for  this  fine  sand 
case,  and  the  same  is  true  for  the  other  three  grain  sizes.  In  terms  of  the  sediment 
continuity  equation,  the  small  concentration  and  relatively  large  eddy  viscosity 
(since  u.c  is  relatively  large)  mean  that  the  concentration  gradient  is  relatively 
small,  regardless  of  the  fall  velocity.  The  stability  parameter,  which  lowers  the 
value  of  the  eddy  viscosity,  and  is  a  function  of  the  concentration  gradient,  stays 
fairly  small. 

Physically,  this  corresponds  to  a  situation  where  the  currents  are  strong  enough 
to  suspend  enough  of  the  small  boundary  concentration  to  preclude  a  strong  con- 


centration  gradient.  In  the  calculation  procedure,  the  dependence  on  C  (^,„)  can  be 
seen  in  the  stratification  terms  of  the  solutions  for  the  velocity  and  concentration 
profiles  (Equations  3.44  and  3.45).  Both  depend  on  the  integral  of  the  concentra¬ 
tion.  A  small  boundary  concentration  therefore  leads  to  small  stratification  effects. 

Low  current  velocity 

A  small  current  relative  to  the  wave  results  in  strong  stratification  effects  at  the  top 
of  the  wave  boundary  layer.  Load  and  transport  may  drop  by  an  order  of  magnitude 
or  more.  The  effect  is  intensified  with  increasing  C[6y,)  and  with  decreasing  grain 
size.  In  Figure  5.2(a]  and  (b),  for  ^  =  0.4  (u,  =  16^),  stratification  leads  to  a 
decrease  in  load  of  approximately  50%  when  C{6J)  =  10“^.  The  values  on  which 
these  comparisons  are  based  are  marked  by  boxes  on  the  plots  and  enumerated  in 
Table  5.3(a)  and  (b).  Moving  to  the  right  across  the  plots  to  the  points  marked  at 
C(Sn)  =  10”’,  the  decrease  in  load  due  to  stratification  is  about  70%.  At  the  same 
current  level  (^  =  0.4),  the  decrease  in  load  for  very  fine  sand  ranges  from  67% 
to  90%  for  C(6v)  =  10“^  - 10"’  (marked  on  Figure  5.3(a)  and  (b)).  For  coarse  silt, 
the  same  range  of  concentrations  at  the  samae  velocity  level  yields  decreases  in  load 
of  83-96%  (marked  on  Figure  5.4). 

The  effect  on  transport  calculations  is  even  more  pronounced  since,  in  stratified 
flow,  the  sediment  that  madces  up  the  load  is  concentrated  near  the  bottom  where 
velocities  are  lowest.  In  the  same  four  figures  (mau’ked  on  Figures  5. 2-5.5,  plots 
(c)  and  (d))  at  the  same  current  level  (^  =  0>4)  for  the  same  range  of  C(Sy,) 
(10"’  —  10"^),  the  transport  of  fine  sand  drops  by  83-97%  due  to  stratification; 
transport  of  very  fine  sand  drops  by  89-98%;  and  transport  of  silt  drops  by  90-99% 
(again,  compare  top  two  sections  of  Table  5.3). 

The  first  order  effect  of  stratification  in  this  large  wave  environment  where 
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a) Low  current,  low  C{S,u) 

0.4  10-*  0.02 

0.01 

0.005 

0.3 

0.05 

percentage  change 

-50% 

-83% 

0.4  lO-*  0.01 

percentage  change 

0.03 

-67% 

0.009 

0.9 

-89% 

0.1 

0.4  10-*  0.006 

percentage  change 

0.09 

-83% 

0.015 

2.0 

-90% 

0.2 

0.4  10-*  0.0025 

percentage  change 

0.27 

-90% 

0.026 

5.5 

-93% 

0.37 

\))Low  current,  high  C{6yj) 

0.4  10-®  0.02 

percentage  change 

0.11 

-70% 

0.034 

1.8 

-97% 

0.05 

lol  l0^^5  0.01 

percentage  change 

0.25 

-90% 

0.025 

5.0 

-98% 

0.08 

0.4  10-®  0.006 

percentage  change 

0.45 

-96% 

0.018 

10.0 

-99% 

0.08 

"Ol  10^^^  0.0025 

percentage  change 

1.33 

-99% 

0.018 

27.4 

-99% 

0.06 

c)High  current,  high  C(6„) 

1.0  10"®  0.02 
percentage  change 

0.28 

-77% 

0.065 

16 

-89% 

1.8 

To  0.01 

percentage  change 

0.70 

-88% 

0.08 

41 

-92% 

3.5 

To  10^  0.006 

percentage  change 

1.2 

85% 

0.18 

60 

-83% 

10 

percentage  change 


300.0 


-20% 


Table  5.3:  Values  of  load  and  transport  in  a  moderate  storm,  picked  from  Figures  5. 2-5.5 
for  ranges  of  conditions  as  discussed  in  text.  For  each  current  category,  results  are  shown 
for  (top-to-bottom)  fine  sand,  very  fine  sand,  and  coarse  silt.  Magnitudes  of  decreases  due 
to  stratification  are  shown.  Load  units  are  transport  units  are  cm^/cm/aec. 


the  current  is  small  is  such  a  drastic  decrease  in  predicted  load  and  transport 
that  the  second-order  effects  of  grain  size  and  reference  concentration  might  seem 
gratuitous.  For  example,  there  appears  to  be  little  difference  between  97%,  98%, 
and  99%  decreases  in  transport  for  fine  sand,  very  fine  sand,  and  silt  in  low  current, 
high  conditions  (Table  5.3(b)).  The  difference  is  made  obvious,  however,  by 

observing  that  instead  of  an  order  of  magnitude  increase  in  transport  from  fine 
sand  to  silt,  transport  predictions  for  all  three  grain  sizes  in  a  16^  current  with 
relatively  high  reference  concentration  are  nearly  the  same.  For  the  same  wave, 
current  and  reference  concentration  conditons,  the  predicted  load  for  silts  is  a  factor 
of  two  smaller  than  the  predicted  load  for  fine  sand.  Thus  stratification  may  nullify 
or  reverse  the  expected  result  that,  for  a  given  wave  and  current  condition,  finer 
sediment  leads  to  greater  suspended  load  and  transport. 

This  counter-intuitive  result  is  caused  by  the  feedback  between  stratification 
2md  shear  velocity.  If  the  holding  capacity  for  suspended  sediment  of  the  current¬ 
generated  turbulence  is  much  less  than  that  due  to  the  wave-current  enhzmced  shear 
stress  in  the  wave  boundary  layer,  then  the  stratification  due  to  the  suspended 
sediment  concentration  gradient  at  the  top  of  the  wave  boundary  layer  acts  to 
damp  the  mean  shear  stress  (but  not  the  maximum  shear  stress,  because  it  is 
dependent  primarily  on  the  wave),  thereby  enhancing  the  stratification.  The  effect 
increases  with  smaller  grains  (for  the  same  wave  and  current  conditions)  and  with 
increased  reference  concentration  (for  the  same  grain  size)  for  the  same  reason: 
the  stratification  terms  depend  on  the  integral  of  the  concentration,  and  a  large 
boundary  concentration  introduces  the  possibility  of  a  large  concentration  gradient. 

In  terms  of  the  sediment  continuity  equation,  this  translates  to  a  relatively  small 
eddy  viscosity  (since  u.e  is  small)  requiring  a  relatively  large  concentration  gradient 
to  balance  the  product  of  fall  velocity  with  concentration.  The  larger  gradient 
generates  a.  larger  stability  parameter,  which  leads  to  a  smaller  eddy  viscosity,  and 


Figure  5.10:  Schematic  diagram  of  alternative  solutions  to  stratified  velocity  profile  for 
neru-bottom  wave  and  current  boundary  layer  flow.  Solid  line  is  low  current  shear  stress 
case  with  strong  stratification;  dashed  line  is  high  current  shear  stress  case. 


therefore  the  gradient  must  be  still  larger.  This  is  the  reason  for  the  drastic  drop 
in  expected  load  and  transport  when  the  current  is  small. 

The  medium  silt  results  present  a  special  case  for  the  stratified  results,  as  a  look 
at  the  contours  for  the  stratified  case  suggests  (Figure  5.5  (b)  and  (d)).  The  small 
fall  velocity  and  large  concentration  require  a  large  sediment  gradient  when  the 
current  eddy  viscosity  is  small,  in  order  to  balamce  terms  in  the  sediment  continu¬ 
ity  equation.  The  current  eddy  viscosity  decreases  with  an  increased  concentration 


gradient,  leading  to  a  stronger  gradient  to  balance  the  sediment  continuity  equa¬ 
tion,  and  so  on  until  convergence. 

The  velocity  profile  is  likewise  affected  by  the  enhanced  stratification.  The 
velocity  profile  is  the  sum  of  a  stratification  term  amd  a  logarithmic  term  (Equar 
tion  3.44).  If  the  stratification  term  increases  more  rapidly  in  the  strongly  stratified 
case  than  the  logarithmic  term  does  in  the  less-strongly  stratified  case,  then  the 
velocity  profiles  can  overlap  at  some  heights  above  the  bottom  (Figure  5.10).  If 
the  reference  velocity  is  specified  in  the  vicinity  of  these  heights,  then  there  is  more 
than  one  solution,  amd  predictions  of  load  and  transport  for  the  solutions  can  vary 
a  lot. 

While  the  results  for  relatively  low  reference  concentrations  are  similar  to  those 
for  larger  grain  sizes,  stratification  effects  are  increased  when  grains  are  small,  the 
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(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

0.175 

0.15 

0.125 

0.10 

0.075 

0.05 

4.57 

4.47 

4.37 

4.27 

4.17 

4.06 

1.84 

1.70 

1.55 

1.39 

1.20 

0.97 

35.2 

34.7 

34.4 

34.3 

34.5 

34.2 

22.1 

20.1 

17.9 

15.5 

12.8 

10.0 

13.1 

14.6 

16.5 

18.8 

21.6 

24.2 

0.97 

0.79 

0.63 

0.47 

0.32 

0.18 

70.7 

54.5 

39.5 

26.2 

15.0 

6.5 

Ur,  str.  contrib. 
N.B.  Load 
N.B.  Trans. 


Table  5.4;  Six  sets  of  model  parameters  representing  the  six  load  and  transport  predic¬ 
tions  in  the  boxed  region  of  Figure  5.5(b)  and  (d).  The  matching  parameter  for  model 
convergence,  Ur,  is  effectively  the  same  for  all  six  cases,  but  the  near-bottom  (N.B.)  load 
and  transport  vary  by  factors  of  5  and  10,  respectively.  The  model  solution  is  not  uniquely 
determined  in  this  case. 


concentration  at  the  top  of  the  wave  boundary  layer  is  greater  than  10~®  and  the 
mean  shear  stress  is  low.  To  illustrate  this  point  a  detailed  examination  of  the 
conditions  defining  the  boxed  region  in  Figure  5.5(b)  will  be  undertaken.  The  six 
points  in  the  box  represent  parameter  ranges  as  listed  in  Table  5.4. 

The  values  listed  in  Table  5.4  for  the  reference  velocity  Ur  calculated  by  the 
model  are  nearly  identical:  they  vary  by  less  than  the  resolution  of  field  measure¬ 
ments,  certainly,  and  are  within  the  one  per  cent  convergence  criterion  for  the 
model.  The  load  and  transport  predictions,  however,  vary  by  factors  of  five  to 
eleven  over  this  essentially  identically-specified  range  of  solutions.  The  difference 
in  Utcw  values  is  very  small,  but  u.^  values  double  over  the  range.  The  reason 
for  this  ambiguity  is  found  in  lines  5  and  6  of  Table  5.4.  Note  that  the  stratifi¬ 
cation  contribution  to  the  reference  velocity  increatses  at  the  same  rate  that  the 
logarithmic  contribution  decreeises  as  u,c  goes  down.  The  stability  parameter  at  Zr 
combines  with  the  variation  in  u.^  to  give  this  ambiguous  result. 

The  similarity  in  predicted  velocities  leads  to  the  difficulty  in  applying  the 
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model  because  convergence  is  set  on  Ur-  As  Table  5.4  demonstrates,  the  same  u, 
can  be  obtained  from  very  different  mean  shear  stress  values.  The  predicted  velocity 
profiles  for  points  (3)  and  (5)  of  Table  5.4  are  shown  in  Figure  5.11,  plotted  linearly 
to  show  detail.  Obviously,  the  model  could  converge  on  either  of  these  solutions  if 
Zr  were  specified  between  about  0.2  and  1.4  meters,  and  the  difference  in  predicted 
load  and  transport  is  a  factor  of  two.  Similar  plots  could  be  shown  for  the  other 
four  points,  so  that  the  possible  range  in  predicted  load  is  a  factor  of  five,  and 
in  transport  a  factor  of  eleven.  To  choose  between  these  options,  a  velocity  or 
concentration  measurement  somewhere  higher  in  the  water  column  would  have  to 
be  specified.  More  fundamentally,  this  result  calls  into  question  the  representation 
of  stratification  using  a  stability  parameter  to  modify  the  eddy  viscosity,  in  some 
circumstances.  It  is  not  clear  that  the  stability  parauneter  should  be  allowed  to 
grow  so  large,  or  that  it  is  an  accurate  representation  of  the  physical  world  when  it 
does.  This  set  of  wave-current-sediment  predictions  needs  to  be  tested  in  controlled 
field  or  laboratory  conditions. 

Large  C(6„),  large  current 

In  contrast  to  the  high  reference  concentration  and  low  current  scenario  just  dis¬ 
cussed,  a  higher  current  restores  the  expected  pattern  of  increased  load  and  trans¬ 
port  with  decreasing  grain  size  when  observed  wave  and  current  conditons  are  the 
same.  The  predictions  used  to  demonstrate  this  point  are  those  for  a  current  of 
40^  ~  1)>  marked  by  boxes  on  Figures  5.2-5.5  amd  listed  in  Table  5.3(c). 

Again,  the  first  order  effect  on  load  and  transport  values  is  stratification:  the 
stratified  load  predictions  are  fifty-two  to  eighty-eight  per  cent  lower  than  the 
neutral  ones.  In  looking  at  the  percentage  decreases,  however,  note  that  the  largest 
change  in  the  stratification  effect,  with  respect  to  the  low  current  case,  occurs  in 
the  fine  grain  size  cases.  For  medium  silt,  instead  of  a  ninety-six  to  ninety-nine  per 
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Figure  5.11:  Velocity  profiles,  plotted  linearly  to  show  detail,  illustrating  how 
stratification  effects  can  lead  to  non-unique  model  solutions.  Specification  of  the 
reference  velocity  between  0.2  amd  1.4  m  could  yield  either  of  these  velocity  profile 
solutions.  These  profiles  represent  velocity  profile  predictions  (3)  and  (5)  for  the 
wave-current-velocity  results  shown  in  Table  5.4. 


cent  decrease  in  load  and  transport  due  to  stratification,  there  is  an  eighty-three  to 
eighty-five  per  cent  decrease,  a  change  of  eleven  to  sixteen  per  cent.  The  increase 
in  load  and  transport  relative  to  the  neutral  case  for  medium  silt  is  even  greater. 

What  these  changes  mean  is  that  the  load  increases  with  decreasing  grain  size, 
as  expected.  Even  more  pronounced  is  the  change  in  transport  patterns.  Rather 
than  being  practically  the  same,  transport  of  coarse  silt  is  an  order  of  magnitude 
greater  than  transport  of  fine  sauid  (1.8  vs.  10  cm*/cm/sec)  and  transport  of 
medium  silt  is  an  order  of  magnitude  higher  still  (241  cm^ / cm f  sec). 

The  shapes  of  the  contours  in  the  plots  of  stratified  predictions  result  from  the 
reconciliation  of  the  tendency  of  the  stratification  effects  to  grow  with  increasing 
C{6^)  and  to  shrink  with  increasing  current.  The  regular,  hyperbolic  shapes  of  the 
neutral  contours  are  suggested  in  the  small-concentration  regions  at  the  left-hand 
edge  of  some  of  the  plots.  The  curves  begin  to  flatten  in  all  cases,  since  stratification 
makes  the  flow  less  competent  to  hold  sediment.  In  some  cases  the  curves  remain 
flat,  indicating  that  the  effects  of  increasing  reference  concentration  and  increasing 
velocity  balance.  For  smaller  grain  sizes,  the  curves  rise  to  the  right,  suggesting 
that  an  increase  in  available  sediment  leads  to  a  decrease  in  total  load  or  transport. 
As  discussed  earlier,  this  results  from  the  larger  concentration  gradient  at  the  top 
of  the  wave  boundary  layer. 

Grain  size  effects 

The  influence  of  grain  size  on  stratification  effects  was  discussed  above.  As  a  gen¬ 
eral  rule,  smaller  grain  sizes  enhance  the  effect  of  stratification  due  to  the  gradient 
in  eddy  viscosity  at  when  the  wave  is  large  and  the  current  small.  Grain  size  also 
influences  the  value  of  C{Su)  in  three  ways:  1)  the  ^  ratio  determines  the  con¬ 
centration  gradient  (Equation  3.35);  2)  the  critical  Shields  parameter  for  initiation 
of  motion  is  a  factor  in  the  reference  concentration;  and  3)  the  bottom  rough- 
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ness  is  determined  by  the  grain  size,  whether  through  the  size  of  ripples,  sediment 
transport  layer  or  grain  roughness  itself  (Section  3.3).  This  is  very  important  in 
the  calculation  of  boundary  shear  stress.  Ultimately,  grain  size  is  the  determining 
factor  in  the  quantity  of  sediment  that  can  be  carried  by  a  flow. 

The  variation  of  the  concentration  at  the  top  of  the  wave  boundary  layer  with 
grain  size  is  most  easily  seen  in  the  x-axes  of  Figures  5.2-5. 5.  The  x-axes’  maxima, 
in  decreasing  order  of  grain  size,  are:  10“*,  1.5  x  10"*,  3  x  10"*  and  6  x  10"*.  Since 
all  of  these  plots  represent  approximately  the  same  ranges  of  wave  and  current 
velocities,  the  variation  in  C{6^)  results  entirely  from  initiation  of  motion  criteria, 
movable  bed  roughness,  and  vertical  decay  inside  the  wave  boundary  layer.  In 
these  cases,  where  wave  shear  stress  is  very  large,  the  differences  are  due  primarily 
to  initiation  of  motion  criteria  and  roughness. 

The  effect  of  roughness  on  C{6J)  for  a  single  grain  size  is  demonstrated  in  the 
plots.  If  one  ignores  the  contours,  the  data  points  (xs)  represent  an  x-y  plot  of  the 
variation  of  C(Syj)  with  reference  velocity.  If  roughness  had  no  effect  on  u,**,,  and 
therefore  on  C(6,„),  the  data  points  in  each  plot  would  be  a  set  of  six  near- vertical 
lines,  curving  slightly  to  the  right  in  response  to  the  slight  increase  in  boundary 
shear  stress  (ro.ew)  due  to  incremental  increases  in  current  contribution.  However, 
with  increases  in  current  and  shear  stress,  the  bottom  roughness  either  increases 
because  bedforms  and  sediment  transport  layers  grow,  or  decreases  because  rip¬ 
ples  are  washed  out.  If  the  roughness  increases,  then  the  shear  stress  in  the  wave 
boundary  layer  increases  and  the  concentration  at  the  top  of  the  wave  boundary 
layer  increases.  The  data  points  in  Figures  5.3-5.5  show  increases  in  C(5«,)  with 
increases  in  current  (i.e.  they  curve  to  the  right).  This  reflects  the  bedform  rough¬ 
ness  element  due  to  the  building  of  ripples  in  the  sand  ceise,  and  the  sediment 
trzmsport  roughness  (kt^j)  due  to  formation  of  a  transport  layer  in  the  silts.  The 
decrease  in  C(6«,)  with  increasing  current  in  the  fine  sand  case,  demonstrated  in 


the  data  points’  curve  to  the  left  in  Figure  5.2,  indicates  that  ripples  are  being 
washed  away  because  of  the  increasing  bottom  shear  stress. 

These  effects  demonstrate  the  interconnections  between  waves,  currents,  grain 
diameter,  bottom  roughness,  and  concentration.  Since  the  predictions  of  stratified 
load  and  transport  for  a  given  grain  size  depend  only  weakly  on  changes  in  C(^,u), 
its  variation  is  more  significant  between  grain  sizes  thzm  within  one  grain  size  when 
the  wave  is  large. 

Bedload  and  Outer  Ekman  Layer  Load  and  Transport 

In  a  storm  wave  environment  the  effects  of  stratification  on  the  values  of  the  outer 
Ekman  layer  load  and  transport  are  more  pronounced  than  the  stratification  effects 
on  bedload  transport  in  a  storm  wave  environment.  The  proportion  of  bedload 
transport  compared  with  near-bottom  transport  in  the  stratified  storm  wave  case 
is  similar  to  that  in  the  neutral  case  (Figures  5.6(b)-5.9(b)).  The  smaller  total 
transport  values  make  the  bedload  proportion  slightly  higher,  but  the  maximum 
is  still  only  ten  per  cent  of  the  near-bottom  transport  (Figure  5.6(b)).  In  general, 
bedload  transport  is  less  than  five  per  cent  of  the  stratified  near-bottom  suspended 
transport. 

Stratification  makes  a  big  difference  in  the  predicted  outer  Ekman  layer  load 
and  transport.  Most  important  for  the  application  of  this  model,  the  predicted 
transports  are,  in  general,  less  than  or  equal  to  the  near-bottom  transports  even 
for  the  finer-grained  sediments  (Figures  5.6(d)-5.9(d)).  Since  the  focus  of  this  work 
is  the  near-bottom  trzmsport,  and  the  predictions  for  the  region  farther  from  the 
bottom  are  more  poorly  constrained  (Section  3.2.4),  this  pattern  is  reassuring  with 
respect  to  the  usefulness  of  the  present  results. 

In  addition  to  the  changes  in  the  magnitudes  of  the  contours,  the  patterns  of  the 
contours  are  different  for  the  stratified  vs.  neutral  outer  Ekman  to  near  bottom 


transport.  The  ratio  increases  with  increasing  current  in  the  stratified  case  as 
it  does  in  the  neutral  case.  However,  since  the  stratified  concentration  profile  is  a 
nonlinear  function  of  the  concentration  at  the  top  of  the  wave  boundary  layer,  with 
effects  of  stratification  increasing  as  increcises,  the  proportion  of  sediment 

transport  in  the  outer  Ekman  layer  decreases  with  increasing  C{6^). 

What  this  means  is  that  for  large  waves,  the  transport  in  the  outer  Ekman 
layer  is  relatively  unimportant  except  when  both  large  currents  and  small  grains 
are  present.  These  cases  must  be  approached  with  caution  for  a  number  of  other 
reasons,  as  discussed  in  relation  to  the  neutral  case;  time-dependence  of  the  model, 
merging  of  the  surface  and  bottom  Ekman  layers,  armoring  of  mixed  grain  sizes, 
and  variations  of  fall  velocity  and  critical  Shields  parameter  due  to  adhesion  and 
cohesion.  However,  the  unrealistically  high  sediment  concentrations  above  the  wave 
boundary  layer  characteristic  of  the  neutral  case  are  not  found  in  the  stratified 
results.  The  inclusion  of  stratification  thus  provides  one  physically  reasonable 
approach  to  eliminating  unsatisfactory  results  from  the  neutral  model.  For  most 
cases,  the  proportion  of  the  total  load  in  the  outer  Ekman  layer  is  small  enough 
that  error  in  the  method  used  to  calculate  the  estimate  of  its  volume  is  acceptably 
low. 

Total  load  and  transport 

The  discussion  of  both  neutral  and  stratified  results  so  far  has  focused  on  the 
relative  importance  of  various  parameters  and  conditions  to  the  calculation  of  sed¬ 
iment  transport  values.  The  amount  of  sediment  involved  is  most  easily  grzisped 
by  converting  the  volumetric  load  values  into  estimates  of  the  reworking  depth,  by 
dividing  the  load  by  the  bed  concentration,  here  taken  to  be  0.6. 

The  reworking  depth  is  very  grain-size  dependent.  For  the  neutral,  near-bottom 
load,  each  step  up  in  grain  size  corresponds  to  a  decrease  in  load  and  transport 


Ranges  of  reworking  depths.  Moderate  Storm 
Near-bottom  load  only 

d{cm) 

neutral  (cm) 

stratified  (cm) 

0.02 

0.01  -  0.53 

0.007  -  0.175 

0.01 

0.02  -  1.6 

0.005  -  0.35 

0.006 

0.055  -  4.6 

0.008  -  0.97 

0.0025 

0.22  -  14.2 

0.03  -  7.6 

Table  5.5:  Reworking  depth  ranges  for  the  moderate  storm  wave  showing  how  stratifica¬ 
tion  effects  are  amplified  with  decreasing  grain  size.  Note  that  the  neutral  and  stratified 
predictions  for  the  same  grain  sizes  represent  the  same  range  of  wave-current  shear  stress 
values,  but  not  the  same  range  of  reference  current  values. 

by  a  factor  of  ~  2  —  3.  Eight  to  10  times  as  much  coarse  silt  as  fine  sand  would, 
therefore,  be  lifted  from  the  seafloor  and  transported  under  these  neutral  wave  and 
current  conditions  (Table  5.5). 

In  the  stratified  cases,  the  erosion  and  transport  values  are  lower,  but  the 
pattern  of  increase,  at  least  for  the  high  velocity  and  high  C((5«,)  cases,  is  similar 
to  the  neutral.  For  each  increase  in  grain  size,  the  values  of  load  and  transport 
predictions  drop  by  a  factor  of  —  1.5—2  (Figures  5.2-5. 5).  For  the  low  velocity  cases, 
however,  the  stratification  effects  are  stronger  for  smaller  grain  sizes,  so  there  may 
be  little  or  no  change  in  predicted  load  due  to  grain  size.  The  ranges  in  predicted 
reworking  depths  are  approximately  halved  relative  to  the  neutral  case  (Table  5.5), 
but  higher  reference  current  velocities  are  needed  to  generate  the  smaller  loads. 

5.3.3  Summary 

While  a  wave  with  a  bottom  velocity  of  40  cm/sec  can  put  much  more  sediment  in 
motion  than  could  a  current  alone,  the  stratification  generated  by  the  gradient  in 
turbulent  energy  beyond  the  wave  boundary  layer  tempers  drastically  the  erosive 
potential  of  the  combined  wave-current  force.  The  load  and  transport  are  depen- 


dent  on  grain  size,  current  velocity  and  reference  concentration,  but  the  maximum 
expected  erosion  due  to  a  moderate  storm  appears  to  be  on  the  order  of  1-5  cm, 
even  in  very  strong  currents.  Neutral  predictions  for  load  and  transport  under 
large  waves  and  currents  appear  to  overpredict  concentrations  by  at  least  an  order 
of  magnitude.  This  effect  is  countered  by  stratification,  which  may  be  an  important 
factor  in  limiting  transport. 

5.4  Results:  Low  Swell 

On  some  continental  shelves,  low,  long  period  swell  waves  occur  as  forerunners  of  an 
approaching  storm  or  as  remnants  of  a  distant  storm.  Swell  is  characterized  by  its 
long  period,  monochromatic  character.  Periods  can  be  up  to  25  or  30  seconds.  This 
set  of  model  runs  treats  a  20  second  wave  generating  bottom  velocities  of  20  cm/sec. 
Such  a  wave  would  have  a  height  of  1.1  m  in  50  m  water  depth.  Figures  5.12-5.15 
show  the  variation  in  load  and  transport  predictions  with  C[S.u,)  and  in  the 
same  format  as  discussed  in  detail  at  the  end  of  Section  5.2. 

These  contours  are  much  simpler  than  those  for  the  cases  with  larger  wave 
bottom  velocities.  The  most  striking  feature  of  these  plots,  though,  is  the  blank 
space  in  the  lower  half  of  each  plot.  Before  examining  the  load  and  transport 
results,  this  phenomenon  will  be  explained. 

Unlike  the  last  wave  cztse,  the  small  wave  considered  here  does  not  generate 
sufficient  boundary  shear  stress  to  put  sediment  into  motion  except  in  the  finest 
sediment  class  considered  (medium  silt).  In  the  fine  sand  czise,  for  example  (Fig¬ 
ure  5.12),  only  when  ^  is  approximately  equal  to  one  does  sediment  begin  to  move. 
At  lower  velocities,  the  data  points  plot  onto  the  y-axis  to  signify  a  reference  con¬ 
centration  of  zero.  The  blank  areas  on  the  plots,  then,  signal  conditions  under 
which  no  motion  is  predicted:  for  example,  on  a  seabed  of  very  fine  sand  in  low 
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Figure  5.12:  Predicted  near-bottom  load  (^r)  »nd  transport  (cm*/cm/sec)  values  for 
fine  sand  in  low  swell  waves.  T=20  sec,  and  the  wave  height  in  50  m  of  water  would  be 
1.1m.  Ur  is  reference  current  velocity.  ranges  from  3-21  cm;  ^  ranges  from  1  to  28  m. 
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Figure  5.15:  Predicted  near-bottom  load  (|^)  and  transport  {cm^/ cm  fate)  values  for 
medium  silt  in  low  swell  waves,  fi*  ranges  from  7.7-8. 1  cm;  ^  ranges  from  3  to  10  m. 


-1.  . N*.% 


swell  of  this  sort,  with  a  current  of  10  cm/sec  =  0.5),  one  would  not  expect  to 
find  C{6n)  =  10~^  (Figure  5.13). 

When  the  moveable  bed  roughness  model  is  used,  the  only  roughness  element 
under  conditions  of  no  motion  is  the  grain  size  roughness  Such  small  rough¬ 
nesses  generate  very  little  turbulence,  so  the  model  predicts  small  Utew  a'Ud  u»e 
values  and,  therefore,  very  low  current  values.  Once  the  skin  friction  shear  stress 
exceeds  the  critical  value,  however,  the  movable  bed  roughness  model  generates 
roughness  estimates  based  on  equilibrium  bedform  and  sediment  transport  mod¬ 
els,  causing  a  discontinuity  between  the  grain  roughness  when  there  is  no  sediment 
motion  and  the  bedform  roughness  which  is  as  much  ;as  an  order  of  magnitude  or 
more  greater. 

The  roughness  discontinuity  leads  to  a  discontinuity  in  predicted  reference  ve¬ 
locity:  the  increased  roughness  leads  to  increased  turbulence  and  a  jump  in  the 
predicted  velocity.  Note,  for  example,  that  ^  =  0.45  in  Figure  5.12a  is  the  largest 
velocity  at  which  there  is  no  sediment  motion  (data  point  plots  on  the  y-axis).  The 
non-dimensionalized  reference  velocity  jumps  to  ^  =  0.98  (the  next  highest  data 
point)  with  the  next  incremental  increase  in  ^  because  sediment  motion  has  begun 
and  the  moveable  bed  roughness  model  is  predicting  ripple  formation.  This  large 
change  in  predicted  reference  velocity,  however,  represents  only  an  incremental 
change  in  the  skin  friction  component  of  bottom  shesu-  stress.  The  ripples  increaise 
the  turbulent  intensity,  so  that  the  current  feels  more  drag  (u.c  is  greater),  and  the 
transport  capacity  jumps  dramatically. 

This  signals  a  logistical  problem  for  the  model:  what  if  the  input  condition 
it  is  trying  to  match  is  a  reference  velocity  of  15^  (^  =  0.75)  when  waves  and 
currents  are  as  shown  in  Figure  5.12.  The  model’s  iterations  would  skip  back  and 
forth  between  a  no-motion  case  (which  would  give  a  velocity  too  low)  and  a  rippled 
bed  ceise,  which  gives  a  velocity  that  is  too  high.  There  is  no  value  of  ^  that  would 
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give  the  specified  reference  velocity. 

There  is  a  simple  solution:  include  a  bed  roughness  as  input,  based  on  observed 
or  expected  animal  mounds  and  ripples,  rather  than  using  the  moveable  bed  rough¬ 
ness  model.  This  is  a  realistic  solution,  since  the  ocean  floor  would  not  be  expected 
to  be  flat,  zmd  the  small  wave-current  velocities  under  which  this  problem  occurs 
would  not  likely  be  the  factors  controlling  bed  roughness. 

Like  the  storm  wave  cases,  the  effects  of  stratification  grow  pronounced  with 
decrease  in  grain  size  and  with  increase  in  concentration.  Because  the  model  pre¬ 
dicts  no  motion  at  small  current  velocities  for  all  but  the  medium  silt,  it  is  the 
only  case  to  display  the  effects  of  the  gradient  in  turbulent  energy  at  the  top  of 
the  wave  boundary  layer.  Note  the  lowest  contours  flattening  toward  horizontal 
(Figure  5.15);  the  effect  of  stratification  on  load  and  transport  of  medium  silt  is  a 
decrease  of  as  much  as  80%.  Since  the  total  transport  is  so  small,  however,  that 
translates  to  a  maximum  difference  in  reworking  depth  of  only  2  mm  even  for  silt. 

In  general,  the  concentrations  at  the  top  of  the  wave  boundary  layer  are  very 
small  compared  with  those  generated  by  the  large  wave  examined  above.  These 
low  concentration  and  the  low  wave-current  shear  stresses  lead  to  small  differences 
between  neutral  and  stratified  predictions,  generally  30%  or  less.  In  all  ceises,  the 
difference  is  beyond  the  resolution  of  field  measurements. 

The  outer  Ekman  layer  transport  is  negligible  in  fine  sand  applications,  reflect¬ 
ing  the  reduced  mobility  of  the  sediment  under  the  small  wave.  Silt  and  very  fine 
sand  have  outer  Ekman  layer  transport  approximately  equal  to  the  near-bottom 
values  for  all  currents  greater  than  20  cm/sec.  Bedload  transport  remains  less  than 
10%  of  near-bottom  transport  in  all  silt  cases,  but  forms  most  of  the  (very  small) 
sand  transport  in  currents  smaller  than  about  20  cm/sec. 

The  ranges  of  reworking  depths  resulting  from  the  near-bottom  load  for  this 
wave  are,  not  surprisingly,  much  smaller  than  for  the  larger  wave  (Table  5.6.  In- 


Ranges  of  reworking  depths,  Low  swell 
Near-bottom  load  only 


d{cm) 

neutral  (cm) 

stratified  (cm) 

0.02 

0.0  -  0.01 

0.0- 

0.008 

0.01 

0.0  -  0.10 

0.0- 

0.075 

0.006 

0.0  -  0.35 

0.0- 

0.23 

0.0025 

0.005  -  0.87 

0.003 

-  0.57 

Table  5.6:  Reworking  depths  for  low  swell,  showing  that  stratification  reduces  the 
near-bottom  load  by  less  than  a  factor  of  two  in  all  cases. 

eluding  the  outer  Ekman  layer  load  increases  the  meocimum  reworking  depth  to 
0.01  cm  for  fine  sand,  to  0.1  cm  for  very  fine  sand,  to  0.6  cm  for  coarse  silt,  and  to 
1.1  cm  for  medium  silt.  Maximum  reference  currents  are  about  30  cm/sec;  higher 
currents  would  of  course  result  in  higher  loads. 

In  summary,  low  swell  in  the  presence  of  only  small  currents  is  relatively  inef¬ 
fective  as  a  means  of  sediment  transport.  If  only  medium  and  coarser  sand-sized 
grains  are  present,  transport  is  likely  to  be  limited  to  bedload.  Silt-sized  grains 
may  be  transported  in  moderately  large  quantities  if  currents  are  strong. 

5.5  Results:  Large  swell 

Occasionally  on  the  continental  shelf,  the  remnants  of  a  large,  distant  storm  will 
bring  high,  long  period  waves  that  generate  substantial  bottom  velocities  at  all 
shelf  depths.  The  case  examined  here  is  a  20  second  wave  generating  bottom 
velocities  of  60  cm/sec;  such  a  wave  would  be  3.3  m  high  at  50  m  water  depth 
(Figures  5.16-5.19). 

The  most  notable  thing  about  these  plots  is  their  striking  similarity  to  the 
corresponding  plots  for  the  moderate  storm  wave  case  (Figures  5.2-5. 5).  Since 
all  stratification  and  grain  size  effects  appear  similar  between  the  two,  I  comment 


only  on  the  relative  magnitudes  of  the  sediment  concentrations,  load,  and  transport 
values. 


The  C(6u,)  values  are  in  general  twice  as  high  as  those  for  the  smaller  wave, 
reflecting  the  increased  turbulent  energy  in  the  wave  boundary  layer.  As  in  the 
storm  wave  case,  the  effects  of  stratiflcation  dampen  or  wipe  out  most  of  the 
variation  in  load  and  transport  predictions  due  to  changes  in  C(6„). 

In  all  cases,  the  higher  wave  bottom  velocity  generates  larger  load  and  transport 
predictions  for  the  neutral  case  than  for  the  moderate  storm  wave  with  a  similar 
current  velocity.  This  follows  from  the  increased  shear  stress  associated  with  the 
wave.  In  the  stratified  case,  however,  the  increjise  in  wave  bottom  velocity  and 
excursion  amplitude  does  not  necessarily  correspond  to  an  increase  in  total  load 
when  the  reference  velocity  is  the  same.  (In  comparing  the  plots,  keep  in  mind 
that  the  velocity  scales  are  not  the  same.  In  Figures  5. 2-5. 5,  the  top  of  the  scale 
represents  Ur  =  64  cm! sec.  In  Figures  5.16-5.19,  the  top  of  the  scale  represents 
Ur  =  72  cm! sec.)  The  change  in  load  from  the  moderate  storm  case  to  the  large 
swell  case  varies  with  grain  size,  and  is  strongly  influenced  by  the  stratification 
effect  at  the  top  of  the  wave  boundary  layer.  For  fine  sand  (Figure  5.16  and 
Figure  5.2),  the  larger  wave  causes  an  increase  in  load  of  about  50%  when  currents 
are  strong.  This  is  about  the  same  increase  seen  for  very  fine  sand  (Figure  5.17  and 
Figure  5.3).  The  increase  in  load  for  coarse  silt  at  high  current  velocities  is  about 
a  factor  of  2  (Figure  5.18  and  Figure  5.4).  The  medium  silt  shows  an  increase  of 
about  a  factor  of  about  4  (Figure  5.19  and  Figure  5.5). 

The  magnitude  of  the  increases  gradually  declines  with  lower  current  velocities 
in  all  cases.  In  the  three  coarsest  sediments,  the  difference  becomes  very  small  at 
low  velocities;  in  some  cases  the  predicted  load  and  transport  are  the  same.  The 
differences  in  predictions  between  grain  sizes  is  also  very  small  or  nonexistent;  the 
problem  has  again  become  dominated  by  stratification  effects. 
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Figure  5.16:  Predicted  near-bottom  load  (^)  *nd  transport  {cm*/cm/«ec)  values  for 
fine  sand  in  large  swell  waves.  T=20  sec,  and  the  wave  height  in  50  m  of  water  would  be 
3.3m.  Ur  is  reference  current  velocity.  6^  ranges  from  29-31  cm;  ^  ranges  from  11  to  39 
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Figure  5.17;  Predicted  near>bottom  load  (^)  end  transport  (cm*/cm/aec)  values  for 
very  fine  sand  in  large  swell  waves.  ranges  from  22-24  cm;  ^  ranges  from  7  to  32  m. 
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Figure  5.18:  Predicted  near-bottom  load  (j^)  and  transport  (cm*/cm/*ec)  vedues  for 
coarse  silt  in  large  swell  waves.  5,  ranges  from  17*2 1  cm;  ^  ranges  from  6  to  29  m. 
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Figure  5.19:  Predicted  near-bottom  load  and  transport  (cm*/cm/aec)  values  for 
medium  silt  in  large  swell  waves.  The  irregular  contour  patterns  for  the  stratified  plots 
occur  in  regions  where  the  model  gives  non-unique  solutions,  as  discussed  in  text,  ranges 
from  14-21  cm;  ^  ranges  from  5  to  28  m. 
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Ranges  of  reworking  depths,  Large  swell 
Near-bottom  load  only 


d{cm) 

neutral  (cm) 

stratified  (cm) 

0.02 

0.027  -  1.29 

0.015  -  0.32 

0.01 

0.062  -  6.44 

0.015  -  0.95 

0.006 

0.22  -  20.45 

0.025  -  3.89 

0.0025 

0.92  -  53.0 

0.014  -  33.0 

Table  5.7:  Reworking  depths  for  large  swell,  showing  large  decrease  in  ranges  due  to 
stratification  effects.  The  maximum  reworking  depths  for  the  neutral  silt  cases  reflect  the 
unrealistically  high  predicted  concentrations  discussed  in  the  text. 

Like  the  storm  wave  case,  bedload  transport  is  insignificant  in  this  wave  en¬ 
vironment,  approaching  ten  percent  of  the  suspended  transport  volume  only  for 
the  coarsest  grains  and  the  smallest  currents.  The  larger  reference  concentrations 
yield  even  larger  stratification  effects  on  the  outer  Ekman  layer  load  and  transport 
estimates,  so  that  for  all  but  the  smallest  grain  size,  outer  Ekman  layer  transport  is 
less  than  0.8  times  the  near-bottom  transport,  and  outer  Ekman  layer  load  is  less 
than  0.4  times  the  near-bottom  load.  The  predicted  outer  Ekman  layer  load  and 
transport  values  for  medium  silt  are  a  factor  of  two  to  three  greater  than  the  near 
bottom  values,  however.  The  large  values  may  be  limited  by  factors  such  as  water 
depth  and  bed  armoring,  so  the  model’s  results  must  be  interpreted  in  light  of  local 
conditions  especially  in  applications  involving  fine  sediments  in  large  storms. 

The  reworking  depths  due  to  near-bottom  suspended  load,  including  the  effects 
of  stratification,  for  the  four  grain  sizes  in  this  wave  environment  are  very  small 
for  sands,  always  less  than  one  centimeter  (Table  5.7).  In  both  the  silt  cases,  the 
predicted  concentrations  are  within  expected  limits  in  all  but  the  strongest  current 
velocities,  which  overcome  much  of  the  effect  of  stratification. 
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5.6  Results:  General  Wave  Effects 

The  results  discussed  so  far  cover  most  of  the  patterns  observed  in  the  model  runs 
for  all  eight  wave  cases.  Plots  of  the  other  load  and  transport  predictions  are 
included  in  Appendix  C.  The  present  discussion  summarizes  those  results. 

The  effects  of  waves  on  sediment  transport  can  be  most  generally  classified 
by  bottom  wave  velocity,  uj.  The  ranges  of  predicted  loads  and  transports  (Ta¬ 
bles  5.8-5.11)  show  that,  for  each  20  cm/sec  increase  in  u^,  the  maximum  and 
minimum  ranges  for  load  and  transport  predictions  increase  by  an  order  of  magni¬ 
tude  when  grain  size  is  constant.  The  exceptions  to  this  pattern  of  increase  are  due 
to  strong  stratification  effects,  also  controlled  by  the  wve  bottom  velocities.  For 
example,  a  current  of  30  cm/sec  in  wave  conditions  with  a  10  second  wave  period 
and  20  cm/sec  wave  bottom  velocities  can  suspend  0.02-0.1  cm®/cm*  from  a  silt 
bottom,  depending  on  reference  concentration  (Figure  C.3b);  if  the  wave  bottom 
velocity  is  40  cm/sec,  the  load  is  0.04-0.06  cm®/cm*  (Figure  flgsteSb).  The  approx¬ 
imate  equality  of  the  ranges  is  due  to  stratification  effects:  the  ranges  predicted  by 
the  neutral  model  are  0.02-0.2  cm®/cm*  for  the  smaller  bottom  wave  velocity  and 
0. 2-2.0  cm®/cm*  for  the  larger  one. 

The  effect  of  period  on  sediment  transport  is  smaller  and  less  direct.  The  wave 
bottom  velocity  has  direct  control  over  shear  velocity  which  controls  the 

sediment  motion  directly  (Equation  3.22).  The  wave  period  effect  comes  through 
the  relation  of  the  wave  period  to  the  height  of  the  wave  boundary  layer  (5^,; 
Equation  2.4),  and  the  relation  of  to  the  roughness  length  for  the  mean  flow 
[zoc]  Equation  3.30).  A  shorter  wave  period  generates  a  thinner  wave  boundary 
layer,  which  means  less  roughness  is  felt  by  the  flow  and  the  sediment  load  is 
expected  to  be  slightly  smaller.  The  response  to  change  in  period  is  barely  visible 
in  predicted  load  for  silts  for  the  neutral  case  (Figure  C.7a  versus  Figure  C.lla,  as 
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Ub  Ab  T(sec)  d(cm) 

20.00  32.00  10.053  0.020 

20.00  64.00  20.106  0.020 

40.00  64.00  10.053  0.020 

40.00  96.00  15.080  0.020 

40.00  128.00  20.106  0.020 

60.00  144.00  15.080  0.020 

60.00  191.00  20.001  0.020 


load  transport 

(max)  (min)  (max)  (min) 

.7000E>02  .0000E400  .2550E400  .0000E-^00 
.5000E-^2  .0000E400  .1840E400  .0000E400 

.7700E-01  .3000E-02  .4664E401  .4000E-02 
.1050E400  .4000E-02  .6329E401  .6000E-02 
.132OE400  .  4000E-02  .7992E401  .7000E<-02 

.1740E+00  .7000E-02  .1207E+02  .1100E-01 
.1900E-f00  .9000E-02  .1233E402  .1500E-01 


100.00  319.00  20.043  0.020  .1048E401  .2300E-01  .1321E403  .5600E-O1 


Table  5.8:  Ranges  in  predicted  values  of  stratified  near  bottom  load  and  transport  for 
fine  sand  (d  =  0.02  cm).  Predictions  are  grouped  by  wave  bottom  velocity.  Load  units  are 
transport  units  are  emblem  face 


20.00 

20.00 


Ab  T(sse)  d(cffl) 

32.00  10.053  0.010 
64.00  20.106  0.010 


load  transport 

(max)  (min)  (max)  (min) 

.7100E-01  .O000E40O  .4515E401  .0000E400 
.4500E-01  .O000E400  .2462E401  .000OE40O 
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40.00  64.00  10.053  0.010 
40.00  96.00  15.080  0.010 
40.00  128.00  20.106  0.010 

60.00  144.00  15.080  0.010 
60.00  191.00  20.001  0.010 


.1980E40e  .2O00E-02  .1552E402  .709OE-02 
.2130E40e  .3000E*-«2  .1599E402  .1000E-^1 
.2340E400  .4O00E-02  .1723E402  .130OE-01 

.6280E40O  .7O0OE-«2  .7076E4«2  .2000E-01 
.5710E40O  .9000E-02  .5865E402  .2300E-01 


|J* 

V  % 


100.00  319.00  20.043  0.010  .4718E+01  .3900E-01  .1162E404  .1730E-^0O 

Table  5.9:  Ranges  id  predicted  values  of  stratified  near  bottom  load  and  transport  for 
very  fine  sand  (d  =  0.01  cm).  Predictions  are  grouped  by  wave  bottom  velocity.  Load 
units  are  transport  units  are  cm^/cm/$ec 


Ub 

Ab 

T(8ec) 

d(cm) 

load  transport 

(max)  (min)  (max)  (min) 

20.00 

32.00 

10.053 

0.006 

.1330E+00 

.0000E+00  .7377E+01 

.1000E-02 

20.00 

64.00 

20.106 

0.006 

.1380E400 

.0000E+00  .7216E+01 

.0000E400 

40.00 

64.00 

10.053 

0.006 

.6750£-f00 

.4000E-02  .7065E+02 

.1900E-01 

40.00 

96.00 

15.080 

0.006 

.5800E40O 

.5O00E-02  .5495E+02 

.2200E-01 

40.00 

128.00 

20.106 

0.006 

.5530E400 

.6000E-02  .4966E+02 

.2500E-01 

60.00 

144.00 

15.080 

0.006 

.2814E401 

.1400E-01  .5079E+03 

.8300E-01 

60.00 

191.00 

20.001 

0.006 

.2332E401 

.1500E-01  .3798E+03 

.7600E-01 

100.00 

319.00 

20.043 

0.006 

.2134E+02 

.9300E-01  .8705E+04 

.1297E+01 

Table  5.10;  Ranges  in  predicted  values  of  stratified  near  bottom  load  and  transport  for 
coarse  silt  (d  =  0.006  cm).  Predictions  are  grouped  by  wave  bottom  velocity.  Load  units 
we  transport  units  are  cm^Jcm/aec 


Ub  Ab  T(a«c)  d(cm) 

20.00  32.00  10.053  0.002 
20.00  64.00  20.106  0.002 


load 

(max)  (min) 
.3870E400  .2000E-02 
.3440E+00  .2000E-02 


transport 
(max)  (min) 
.2078E+02  .1200E-01 
.1622E+02  .8000E-02 


40.00  64.00 
40.00  96.00 
40.00  128.00 


10.053  0.002 
15.080  0.002 
20.106  0.002 


.5961E401 

.4582E401 

.3881E+01 


.2300E-01 

.1900E-01 

.1800E-01 


.9591 E+03 
.6652E403 
.5217E+03 


.2990E400 

.2160E+00 

.1830E400 


60.00  144.00  15.080  0.002 
60.00  191.00  20.001  0.002 


.2345E402 

.1978E+02 


.1020E400 

.8800E-01 


.6220E-f04 

.4945E+04 


.2464E+01 

.1845E+01 


100.00  319.00  20.043  0.002 


.1274E+03  .7550E+00 


.8944E+04  .4288E+02 


Table  5.11;  Ranges  in  predicted  values  of  stratified  near  bottom  load  and  transport  for 
medium  silt  (d  =  0.0025  cm).  Predictions  are  grouped  by  wave  bottom  velocity.  Load 
units  are  transport  units  are  cm’/cm/see 
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marked).  The  predicted  load  rises  somewhat  for  the  stratified  case  with  increasing 
period. 


5.7  Summary 

The  sensitivity  analysis  has  brought  to  light  a  number  of  considerations  which 
should  be  kept  in  mind  when  applying  the  Grant-Madsen-Glenn  model  to  questions 
of  sediment  transport  and  when  interpreting  the  model’s  results.  With  these  results 
in  hand,  the  boundary  layer  model  can  be  a  useful  geological  tool;  if  the  model  is 
applied  blindly  it  may  be  misleading  and,  in  some  cases,  simply  will  not  work. 

The  conclusions  can  be  cast  in  terms  of  physical  conditions  at  the  seafloor  under 
which  caution  should  be  used  in  applying  the  model: 

•  Large  wave  bottom  velocities  (>~  40^),  especially  when  near-bottom  cur¬ 
rents  are  small  (<'^  stratification  effects  which  reduce  load 

predictions  by  50%  or  more.  These  predictions  require  field  verification. 

•  Large  wave  bottom  velocities  occurring  where  the  seafloor  surface  grains  are 
silt  sized  or  finer  may  lead  to  nonunique  solutions  to  the  stratified  bound¬ 
ary  layer  model,  depending  on  specified  reference  velocity.  The  predicted 
sediment  load  and  transport  may  vary  by  an  order  of  magnitude  for  the 
different  solutions.  If  the  model’s  representation  of  stratification  is  correct, 
this  nonuniqueness  requires  specification  of  another  boundary  condition,  ei¬ 
ther  a  current  velocity  higher  in  the  water  column  or  a  suspended  sediment 
concentration  mezisurement. 

•  Large  bottom  currents  {>~  30^  at  z  =  Im)  can  suspend  sediments  to 
heights  in  the  water  column  greater  than  the  near-bottom  layer  on  which 
this  model  focuses.  When  estimates  of  load  and  transport  volumes  in  the 
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outer  Ekman  layer  are  greater  than  or  equal  to  the  near  bottom  volumes,  the 
model’s  predictions  of  total  transport  are  less  certain. 

•  In  a  related  concern,  the  predicted  bottom  Ekman  layer  depth  based  on  mean 
shear  velocity  can  exceed  the  water  column  depth  in  strong  current  flows.  If 
so,  the  outer  Ekman  layer  predictions  should  be  discarded. 

•  If  wave  shear  stress  is  insufficient  to  move  sediment  without  current  input, 
then  there  exists  a  set  of  input  conditions  on  which  the  model  will  not  con¬ 
verge  when  grain  roughness  is  used  as  the  roughness  length.  In  this  case,  a 
more  realistic  representation  of  roughness,  bzised  on  existing  bedforms,  should 
be  used  instead  of  the  moveable  bed  roughness. 

•  When  grain  sizes  are  large,  bedload  may  be  a  significant  factor  in  total  trans¬ 
port.  This  model  provides  bedload  estimates,  but  they  should  be  used  only 
for  estimating  the  significance  of  bedload  transport  on  a  given  shelf  under 
given  circumstances.  If  bedload  is  determined  to  be  the  primary  transport 
mechanism,  a  transport  model  which  focuses  on  bedload  should  be  used. 

The  controlling  influence  of  stratification  on  sediment  load  and  transport  that 
is  predicted  for  large  wave  bottom  velocities  has  important  implications  for  sed- 
imentology.  In  limiting  upward  transport  of  sediment  from  the  seafloor,  it  could 
explain  why  storms  do  not  move  more  sediment  than  they  do.  The  stratification 
model  requires  testing  and  verification,  however,  before  extensive  application  to  ge¬ 
ological  questions  is  warranted.  The  next  chapter  provides  some  examples  of  ways 
the  model  can  be  used  to  explore  patterns  of  sediment  transport  and  deposition. 


deposition;  if  more  sediment  is  leaving,  there  is  net  erosion.  The  deposition  rate 
for  the  box  is  determined  by  averaging  the  volume  of  net  transport  over  the  area 
of  the  grid  square.  The  average  reworking  depth  is  the  average  of  the  values  at  the 
four  corners,  including  both  the  near-bottom  and  outer  Ekman  layer  predictions. 
This  quantity  represents  the  total  volume  of  sediment  moved  by  the  waves  and 
currents.  Adding  the  reworking  depth  to  the  net  transport  for  the  duration  of  the 
storm  gives  the  total  depth  to  which  sediment  is  affected  by  the  storm. 

In  order  to  gauge  the  effects  of  specific  parameters  on  deposition,  the  physical 
scenario  is  simplified.  The  continental  shelf  is  considered  in  cross  section,  con¬ 
ceptualized  as  a  series  of  boxes  at  intervals  of  ten  meter  changes  in  water  depth 
(Figure  6.1).  The  x-direction  is  taken  to  be  cross-shelf,  and  the  unit  length  (Ai)  is 
determined  by  the  slope  of  the  shelf.  The  y-direction  is  alongshelf,  and  the  length 
of  the  box  side  in  the  y-direction  (Ay)  is  taken  to  be  the  same  as  that  in  the  x- 
direction.  Note  that  adjacent  boxes  are  autonomous;  there  is  no  assumption  that 
what  happens  in  one  square  in  any  way  affects  an  adjacent  square. 

The  directions  of  the  wave  and  current  (fl*  and  0cur)  are  given  with  respect  to 
the  alongshelf  direction,  which  is  taken  to  be  zero  degrees.  The  angle  between  the 
wave  and  the  current  {<pc)  is  the  smallest  angle  between  the  direction  of  the  wave 
and  the  direction  of  the  mean  bottom  flow,  so  <pe  is  at  most  ninety  degrees.  The 
effects  of  wave  refraction  on  9^  are  not  taken  into  account,  so  that  the  effects  of 
wave  shoaling  on  bottom  velocity  and  excursion  will  be  isolated;  the  influence  of 
wave  direction  is  small.  The  value  of  70  is  taken  to  be  0.002  in  all  cases. 

The  boxes  in  Figure  6.1  begin  at  a  depth  of  thirty  meters.  This  is  because,  for 
the  storm  cases  considered  here,  shallower  depths  tend  to  have  full  Ekman  depths 
greater  than  thirty  meters;  certainly  the  surface  and  bottom  Ekman  layers  overlap, 
and  this  is  a  complication  not  addressed  in  this  simple  model.  The  offshore  depth 
limit  is  constrained  only  by  the  depth  at  which  sediment  is  no  longer  mobilized, 


which  is  in  lurn  controlled  by  the  depth  of  influence  of  the  wave  (Section  5.1). 

Depth  contours  are  assumed  to  run  alongshelf,  so  there  is  no  variation  in  depth 
in  the  y-direction.  This  translates  to  no  variation  in  the  bottom  wave  parameters 
with  y.  The  ten  meter  variation  in  water  depth  introduces  a  significant 
change  in  bottom  wave  parameters  in  the  x-direction,  however,  and  the  net  trans¬ 
port  introduced  by  this  variation  will  be  the  focus  of  this  discussion.  Variation  in 
wave  parameters  introduces  no  net  transport  in  and  of  itself,  of  course,  since  the 
GMG  model  is  based  on  a  linear  wave  assumption;  transport  is  carried  out  by  the 
mean  current,  so  if  there  is  no  current,  there  is  no  net  transport.  If  the  only  varia¬ 
tion  in  parameters  is  in  the  x-direction  (e.g.  depth  change,  variation  in  grain  sizes 
across  shelf)  and  all  flow  is  alongshelf,  there  is  transport  but  no  net  transport,  since 
the  same  volume  of  sediment  which  flows  into  the  square  also  flows  out  at  the  other 
end.  If  hydrodynamic  and  sedimentological  parameters  vary  only  perpendicular  to 
the  flow,  then  their  variation  will  not  contribute  to  the  net  transport  of  sediment. 

In  the  presence  of  a  cross-shelf  component  to  the  near-bottom  meam  flow,  the 
variation  in  bottom  wave  parameters  with  depth  can  effect  net  transport  by  sev¬ 
eral  means.  The  most  obvious  is  the  increasing  bottom  velocity  of  waves,  due  to 
shoaling.  The  maximum  shear  stress  is  expected  to  decline,  so  the  sediment  made 
available  for  transport  by  the  current  would  also  decline.  The  expected  transport 
pattern  in  the  presence  of  an  offshore  component  of  current  would  be  a  progres¬ 
sively  declining  offshore  transport  resulting  in  a  net  deposition  of  sediment  at  all 
locations  where  sediment  was  put  in  motion.  To  illustrate  this  point,  consider  Fig¬ 
ure  6.1.  For  the  same  crosshelf  current  component  at  all  depths,  assumed  directed 
offshore,  the  average  offshore  transport  into  a  box  at  a  depth  of  thirty  meters  would 
be  greater  than  the  transport  out  of  the  box  at  forty  meters,  since  the  wave  bot¬ 
tom  shear  stress  is  smaller  at  the  deeper  location.  Since  more  enters  than  leaves, 
this  situation  would  result  in  deposition.  The  same  would  be  true  for  each  of  the 
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Figure  6.1:  Schematic  diagram  of  grid  squares  used  to  represent  the  continental  shelf  in 
cross  section,  beginning  at  a  depth  of  thirty  meters  and  extending  across  the  shelf  at  depth 
intervals  of  ten  meters.  Scale  in  the  x-direction  (cross-shelf)  determined  by  slope  of  shelf; 
scale  in  y-direction  determined  by  alongshelf  changes  in  input  parameters,  or  by  x-scale. 
Wave  and  current  directions  specified  with  respect  to  alongshelf  direction. 


other,  deeper,  boxes.  In  this  situation,  the  only  question  is  the  amount  of  expected 
deposition.  In  the  discussion  which  follows,  this  pattern  of  deposition  is  referred 
to  as  the  ‘expected’  pattern. 

As  the  sensitivity  analysis  showed,  waves  affect  two  other  parameters  impor¬ 
tant  to  sediment  transport:  bed  roughness  and  sediment  stratification.  Wave  shear 
stresses  can  create  ripples;  stronger  stresses  can  wash  them  out.  The  lower  rough¬ 
ness  under  larger  wave  velocities  can  lead  to  situations  where  the  enhancement  of 
sediment  transport  by  waves  is  smaller  at  shallower  depths  than  deeper.  Likewise, 
stratification  effects  can  be  stronger  under  larger  waves  (see  Section  5.3),  which 
would  also  lead  to  a  decrease  in  suspended  sediment  transport  at  smaller  depths. 
These  factors  tend  to  counter  the  ‘expected’  pattern  of  wave  effects  described  in 
the  last  paragraph.  The  significance  of  each  one  for  deposition  and  erosion  patterns 
is  not  obvious,  however.  The  applications  below  explore  some  possible  variations. 

To  exercise  the  model  in  a  quasi-realistic  way,  net  transport  is  estimated  for 
areas  of  the  continental  shelf  modelled  on  the  mid- Atlantic  region  of  the  east  coast 
of  the  U.S.  and  a  region  of  the  northern  California  coast  where  wave,  current,  and 
sediment  conditions  are  reasonably  well  understood.  These  two  shelves  we  used 
because  they  constitute  a  pair  of  very  different  shelf  environments  themselves,  but 
contain  features  which  are  common  to  many  shelves  worldwide;  each  has  been 
closely  studied  in  recent  years,  with  data  including  extensive  near-bottom  mea¬ 
surements  of  wave,  current  and  sediment  conditions.  That  data  is  the  source  for 
the  storm  conditions  used  in  the  model  runs. 

6.1  Mid- Atlantic  Bight  shelf  type 

The  Mid- Atlantic  Bight,  encompassing  the  area  north  of  Cape  Hatteras  smd  south 
of  Cape  Cod  (Figure  6.2),  is  typical  of  continental  shelves  on  passive  continental 
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Figure  6.2:  The  AtUntic  Continental  shelf,  showing  the  distribution  of  grain  sizes  and 
types.  The  Mid- Atlantic  Bight  is  marked,  and  the  location  of  the  silty  ‘Mud  Patch*  area  is 
shown.  The  200m  coutour  shows  the  approximate  location  of  the  shelf  edge  (from  Milliman 
et  al.,  1972). 

margins.  It  is  wide  in  most  places,  widening  from  about  23  km  at  Cape  Hatteras 
to  150  -170  km  off  New  Jersey  and  Cape  Cod.  Depth  contours  are  generally  shore- 
parallel.  The  shelf  break  is  close  to  the  world  average  of  130  m,  and  the  average 
bottom  slope  ranges  from  about  5  x  10"’  to  3  x  10"^.  Outside  of  a  few  estuarine 
areas,  the  entire  Atlantic  margin,  from  Camada  to  Florida,  is  covered  with  sand¬ 
sized  materiad,  mostly  relict  from  lower  stands  of  sea  level.  The  single  exception  is 
a  region  south  of  Cape  Cod  in  water  depths  of  sixty  to  one  hundred  meters  where 
silt-  and  clay-sized  grains  make  up  25  to  75  per  cent  of  the  surface  bed  material 
(Shepard  and  Cohee,  1936;  Bothner,  et  al,  1981).  It  is  referred  to  locally  as  the 
‘Mud  Patch*. 
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Under  everyday,  non-storm  conditions  on  the  mid-to-outer  Atlantic  continental 
shelf,  sediment  does  not  move.  Away  from  nearshore  regions,  tidal  currents  and 
mean  flows  are  typically  less  than  10  cm/sec  (Butmsin,  1987a).  Waves  are  locally 
wind-generated;  exclusive  of  storms  they  are  too  small  to  generate  bottom  velocities 
sufficient  to  move  the  coarse  sediment  typical  of  this  margin.  Using  bottom  pressure 
observations  collected  between  1975  and  1980  at  depths  of  60  and  80  meters  on 
Georges  Bank  and  at  80  meters  in  the  Mud  Patch,  Butman  (1987a)  determines 
that  maximum  wave  bottom  velocities  are  almost  never  greater  than  10  cm/sec 
during  the  summer  months,  and  reach  that  velocity  only  one  per  cent  of  the  time 
in  the  fall  and  five  per  cent  of  the  time  in  spring.  The  winter  is  characterized  by 
about  two  storms  per  month,  however,  and  Butman  suggests  that  wave-generated 
bottom  velocities  are  greater  than  10  cm/sec  up  to  15  per  cent  of  the  time  in  winter. 
Not  even  the  silts  of  the  Mud  Patch,  and  certainly  not  the  medium  to  co2urse  sands 
which  make  up  the  bulk  of  the  surface  sediment  on  the  Atlantic  continental  margin 
respond  to  wave  and  current  velocities  less  them  10  cm/sec.  We  therefore  focus  on 
representative  storm  conditions. 

In  studies  from  1976  through  the  present,  Butman  has  been  measuring  near- 
bottom  conditions  on  the  mid-Atlantic  shelf,  often  in  conjunction  with  physical 
oceanographic  studies  of  shelf  circulation  (see,  e.g.  Butm2m,  1987a,  1987b;  Butman, 
et  al,  1979).  Butman’s  data  suggest  that  sediment  resuspension  is  highly  correlated 
with  storm  activity  (Butman,  1987a).  Sediment  as  coarse  as  fine  sand  has  been 
found  in  sediment  traps  as  high  as  three  meters  above  the  bottom  (Parmenter, 
Bothner  and  Butman,  1983),  and  suspended  sediment  concentrations  of  at  least 
10  and  perhaps  50  mg/1  (equivalent  to  volumetric  concentrations  of  4  x  10“®  to 
2  X  10“®)  have  been  observed  at  heights  of  two  meters  above  the  bottom  (Butman, 
1987a). 

Storms  affecting  the  mid-Atlantic  shelf  generally  follow  one  of  two  patterns: 


1)  winds  blow  from  the  south,  shifting  to  southwest  over  the  course  of  the  storm, 
setting  up  bottom  currents  which  are  generally  alongshelf  toward  the  north  with  a 
slight  offshore  component;  or  2)  winds  blow  from  the  northeast,  shifting  to  north¬ 
west,  setting  up  bottom  currents  which  flow  generally  alongshelf  toward  the  south, 
with  a  small  offshore  component  (Butman,  1987b).  Wave  direction  is  likewise 
controlled  by  local  winds.  Storm  effects  last  up  to  about  three  days. 

These  observations  of  shelf  and  storm  conditions  provide  guides  to  the  input 
conditions  to  this  model  and  interpretations  of  results  for  the  first  application. 
It  is  to  be  emphasized  that  these  are  not  specifically  results  for  any  location  on 
the  Atlantic  continental  margin,  but  are  meant  to  illustrate  the  application  of  the 
techniques  developed  here  as  they  can  be  applied  to  that  margin  or  one  like  it. 

Application  One:  Large  storm,  single  grain  size  bed 

The  first  application  of  the  boundary  layer  model  to  a  continental  shelf  uses  wave 
and  current  conditions  observed  during  a  large  storm  on  the  Atlantic  continental 
shelf.  The  seafloor  is  assumed  to  be  covered  with  sediment  of  a  single  grain  size. 
The  range  of  water  depths  covered  is  thirty  to  one  hundred  meters,  in  a  series  of 
seven  grid  squares,  each  eighteen  kilometers  on  a  side  and  each  spanning  a  ten 
meter  change  in  water  depth  (Figure  6.3).  This  box  configuration  corresponds  to 
a  slope  of  0.055  %.  The  wave  and  current  conditions  used  are  modeled  after  data 
collected  in  February,  1978,  by  Butman  (1987a).  This  period  represents  an  extreme 
winter  storm  for  this  region;  the  storm  is  locally  referred  to  (still)  as  the  ‘Blizzard 
of  ’78’  (always  capitalized).  The  mean  current  was  approximately  32  cm/sec  for 
the  duration  of  the  storm,  and  had  a  cross-shelf  component  near  the  bed  of  about 
5  cm/sec.  This  translates  to  a  current  direction  of  161°,  where  0°  is  the  alongshelf 
direction.  Since  waves  are  wind-generated  and  the  wind  in  this  storm  blew  initially 
out  of  the  northeast,  a  direction  of  225°  is  used  for  the  wave.  The  wave  conditions 
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Figure  6.3:  Schematic  diagram  of  grid  representing  the  Mid-Atlantic  Bight  continental 
shelf,  emd  wave  and  current  conditions  used  in  applications  described  in  text 

over  approximately  three  days  sustained  heights  of  five  meters  with  twelve  second 
periods.  The  model  was  run  for  each  of  four  grain  sizes  ranging  from  medium  silt 
to  medium  sand. 

The  results  are  shown  in  Figure  6.4.  A  lot  of  information  is  packed  into  each 
graph,  so  the  format  is  explained  in  some  detail  here.  Each  part  (a,b,c,d)  of  the 
figure  shows  the  predictions  for  a  single  grain  size  for  a  row  of  seven  boxes  across 
the  shelf,  from  thirty  to  one  hundred  meters  water  depth.  The  total  deposition  rate 
(or  erosion  rate  if  the  value  is  negative),  in  cm/hr,  is  shown  for  each  grid  square  in 
the  left-hand  bar  graph,  and  repeated  in  the  right-hand  list.  The  other  three  bar 
graphs  show  the  proportion  of  the  total  deposition  rate  attributable  to  each  of  the 
three  contributions  (bedload,  near-bottom,  and  outer  Ekman  layer  transports). 
A  transport  rate  of  4  X  10~*  cm/hr  tramslates  to  a  net  accumulation  of  about 
0.1  mm/day.  The  average  reworking  depth  is  listed,  suggesting  gross  transport 
through  the  square,  as  opposed  to  the  net  transport  represented  by  deposition 
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The  variation  of  wave  parameters  with  water  depth  can  introduce  net  cross¬ 
shelf  transport  in  the  presence  of  a  steady  current  with  even  a  small  cross-shelf 
component.  However,  the  depositional  patterns  vary  dramatically  from  the  intu¬ 
itively  ‘expected’  pattern  of  net  deposition  in  all  boxes.  Particularly  in  the  sand 
ca.ses,  the  decrease  in  bottom  wave  velocity  with  depth  leads  to  increase  in  rough¬ 
ness,  because  waves  strong  enough  to  wzish  ripples  out  at  one  depth  are  only  strong 
enough  to  create  ripples  in  water  ten  meters  deeper.  Predicted  transport  out  of 
the  box  at  the  deeper  end  is  thus  increased  relative  to  input  at  the  shallower  end. 
This  can  be  augmented  by  a  larger  sediment  stratification  effect  at  the  shallower 
end  of  the  box,  which  limits  input.  The  stratification  effect  dominates  in  silt  cases. 
Recall  that  the  deposition  rates  represent  only  cross-shelf  transport;  exactly  as 
much  sediment  enters  the  square  from  the  north  as  leaves  toward  the  south,  since 
there  is,  by  assumption,  no  alongshelf  variation  in  any  parameter. 

To  illustrate  the  effects  of  stratification  and  roughness  on  the  net  transport 
predictions,  two  sets  of  results  will  be  examined  in  more  detail  and  related  to  the 
sensitivity  results.  First,  we  examine  the  strip  from  fifty  to  seventy  meters  depth 
in  the  medium  sand  case  (Figure  6.4a).  In  the  box  covering  fifty  to  sixty  meters 
depth,  a  net  erosion  of  1  x  10“^  cm/hr  (0.02  mm/day)  is  predicted;  deposition  of 
1.3  X 10“^  cm/hr  is  predicted  for  the  sixty  to  seventy  meter  box.  This  change  results 
primarily  from  moveable  bed  roughness  effects  due  to  the  change  in  mziximum  wave 
bottom  velocity  from  59  cm/sec  to  36  cm/sec  from  50  to  70  meters,  and  in  wave 
excursion  amplitude  Ai,  from  113  to  68  cm  over  that  distance.  The  skin  friction 
shear  stress  is  approximately  equal  to  the  breakoff  value  (Equation  3.55)  at  fifty 
meters,  but  lower  than  the  breakoff  value  at  sixty  meters;  wave-generated  ripples 
are  therefore  wzished  out  at  the  shallower  depth,  but  not  at  the  deeper  location 
(see  Section  3.3).  The  value  of  the  boundary  roughness  zq  is  0.8  cm  at  fifty  meters 
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Figure  6.4:  Predicted  depoaition/erosion  rates  aod  reworking  depths  for  single-grain-size 
shelves,  with  shelf  and  storm  description  derived  from  conditions  typical  of  continental 
shelf  of  the  Mid- Atlantic  Bight.  Grid  is  shown  in  Figure  6.3. 

H  =  5m,  T  =  12s,  =  225*,  u,  =  32  em/»ee,$cuf  =  161* 

(a)  Results  for  medium  sand  show  variable  deposition  ratss  across  the  shelf,  ranging  no 
higher  than  O.OS  mm/day.  Variation  is  due  in  large  measure  to  variation  in  bed  roughness. 
Very  little  variation  in  reworking  depths  is  evident.  Motion  ceases  at  100  m. 


b)  Deposition, fine  sand, storm,  neutral 
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(b)  Results  for  fine  sand  show  deposition/erosion  rates  and  reworking  depths  no  higher 
than  0.3  mm/day,  with  variation  due  to  roughness  and  stratification.  Motion  ceases  at 
100  m. 


Figure  6.4:  (c)  Results  for  coarse  silt  show  net  transport  rates  generally  increasing  with 
depth,  with  erosion  as  high  as  1.5  mm/day  on  the  central  shelf  and  deposition  as  high  as 
0.5  mm/day  only  10'20  m  deeper,  as  effects  of  stratification  decline  with  wave  strength 
in  deeper  water.  Motion  ceases  at  100  m,  leading  to  the  large  (and  questionable)  outer 
Ekman  layer  deposition  in  the  deepest  box. 
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(d)  Results  for  medium  silt  show  similar  effects  of  depth  and  stratification  as  the  coarse  silt 
case,  with  erosion  as  high  as  2.5  mm/day.  The  largest  predicted  erosion  results  from  an 
ambiguous  stratification  effect,  discussed  in  text.  Motion  ceases  at  110  m.  Stratification 
effects  limit  reworking  and  net  erosion  to  relatively  small  volumes  on  inner  shelf. 
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and  1.8  cm  at  sixty  meters,  representing  a  change  in  ripple  steepness  (^)  from  0.08 
to  0.13  as  predicted  using  the  Grant  and  Madsen  (1982)  moveable  bed  roughness 
model. 

Because  of  the  increased  roughness,  the  wave-current  shear  velocity  Utew  is 
greater  at  sixty  meters  than  at  fifty  meters,  even  though  the  maximum  wave  bottom 
velocity  ut,  drops  from  59  to  46  cm/sec.  The  larger  u,e«  increases  the  predicted 
load  and  transport  in  two  ways.  First,  more  sediment  can  be  transported  to  the 
top  of  the  wave  boundary  layer,  so  the  current  has  more  sediment  available  to 
suspend.  Second,  more  low-momentum  fluid  is  transported  upward,  so  the  wave- 
enhanced  roughness  felt  by  the  meaji  flow  (zoci  Section  3.3)  is  greater.  In  this 
case,  zoc  increases  from  8  cm  to  10  cm  between  fifty  and  sixty  meters  depth.  The 
increase  in  roughness  at  sixty  meters  depth  thus  generates  the  conditions  which 
lead  to  net  erosion  between  fifty  and  sixty  meters.  By  contrast,  roughness  lengths 
at  sixty  and  seventy  meters  are  approximately  equal,  so  the  expected  pattern  holds: 
the  decreaise  in  maximum  wave  bottom  velocity  leads  to  decreases  in  wave  shear 
velocity,  current  roughness  length,  sediment  concentration  at  the  top  of  the  wave 
boundary  layer,  and  total  load  2md  transport.  The  transport  into  the  box  at  sixty 
meters  depth  is  therefore  greater  than  the  transport  out  at  seventy  meters,  so  there 
is  net  deposition. 

The  net  erosion  predicted  between  sixty  and  seventy  meters  on  a  shelf  covered 
with  medium  silt  (Figure  6.4d)  is  due  to  stratification  effects,  rather  than  roughness. 
The  physical  roughness  lengths  [zq)  at  the  two  depths  are  nearly  equal:  0.028  cm 
at  sixty  meters  vs.  0.022  cm  at  seventy  meters,  reflecting  ripple  heights  of  less  than 
a  centimeter.  The  larger  wave  bottom  velocities  at  the  shallower  location  therefore 
lead  to  a  concentration  at  the  top  of  the  wave  boundeuy  layer  which  is  twice  as 
large  at  sixty  meters  as  at  seventy.  Recall  from  the  sensitivity  analysis  discussion 
of  stratification  effects  (Section  5.3.2)  that  an  increase  in  C{6cv)  of  a  factor  of  two 


could  generate  a  substantial  decrezise  in  predicted  trzmsport  by  a  given  current 
due  to  increzised  stratification  effects  (Figure  5.5d).  Stratification  therefore  limits 
transport  into  the  box  at  sixty  meters  more  than  it  limits  transport  out  at  seventy 
meters,  and  erosion  results. 

Under  a  large  wave,  large  current  and  small  grain  size  conditions  of  the  sen¬ 
sitivity  analysis  results  suggested  caution,  since  results  could  be  ambiguous.  The 
sudden  increase  in  reworking  depth  from  0.24  cm  to  1.3  cm  between  the  60-70  meter 
box  and  the  70-80  meter  box  suggests  that  this  case  should  be  investigated  (Fig¬ 
ure  6.4d).  Examination  of  the  model  results  for  that  box  reveals  that  this  wave 
and  current  case  at  seventy  meters  has  a  non-unique  solution  (Section  5.3.2).  The 
model  converged  on  one  solution  at  seventy  meters  for  the  shallower  box,  and  a 
different  solution  for  the  deeper  box.  At  70  m  in  the  60-70  meter  box,  the  predicted 
current  shear  velocity  u«e  is  0.97  cm/sec,  the  predicted  reworking  depth  is  0.35  cm 
and  total  predicted  transport  is  9.3  cm®/cm/sec.  At  70  m  in  the  70-80  meter  box, 
stratification  is  less  effective  in  reducing  turbulence,  and  the  predicted  parameters 
are  1.3  cm/sec  for  current  shear  velocity,  0.85  cm  predicted  reworking  depth,  and 
31.6  cm*/cm/sec  predicted  transport.  The  net  trzmsport  rates  are  not  affected  by 
the  inconsistency:  either  result  gives  erosion  in  each  box.  However,  this  exeimple 
demonstrates  again  that  caution  must  be  used  in  applying  this  model  to  fine  grain 
sizes  in  storm  situations  until  field  applications  have  resolved  some  of  the  uncer¬ 
tainties  and  refined  the  model,  most  importantly  the  treatment  of  stratification. 

Stepping  bzw:k  from  these  detailed  examinations  of  the  mechajiics  of  the  theoret¬ 
ical  model,  some  useful  observations  can  be  made  about  net  transport  in  storms. 
The  deposition  and  erosion  rates  in  Figure  6.4  would  introduce  net  erosion  and 
deposition  of  a  few  tenths  of  a  millimeter  to  a  few  millimeters  over  the  course  of 
a  single,  three-day  storm.  The  exceptions  are  the  deepest  boxes  in  the  silt  cases, 
where  a  two-day  storm  might  lead  to  one  to  two  cm  of  erosion  or  deposition.  In 


the  coarse  silt  case,  this  results  because  the  deepest  box  has  no  sediment  motion 
predicted  at  its  offshore  edge.  This  is  the  reason  for  large  deposition  rates  at  the 
offshore  edge  of  a  number  of  the  profiles,  particularly  for  fine  sediments. 

Comparison  of  deposition  rates  at  similar  depths  shows  a  tendency  for  net  trams- 
port  rates  to  rise  with  decreasing  grain  size,  especially  in  the  deeper  boxes  where 
stratification  effects  are  less  limiting  to  silt  resuspension.  Stronger  stratification  ef¬ 
fects  where  wave  velocities  are  greater  (i.e.  shallower  depths)  explain  the  consistent 
erosion  of  silts  in  the  shallower  areas.  Stratification  effects  keep  the  magnitudes 
of  the  total  deposition  rates  fairly  similar  between  grain  sizes,  though  there  is  a 
slight  tendency  for  net  transport  rates  to  rise  with  decreasing  grain  size.  Rework¬ 
ing  depth  increases  fairly  consistently  with  decreasing  grain  size.  The  exceptions 
to  this  pattern  follow  from  the  observations  in  the  sensitivity  analysis  that  in  some 
storm  situations,  load  (and  therefore  reworking  depth)  is  nearly  independent  of 
grain  size. 

Bedload  transport  is  insignificant  in  every  case,  and  nezir-bottom  transport  is 
more  significant  than  outer  Ekman  layer  transport  in  most  cases.  When  outer 
Ekman  layer  transport  is  significant,  however,  net  transport  rates  are  suspect.  The 
poor  constraints  on  the  outer  Ekman  transport  direction  could  cause  large  errors. 

For  contrzist,  the  predictions  for  the  neutral  case  are  shown  in  Figure  6.5.  The 
pattern  of  deposition  is  much  closer  to  the  expected  one:  in  general,  there  is  net  de¬ 
position  in  each  box,  decreasing  in  magnitude  with  depth.  This  pattern  is  changed 
only  by  roughness  effects,  most  notably  in  the  50-60  meter  box  in  the  medium 
sand  case  where  erosion  is  predicted  (Figure  6.5a).  This  predicted  erosion  results 
from  the  same  variation  in  roughness  2is  in  the  stratified  case  discussed  above,  al¬ 
though  the  erosion  rate  is  five  times  as  great  in  the  neutral  case.  The  occurrence 
of  an  erosional  region  and  the  irregularity  in  deposition  rates  across  the  shelf  for 
the  sand  czises  testify  to  the  significance  of  moveable  bed  roughness  in  combined 
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wave-current  flows. 

In  general,  the  deposition  rates  are  at  least  a  factor  of  two,  and  more  often  an 
order  of  magnitude  or  more,  greater  when  the  effects  of  stratification  are  neglected. 
The  same  holds  true  for  reworking  depths.  Stratification  effects  generally  increase 
with  smaller  depths  and  smaller  grain  sizes,  as  discussed  in  the  sensitivity  analysis, 
where  smaller  depths  correspond  to  larger  waves. 

A  word  should  be  said  about  the  plausibility  of  some  of  the  neutral  Ctise  pre¬ 
dictions  listed.  They  have  been  included  simply  as  predicted  by  the  model,  using 
it  as  a  ‘black  box’  for  input  and  output;  the  results  must  be  interpreted  in  light 
of  known  physical  conditions.  The  reworking  depths  predicted  for  silt-sized  grains 
are  impossible,  extending  over  a  meter  in  some  cases  and  over  ten  meters  in  one 
case  (Figure  6.5c  and  d).  The  accompanying  transport  rates  would  build  small 
mountains  on  the  seafloor  over  the  course  of  a  storm,  and  predicted  sediment  con¬ 
centrations  as  high  as  twenty  meters  above  the  seafloor  are  on  the  order  of  10~*, 
which  is  impossible.  The  near-bottom  layer  thickness  is  approximately  equal  to 
the  water  depth  in  the  shallowest  water  depths,  so  the  outer  Ekman  layer  trans¬ 
port  should  be  disregarded;  this  resolves  some  of  the  problems  with  the  extreme 
deposition  rates  in  the  30-50  m  range,  since  the  outer  Ekman  layer  transport  is  50 
times  the  near-bottom  transport  at  those  depths,  but  not  with  the  concentrations. 
The  sediment  concentration  parameter,  ')o,  used  in  these  predictions  has  the  value 
0.002;  decreasing  it  even  by  an  order  of  magnitude  would  still  give  concentrations 
that  are  far  too  high.  Comparisons  of  the  neutral  case  predictions  with  observed 
concentrations  demonstrates  again  that  some  modification  to  the  steady  flow  mod¬ 
eling  of  suspended  sediment  transport  must  be  included  when  there  are  strong  wave 
effects. 

The  near-bottom  sediment  concentations,  transport  estimates,  and  reworking 
depths  predicted  by  the  stratified  model,  however,  are  plausible  with  respect  to 
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Figure  6.5:  Predicted  depoeition/eroeion  rates  and  reworking  depths  for  single-grain-size 
shelves,  same  conditions  as  Figure  6.4,  but  disregarding  the  effects  of  stratification.  Re¬ 
working  depths  and  maximum  deposition  rates  are  consistently  higher  in  the  neutral  case, 
as  expected. 

H  =  5m, T  =  12s, =  225“, u,  =  32  cmfttc^tcvr  =  161“ 

(a)  Results  for  medium  sand  show  deposition  rates  varying  from  the  expected  pattern  due 
to  bottom  roughness  variability  (see  text). 


b)  Deposition, fine  sand, storm,  neutral 
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(b)  Results  for  fine  sand  show  deposition  rates  as  high  as  2.5  cm/day  in  shallow  water,  and 
some  variation  due  to  roughness. 
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c)  Deposition,  coarse  silt,  storm,  neutral 
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Figure  6.5:  (c)  Results  for  coarse  silt  sbow  unrealisticedly  high  reworking  depths  and 
transport  rates  in  shallow  water,  dominated  by  outer  Ekman  layer  transport. 


<j)  Deposition,  med.  silt,  storm,  neutral 
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(d)  Results  for  medium  silt  show  similar  effects  of  depth  and  stratification  as  the  coarse 
silt  case,  with  deposition  as  unrealistically  high  as  a  meter  per  day  in  shallow  water.  Most 
of  this  is  due  to  the  outer  Ekman  layer  transport  estimate. 
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storm  observations.  When  strong  storm  waves  are  present  on  a  continental  shelf, 
particularly  when  sediment  is  fine-grained,  the  neutral  case  formulation  for  sedi¬ 
ment  transport  predictions  is  unsatisfactory.  Stratified  effects,  as  modeled  here, 
appear  to  provide  a  reasonable  correction  to  the  simple  neutral  caise,  though  fur¬ 
ther  work  is  necessary  to  see  if  this  stratification  effect  can  be  observed  in  the  field, 
or  whether  a  similarly  satisfactory  correction  can  be  achieved  through  revision  of 
the  eddy  viscosity,  reference  concentration,  or  other  model  features. 

Application  Two:  Small  storm,  single  grain  size  bed 

This  case  is  the  same  as  the  previous  one,  except  that  the  wave  and  current  condi¬ 
tions  represent  a  less  violent  storm.  The  wave  height  is  three  meters,  wave  period 
is  ten  seconds,  and  current  velocity  at  one  meter  off  the  bottom  is  20  cm/sec. 
Wave  and  current  directions  are  the  same:  161°  for  the  current,  225°  for  the  wave. 
The  results  are  shown  in  Figure  6.6;  for  contrast,  the  neutral  results  are  shown  in 
Figure  6.7. 

There  is  no  sediment  motion  at  depths  greater  than  sixty  or  seventy  meters, 
depending  on  grain  size.  The  reworking  depth  is  often  smaller  than  the  grain  size 
for  sand-sized  grains  (Figure  6.6a  and  b),  indication  that  not  even  the  entire  surface 
layer  is  mobilized.  There  is  mostly  slight  erosion  of  the  silts.  This  contrasts  with 
the  neutral  case  (Figure  6.7),  where  there  is  at  least  some  deposition  in  all  regions, 
and  substantial  accumulation  for  medium  silts.  Even  for  this  small  storm,  the 
reworking  depths  and  concentrations  high  in  the  water  column,  predicted  using 
the  neutral  formulation,  are  unreasonably  high,  indicating  that  stratification  or 
some  other  mechanism  must  act  to  limit  resuspension. 

The  principal  result  is  that  most  sediment  motion,  especially  on  the  mid-  to 
outer-shelf  results  only  from  large  storms.  If,  as  Butman  (1987b)  suggests,  along- 
shelf  transport  direction  alternates  between  two  primary  storm  directions,  then  one 
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Figure  6.6:  Predicted  deposition/eroaion  rates  and  reworking  depths  for  single-grain-size 
shelves,  with  shelf  and  storm  description  derived  from  conditions  typical  of  a  small  storm 
on  the  Mid-Atlantic  Bight.  Grid  is  shown  in  Figure  6.3.  Sediment  is  set  in  motion  no 
deeper  than  60-70  m,  and  net  transport  is  no  greater  than  1  mm/day,  demonstrating  the 
ineffectiveness  of  small  storms  at  transporting  sediment. 

H  =  3m,  T  =  10s,  =  225®,  Ur  =  20  cm/aec,$eur  =  161® 

(a)  Results  for  medium  swd  show  insignificant  net  transport  rates  and  reworking  depths, 
following  expected  patterns  of  deposition  variance  with  depth. 
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(b)  Results  for  fine  sand  are  very  small,  but  the  patterns  are  strongly  affected  by  roughness 
and  stratification. 
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Figure  6.6:  (c)  Results  for  coarse  silt  show  the  effects  of  decreasing  stratification  with  de¬ 
creasing  wave  velocity  in  the  reworking  depth  that  increases  with  depth,  and  the  2  mm/day 
erosion  rate  on  the  central  shelf. 
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(d)  Results  for  medium  silt  show  similar  effects  of  depth  and  stratification  as  the  coarse 
silt  case. 
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Figure  6.7:  Predicted  deposition/erosion  rates  and  reworking  depths  for  single-grain-size 
shelves,  with  shelf  and  storm  description  derived  from  conditions  typical  of  a  small  storm 
on  the  Mid- Atlantic  Bight,  neglecting  the  effects  of  stratification.  Sediment  is  set  in  motion 
no  deeper  than  60-70  m,  and  net  transport  is  no  greater  than  1  mm/day,  demonstrating 
the  ineffectiveness  of  small  storms  at  transporting  sediment. 

H  =  Zm,T  =  10s,  =  225®,  Ur  =  20  cm/sec,9cur  =  161“ 

(a)  Results  for  medium  sand  show  insignificant  net  transport  rates  and  reworking  depths, 
following  expected  patterns  of  deposition  variance  with  depth. 
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(b)  Results  for  fine  sand  are  very  small,  but  the  patterns  are  strongly  affected  by  roughness 
and  stratification. 
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(d)  Results  for  medium  silt  show  simil&r  expected  depoeitioD  patterns  as  the  coarse  silt 
case,  and  the  same  depth  limitation  on  the  outer  Ekman  layer  transport  holds  with  respect 
to  the  unreasonably  high  reworking  depths  and  deposition  rates  in  the  shallowest  boxes. 
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Figure  6.8:  Initi&l  grain  bed  concentrations  in  application  with  alongshelf  grain  texture 
changes.  The  sediment  is  silty  sand  with  17%silt  upcurrent,  35%  downcurrent. 


might  expect  the  top  few  millimeters  of  a  sandy  seafloor  on  the  Atlantic  margin  to 
be  tramsported  back  and  forth  alongshelf,  with  a  slow  accumulation  or  erosion  of 
the  sand-sized  grains  due  to  cross-shelf  transport  of  perhaps  a  millimeter  or  two 
per  year,  depending  on  water  depth. 

Application  Three:  Alongshelf  grain  texture  changes 

On  the  Atlantic  continental  shelf  as  one  approaches  the  ‘Mud  Patch’,  there  b  a 
perceptible  fining  of  grain  texture  alongshelf.  In  order  to  test  the  effect  of  alongshelf 
grain  texture  changes,  initial  bed  concentrations  of  83%  fine  sand  {d  ~  0.015  cm) 
and  17%  coarse  silt  (d  =  0.006  cm)  at  the  upcurrent  points  and  65%  fine  sand  with 
35%  coarse  silt  at  the  downcurrent  points  are  assumed  (Figure  6.8).  The  same 
‘Blizzard  of  ’78’  wave  conditions  are  used  [H  =  5m;  T  =  12sec;  Ow  =  225”),  but  the 
32  cm/sec  current  is  taken  as  heading  directly  alongshelf  (0eur  =  180°)  so  that  all 
net  transport  results  from  grain  texture  changes. 

The  results  are  shown  in  Figure  6.1,  both  total  deposition  rates  and  reworking 
depths  (Figure  6.1a3  and  by  grain  size  (Figure  6.1b).  Note  that  despite  the  fact 
that  the  bed  is  predominately  sand,  both  total  trzmsport  (as  signaled  by  reworking 
depth)  and  net  treuisport  are  predominately  silt.  Bedload  is,  again,  insignificant, 
but  outer  Ekman  layer  transport  plays  a  large  role  in  the  net  transport  in  several 
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boxes.  This  is  of  less  concern  here,  however,  since  net  transport  predictions  are  not 
dependent  on  a  small  component  of  the  total  flow  as  they  were  in  the  last  cases; 
some  veering  in  the  outer  Ekman  layer  will  still  leave  the  strongest  component  of 
transport  alongshelf. 

Intuitively,  one  expects  net  erosion,  since  the  silt  would  be  more  likely  to  be 
transported  and  there  is  more  silt  on  the  south  end  of  the  box,  which  is  the  direction 
of  transport.  This  pattern  holds  for  areas  deeper  than  sixty  meters.  Shallower 
than  that,  however,  the  larger  wave  bottom  velocities  generate  larger  suspended 
sediment  concentrations  at  the  top  of  the  wave  boundary  layer,  and  stratification 
effects  are  pronounced.  The  larger  proportion  of  silt  farther  south  is  more  affected 
by  the  stratification  effects,  so  much  so  that  the  net  tramsport  out  at  the  southern 
end  of  the  box  is  lower  than  the  net  transport  in  at  the  northern  end.  This  explains 
the  predicted  deposition  at  depths  less  than  sixty  meters. 

The  predicted  erosion  rates  are  significant,  representing  a  centimeter  or  more 
over  the  course  of  a  large  storm.  A  prevailing  pattern  of  this  sort  would  soon 
winnow  away  all  fine  sediment  if  there  were  no  source  of  silts  from  the  north. 

6.2  Northern  California  shelf  type 

The  northern  California  continental  shelf  is  typical  of  continental  shelves  on  mu^row, 
swell  and  storm  dominated  continental  margins.  The  region  chosen  for  exercising 
the  model  is  north  of  San  Francisco,  between  Pt.  Reyes  at  38°N  and  Pt.  Arena  at 
39°N  (Figure  6.10.  The  shelf  width  ranges  from  about  16  km  at  Pt.  Arena  to  about 
30  km  just  north  of  Pt.  Reyes.  The  shelf  break  is  slightly  deeper  than  the  world 
average,  at  about  150-170  m.  Bottom  topography  is  simple  and  shore  parallel. 
The  region  is  subject  to  large  wind  stress  variations  and  is  exposed  to  open  ocean 
swell.  Sediment  input  comes  primarily  from  a  single  source,  the  Russian  River,  and 
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Figure  6.9:  Predicted  deposition/erosion  rates  and  reworking  depths  for  H=5m,  T=12sec, 
when  grain  size  distribution  changes  alongshelf  (Figure  6.8).  Transport  of  sand  extends  no 
farther  than  90  m  water  depth;  silt  is  transported  to  a  depth  of  100  m.  All  transport  is 
alongshelf;  there  is  no  cross-shelf  component  of  current,  (a)  Total  erosion  is  up  to  1  cm/day 
on  the  niid-to>outer  shelf,  and  deposition  of  up  to  3  mm/day  occurs  in  the  shallow  parts  of 
the  shelf.  Deposition  results  &om  the  combined  increase  in  stratification  effects  of  waves 
at  shallower  depths  and  the  increased  effect  of  stratification  on  the  larger  silt  component 
at  the  southern  end  of  the  box. 
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(b)  Deposition  and  erosion  rates  by  grain  size  show  that  nearly  all  of  the  net 
transport  is  in  the  silt  size  class,  despite  its  small  bed  concentration.  The  reworking 
depths  show  that  nearly  all  of  the  load  is  silt. 
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Figure  6.10:  The  Northern  California  Shelf,  on  which  the  second  set  of  applications  of 
the  sediment  transport  model  are  based,  showing  the  sand-silt-sand  sediment  distribution 
across  the  shelf  (from  Drake  and  Cacchione,  1984). 

consists  of  about  12%  sand  and  about  44%  each  silt  and  clay  (Klise,  1983).  Most 
of  the  sediment  input  occurs  in  winter  as  a  result  of  seasonal  storms. 

In  cross-section,  the  shelf  can  be  divided  into  three  zones,  distinguished  by 
bottom  slope  and  sediment  type  (Cau:chione  et  al,  1983).  The  inner  shelf  zone  has 
a  bottom  slope  of  1  —  2®,  extends  to  a  depth  of  about  60  m,  and  is  covered  with 
fine-  to  coarse-grained  sands.  The  slope  decreases  abruptly  between  60  and  70  m 
to  about  0.2  —  0.3®.  This  central  shelf  region  is  covered  with  sandy  silts,  with  a 
clay  component  of  up  to  15%.  At  depths  of  130-140  m,  the  silty  sediment  thins  to 
less  than  0.5  m  and  the  slope  increases  to  about  0.5  —  1®.  The  outer  shelf-upper 
slope  region  is  char2u:terized  by  patchy,  sandy  sediment  cover,  of  modern  and  relict 


origin,  to  depths  of  about  200  m  (Klise,  1983). 

The  physical  oceanography  and  meteorology  of  this  region  were  studied  inten¬ 
sively  in  the  Coastal  Ocean  Dynamics  Experiment  (CODE)  from  1981-1983  (CODE 
Group,  1983),  and  studies  of  the  near-bottom  flow  and  sediment  transport  have 
been  conducted  both  as  part  of  CODE  and  independently  (Drake  and  Cacchione, 
1984;  Cacchione  and  Drake,  1982;  Cacchione  et  al,  1987),  These  studies  provide 


essential  first  order  data  needed  for  the  model. 


Circulation  on  the  northern  California  shelf  is  seasonal,  primarily  controlled  by 
steady  winds  from  the  north  in  spring  and  summer  (upwelling  season)  and  cyclonic 
storm  winds  from  the  south  interspersed  with  strong  northerlies  in  winter.  From 
about  April  through  October,  the  shelf  is  subject  to  strong,  steady  equatorward 
winds.  These  winds  control  waves  and  currents,  generating  southerly  currents  with 
bottom  velocities  less  than  15  cm/sec  (Drake,1984)  and  windsea  with  periods  less 
than  12  seconds  and  significant  height  less  than  3  meters  (U.S,  Army  Corps  of 
Engineers,  1981-82).  About  October,  winter  wind  and  circulation  patterns  become 
dominant.  Winds  are  generally  lighter  and  more  variable  than  in  summer,  but 
are  interrupted  by  intermittent  cyclonic  storms  whose  strong  winds  lead  to  near¬ 
bottom  currents  up  to  40  cm/sec  (Drake  and  Cacchione,  1984).  Waves  associated 
with  these  storms  commonly  have  periods  of  14-18  seconds  and  significant  heights 
of  3-5  m  (H.  Graber,  unpub.  data).  During  any  .seaison,  swell  generated  by  distant 
storms  may  affect  the  shelf  with  waves  of  periods  up  to  25  seconds  or  more;  these 
waves  may  persist  for  durations  of  12  to  24  hours  (Cacchione,  pers.  comm.). 

Like  the  Mid-Atlantic  Bight,  the  northern  California  shelf  under  non-storm 
conditions,  which  includes  most  of  the  summer  season,  is  characterized  by  little 
or  no  sediment  motion  (Drake  and  Cacchione,  1984).  Winter  storms  are  cyclonic, 
with  winds  initially  from  the  southeast.  The  strong,  poleward  bottom  currents  they 
generate  typically  persist  for  one  to  three  days.  Transmissometer  observations  of 
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Figure  6.11:  Schematic  diagram  of  grid  representing  the  N.  California  continental  shelf 
and  wave  and  current  conditions  used  in  applications  described  in  text 

suspended  sediment,  here  like  the  Atlantic,  are  strongly  correlated  with  storms 
(Drake  and  Cacchione,  1984). 

Application  One:  Large  storm,  single  grain  sixe  bed 

The  first  application  of  the  boundary  layer  model  addresses  the  case  of  a  single 
grain  size  bed  under  wave  and  current  conditions  typical  of  a  large  storm  on  the 
northern  California  shelf.  The  boxes  begin  at  a  depth  of  30  meters  and  extend 
out  to  the  depth  at  which  sediment  no  longer  is  put  in  motion  by  the  waves  and 
currents  being  tested.  However,  in  order  to  approximate  the  slope  on  the  California 
shelf,  the  grid  spacing  is  one  kilometer  in  water  depths  of  30-60  m,  two  kilometers 
from  60-130  m,  and  one  kilometer  at  depths  greater  than  130  m  (Figure  6.11). 
The  smaller  grid  magnifies  the  erosion/deposition  rate  due  to  net  transport,  since 
the  difference  in  sediment  volume  entering  and  leaving  the  box  at  its  boundaries  is 
being  distributed  over  a  smaller  area  inside  the  box. 


The  wave  used  to  represent  a  storm  condition  is  five  meters  high,  with  a  period 
of  15  seconds  and  a  direction  of  45°  with  respect  to  the  alongshelf  direction.  The 
current  at  one  meter  above  the  bottom  is  30  cm/sec  with  a  direction  of  330°  with 
respect  to  the  alongshelf  direction  (Figure  6.11).  This  represents  a  near-bottom 
current  with  a  26  cm/sec  alongshelf  component  and  a  15  cm/sec  velocity  in  the 
offshore  direction.  This  is  typical  of  near-bottom  storm  current  observations  at 
depths  of  about  50  meters,  although  the  cross-shelf  component  generally  diminishes 
to  about  7  cm/sec  with  increasing  water  depth  (D.  Cacchione,  pers.  comm.).  The 
deeper  estimates  of  erosion/deposition  may,  therefore,  be  excessive  by  a  factor  of 
two  or  more.  This  wave  and  current  condition  was  applied  to  the  same  four  single¬ 
grain-size  shelf  scenarios  as  the  Mid-Atlantic  Bight  storm,  with  grain  diameters 
ranging  from  medium  sand  to  medium  silt. 

The  results  are  shown  in  Figure  6.12.  A  detailed  description  of  the  figure’s 
format  was  given  in  conjunction  with  the  Atlantic  Shelf  application.  The  effects 
of  roughness  and  stratification  were  discussed  in  detail  in  the  Mid- Atlantic  Bight 
results,  so  they  will  receive  only  cursory  mention  here.  The  points  to  note  in  these 
results  are  the  increased  transport  rates  resulting  from  the  increased  depth  gradient 
and  the  increased  transport  depths  due  to  the  larger  wave  period. 

Generally,  all  grain  sizes  show  a  similar  depositional  pattern  zis  wave  bottom 
velocities  and  excursion  amplitudes  decrease  across  the  shelf  (Figure  6.12).  On  the 
inner  shelf,  there  is  relatively  small  volume  erosion,  or  near-equilibrium.  A  region 
of  net  cross-shelf  erosion,  increasing  in  volume  with  depth,  is  consistently  found  on 
the  central  shelf.  The  outer  shelf  is  an  area  of  deposition  for  all  grain  sizes,  with 
the  increase  or  decrease  controlled  by  the  influence  of  outer  Ekman  layer  transport. 
These  trends  mirror  trends  in  reworking  depths,  which  begin  fairly  high,  decrease 
on  the  mid-shelf  where  stratification  effects  are  strong,  incre2ise  again  on  the  outer 
shelf  where  stratification  effects  decrease,  and  decrcMe  toward  no  motion  as  wave 


total  bed-  near  O.E.  d-’0.025rfn 
load  bottom 


Figure  6.12:  Predicted  depoeition/erosion  rates  and  reworking  depths  for  single-grain-size 
shelves,  with  shelf  and  storm  description  derived  from  conditions  typical  of  N.  California 
continental  shelf.  Grid  and  wave-current  conditions  shown  in  Figure  6.11. 

H  -  5m,  T  =  15$,  =  45*,  u,  =  30  em/$ec,0eu,  =  330* 

(a)  Results  for  medium  sand  show  variable  deposition  rates  across  the  shelf,  ranging  as 
high  as  2  mm/day,  with  erosion  or  deposition  of  0.5  mm/day  or  less  being  more  typical. 
Motion  ceases  at  130  m  depths. 
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(b)  Results  for  fine  sand  show  deposition  and  erosion  rates  as  high  as  3  mm/day.  Motion 
ceases  at  130  m. 
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Figure  6.12:  (c)  Results  for  coarse  silt  show  net  transport  rates  generally  increasing  with 
depth,  with  erosion  as  high  as  1.5  cm/day  in  the  central  shelf  and  deposition  as  high  as 
3  cm/day  30  m  deeper. 
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(d)  Results  for  medium  silt  show  net  transport  increasing  with  depth,  changing  from  erosion 
to  deposition  as  the  effects  of  stratification  weaken  with  decreasing  bottom  wave  velocities. 
Deposition  of  as  much  as  14  cm/day  is  predicted  for  the  outer  shelf,  though  the  large  outer 
Ekman  layer  contribution  is  unreliable. 


influence  dies  out  with  depth.  The  location  on  the  shelf  at  which  each  of  these 
transitions  occurs  is  a  function  of  grain  size. 

Deposition  and  erosion  rates  on  the  Pacific  shelf  (Figure  6.12)  exceed  those  for 
a  typical  large  storm  on  the  Atlantic  shelf  (Figure  6.4)  by  at  least  a  factor  of  two 
and  generally  an  order  of  magnitude  or  more.  There  are  three  controlling  factors: 
increased  wave  period,  increased  cross-shelf  component  of  velocity,  and  increased 
shelf  depth  gradient.  The  increase  due  to  the  difference  in  shelf  gradient  is  the 
overriding  factor,  and  is  easy  to  determine.  Since  the  change  in  depth  over  18  km 
on  the  Atlantic  shelf  corresponds  to  1  or  2  km  on  the  Pacific  shelf,  the  net  transport 
rate  for  similar  depths  is  immediately  raised  by  a  factor  of  18  or  9  for  1  or  2  km 
grid  spacings,  respectively.  The  increased  cross-shelf  component  of  current  velocity 
adds  to  the  net  transport  slightly,  but  the  changes  in  stratification  and  roughness 
effects  caused  by  the  increased  wave  period  can  temper  the  expected  increases  due 
to  the  other  two  factors. 

The  same  patterns  of  transport  are  evident  on  both  the  Pacific  and  Atlantic 
shelves,  except  that  the  increased  wave  period  causes  the  shallow  water  {large 
u»)  patterns  to  reach  farther  out  on  the  shelf.  Bedload  transport  is,  again,  of 
minor  importance  relative  to  suspended  tremsport  in  storms.  The  outer  Ekman 
layer  contribution  to  the  net  transport  is  small  for  medium  sand  and,  at  relatively 
shallow  depths,  for  smaller  diameter  sediments  as  well.  The  sudden  decrease  in 
stratification  effects  can  be  tracked  at  90-100  m  for  all  three  smaller  grain  sizes, 
where  increased  effects  of  outer  Ekman  layer  transport  become  noticable.  This 
occurs  because  both  the  height  of  the  Ekman  layer  and  the  sediment  transporting 
capacity  of  the  current  boundary  layer  above  the  wave  boundary  layer  are  decreased 
when  stratification  effects  are  strong.  These  combine  to  enhance  erosion  when 
stratification  is  stronger  at  the  shallower  end  of  the  grid  squ2U'e,  as  is  the  case 
between  90  and  100  m  for  these  three  grain  sizes  in  this  wave-current  setting. 
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The  patterns  of  deposition  and  erosion  are  interesting  with  respect  to  the  ob¬ 
served  sediment  distributions  on  the  Pacific  shelf.  The  net  deposition  of  sands, 
decreasing  to  a  depth  of  60  m,  and  erosion  of  silts  in  the  onshore  boxes  predicts 
the  sandy  nearshore  zone  of  the  California  shelf.  The  predicted  small  volumes  of 
erosion  and  deposition  of  sand  out  to  depths  of  130  m  is  likewise  reflected  in  the 
small  percentages  of  sand  found  in  the  mid-shelf  mud  zone.  The  predicted  pattern 
of  distribution  for  silts  likewise  mimics  the  observed:  there  is  nearshore  erosion, 
then  deposition  out  to  a  depth  at  which  all  sediment  motion  suddenly  ends.  The 
predicted  zones  of  deposition  are  shifted  seaward  by  ~  20  m  relative  to  the  ob¬ 
served,  especially  with  respect  to  the  medium  silt.  However,  if  a  smaller  wave  had 
been  used  -four  meters  rather  than  five  meters  in  height  perhaps-  the  zones  would 
shift  landward.  The  results  suggest  that  the  offshore  limits  of  the  mud  belt  may 
be  probably  deflned  by  seasonal  storm  wave  characteristics. 

Application  Two:  Moderate  swell,  single  grain  size  bed 

This  application  corresponds  to  the  occurrence  of  moderate  swell  on  the  continental 
shelf  on  a  day  of  average  currents.  The  wave  height  is  two  meters;  period  is 
twenty  seconds;  direction  is  out  of  the  South  Pacific  at  45°.  The  current  velocity 
is  15  cm/sec  at  one  meter  above  the  bottom  in  a  direction  330°,  corresponding  to 
a  cross-shelf  component  of  7.5  cm/sec  and  an  alongshelf  component  of  13  cm/sec. 

Because  of  the  20  second  wave  period,  the  maximum  depth  of  sediment  motion 
is  only  30  -40  m  less  than  the  five  meter  waves  of  the  first  application  (Figure  6.13). 
The  patterns  of  deposition  across  the  shelf  are  similar:  deposition  or  steady  state 
holds  in  the  shallowest  reaches;  increasing  effects  of  stratification  in  the  mid-shelf 
lead  to  erosion  for  the  three  finer  grain  sizes;  and  decreasing  stratification  causes 
deposition  on  the  outer  shelf.  The  medium  sand  predictions  follow  the  ‘expected’ 
deposition  pattern  entirely,  demonstrating  negligible  effects  of  either  roughness  or 
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stratification  changes  with  depth.  The  negligible  deposition  and  reworking  rates  for 
fine  sand  and  coarse  silt  suggest  that  storm  transport  dominates  mid-  and  outer- 
shelf  deposition  patterns  for  these  grain  sizes.  In  the  mid-shelf  areas,  however, 
the  deposition  rates  for  medium  silt  is  similar  for  the  storm  and  swell  conditions. 
Stratification  effects  under  the  large  waves  limit  the  transport  of  the  fine  grains, 
leading  to  a  case  where  deposition  patterns  are  as  dependent  on  swell  as  on  storm 
conditions. 
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Figure  6.13:  Predicted  deposition/eroaion  rates  and  reworking  depths  for  single-grain-size 
shelves,  with  shelf  and  swell  description  derived  from  conditions  typical  of  N.  California 
continental  shelf.  Grid  is  shown  in  Figure  6.11. 

H  =■  2m, T  =  20s, =  45°, Ur  =  20  em/Bec,0ew  =  330° 

(a)  Results  for  medium  sand  show  expected  deposition  rates,  demonstrating  that  the 
wave-generated  shear  stresses  are  neither  high  enough  to  wash  out  ripples  nor  to  set  up 
stratification  efiects  strong  enough  to  counter  the  expected  pattern.  Deposition  rates  are 
no  more  than  .2  mm/day.  Motion  ceases  at  90  m  depths.  .  ^  1 
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(b)  Results  for  fine  sand  show  deposition  and  erosion  rates  no  more  than  0.2  mm/day. 
Motion  ceases  at  90  m. 
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Figure  6.13:  (c)  Results  for  co&rse  silt  show  net  transport  rates  generally  increasing  with 


depth,  with  erosion  as  high  as  3.0  tnm/day  in  the  central  shelf  and  deposition  as  high  as 
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(d)  Results  for  medium  silt  show  net  transport  increasing  with  depth,  changing  from  erosion 


to  deposition  as  the  effects  of  stratification  weaken  with  decreasing  bottom  wave  velocities. 
Deposition  of  as  much  as  2.1  cm/day  predicted  for  the  outer  shelf  is  comparable  with  the 
storm  case.  Motion  ceases  at  130  m 
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Chapter  7 


Conclusions  and  Future  Work 


The  purpose  of  this  work  has  been  to  apply  a  theoretical  wave-current  boundary 
layer  model  to  the  prediction  of  sediment  transport  on  the  continental  shelf,  to 
test  the  sensitivity  of  the  model  to  some  of  the  many  parameters  and  assumptions 
necessary  to  represent  the  physical  processes,  and  to  use  the  model  to  predict 
patterns  of  sediment  transport  and  deposition  under  conditions  representative  of 
some  found  on  different  continental  shelves.  This  work  may  be  seen  as  a  step 
toward  the  ultimate  sedimentological  goal  of  producing  a  general  predictive  model 
which  reliably  explains  existing  patterns  of  sediment  deposition  and  distribution, 
and  which  can  be  used  to  characterize  the  depositional  environment  of  elements  of 
the  stratigraphic  record. 

A  physical  explanation  and  a  mathematical  treatment  of  the  fluid  dynamic 
theory  were  offered  separately  so  that  a  reader  unfamiliar  with  boundary  layer 
theory  could  follow  the  explanations  of  the  results.  The  equations,  however,  provide 
precise  demonstrations  of  the  relations  of  the  forcing  and  sediment  response. 

A  presentation  of  the  response  of  the  model  to  variation  in  the  input  parameters 
and  some  of  the  assumptions,  presented  in  the  Chapter  5,  demonstrated  a  number 
of  consistent  results  and  limitations  which  were  stated  in  the  sumn.ary  to  that 
chapter.  Some  of  them  are  repeated  here: 


•  Large  wave  bottom  velocities  (>~  40^),  especially  when  near-bottom  cur¬ 
rents  are  small  (<~  10 ^)j  cause  stratification  effects  which  reduce  load 
predictions  by  50%  or  more.  These  predictions  require  field  verification. 

•  Large  wave  bottom  velocities  occurring  where  the  seafloor  surface  grains  are 
silt  sized  or  finer  may  lead  to  nonunique  solutions  to  the  stratified  bound¬ 
ary  layer  model,  depending  on  specified  reference  velocity.  The  predicted 
sediment  load  and  transport  may  vary  by  an  order  of  magnitude  for  the  dif¬ 
ferent  solutions.  If  the  model’s  representation  of  stratification  is  correct,  this 
nonuniqueness  requires  specification  of  another  boundary  condition. 

•  Large  bottom  currents  (>~  30^  at  z  =  Im)  can  suspend  sediments  to 
heights  in  the  water  column  greater  than  the  near-bottom  layer  on  which 
this  model  focuses.  When  estimates  of  load  and  transport  volumes  in  the 
outer  Ekman  layer  are  greater  than  or  equal  to  the  near  bottom  volumes,  the 
model’s  predictions  of  total  transport  aure  less  certain. 

•  In  a  related  concern,  the  predicted  bottom  Ekman  layer  depth  based  on  mean 
shear  velocity  can  exceed  the  water  colunm  depth  in  strong  current  flows.  If 
so,  the  outer  Ekmain  layer  predictions  should  be  discarded. 

To  test  the  wave-current  boundary  layer  model  in  more  realistic  situations,  it 
was  inserted  into  a  simple  numerical  box  model  and  used  to  calculate  net  transport 
under  conditions  representative  of  two  very  different  continental  shelves.  Data  from 
the  Mid- Atlantic  Bight  were  consulted  to  determine  storm  conditions  representative 
of  that  area.  The  model  was  run  to  test  the  effects  on  sediment  load  and  transport 
predictions  of  intensification  of  wave  bottom  velocities  due  to  cross-shelf  slope  on 
the  mid-to-outer  shelf.  One  might  expect  a  uniform  increzise  in  sediment  load  and 
transport  in  the  onshore  direction,  caused  by  the  increase  in  wave  bottom  velocity 


as  the  wave  shoals.  In  the  presence  of  an  offshore  current,  this  would  generate 
a  prediction  of  net  deposition  at  all  depths,  declining  in  volume  with  increasing 
water  depth.  This  pattern  is  not  observed,  however,  for  two  reasons: 

•  Sediment  stratification  effects  intensify  with  increased  wave  bottom  veloc¬ 
ity.  Predicted  transport  is  sometimes  smaller  at  the  onshore  side  of  the  grid 
square,  leading  to  predicted  erosion  in  the  presence  of  an  offshore  component 
of  flow; 

•  Sediment  roughness  responds  to  the  flow;  in  particular,  larger  wave  bottom 
velocities  can  wash  out  ripples  on  sand  beds,  leading  to  lessening  of  wave¬ 
generated  turbulence.  This  can  lead  to  smaller  deposition  than  expected 
under  the  wave. 

Storm  conditions  typical  of  the  Northern  California  shelf  were  applied  to  a 
narrow  continental  shelf.  The  shelf’s  depth  gradient  itself  aimplifled  the  ero¬ 
sion  and  deposion  predictions.  The  predicted  patterns  of  deposition  were  in 
general  agreement  with  the  observed  sand-silt-sand  patterns  on  the  Califor¬ 
nia  margin.  The  results  suggest  that  the  outer  limit  of  the  mud  belt  found  on 
that  shelf  may  be  defined  by  sezisonal  storm  wave  characteristics,  and  in  this 
area  that  storm  transport  plays  a  more  important  part  in  moving  sediment 
through  the  2irea  than  do  the  average  conditions. 

In  conclusion,  the  analysis  performed  using  the  Grant-Madsen-Glenn  wave- 
current  boundary  layer  model  showed  that  the  model  gives  realistic  predic¬ 
tions  of  sediment  resuspension  and  transport  when  stratification  effects  are 
included.  Stratification  is  an  important  effect  in  this  model,  and  field  tests 
are  needed  to  determine  how  realistic  the  model’s  treatment  of  stratification 
is. 
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7.1  Future  Work 

There  are  a  number  of  questions  raised  by  the  model  which  require  field 
experiment  to  answer.  They  are: 

—  The  model  predicts  intense  sediment  stratification  effects  when  wave 
bottom  velocities  are  high,  resulting  from  large  sediment  concentrations 
inside  the  wave  boundary  layer.  Are  these  real,  or  are  they  artifau:ts  of 
the  model?  How  can  they  be  observed  in  the  storm  wave  environment? 

—  The  model  predicts  very  large  sediment  concentrations  in  the  wave 
boundary  layer.  Are  they  real? 

—  The  model  predicts  reworking  depths  of  sediments  up  to  several  cen¬ 
timeters.  Can  those  be  observed? 

—  The  model  predicts  net  erosion/deposition  baised  on  observable  pairam- 
eters  (i.e.  cross-shelf  components  of  currents).  Can  such  deposition  be 
detected,  or  is  it  buried  in  the  noise  of  other  vaxiables  and  uncertainties? 

—  The  model  represents  the  surface  wave  spectrum  as  a  single,  monochro¬ 
matic  wave  velocity  at  the  bottom.  How  is  the  spectrum  best  repre¬ 
sented? 

If  this  model  is  to  be  useful  in  geological  studies,  these  questions  need  answers. 
Measurement  technology  has  lagged  behind  theory,  so  that  most  models  have 
to  base  their  verification  on  gross  measurements  when  fine  ones  are  needed. 
New  instrumentation  arid  techniques  may  be  applicable,  however,  and  field 
verification  is  unquestionably  the  next  step  in  applying  this  model. 
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Appendix  A 


Modelling  approach 


The  diagram  of  Figure  A.l  outlines  the  approach  driving  the  modelling  effort 
described  in  this  dissertation.  Each  box  contains  a  model  element.  The 
guiding  premise  is  that  each  element  is  autonomous  within  the  model,  and  can 
be  replaced  as  advances  in  the  appropriate  field  outmode  the  approach  being 
used.  The  shaded  boxes  indicate  model  elements  included  in  this  dissertation. 

A  complete  continental  shelf  model  requires  input  from  a  wide  range  of  sci¬ 
entific  disciplines  (Figure  A.l).  Work  is  now  underway  attempting  to  char¬ 
acterize  wave  fields  from  satellite  observations  and  shelf  circulation  from  me¬ 
teorological  and  topographic  data.  At  the  same  time,  studies  are  underway 
to  characterize  the  response  of  sediments  to  shear  stresses  caused  by  flows, 
the  influence  of  bottom  benthos  on  those  responses  and  the  influence  of  the 
responses  (i.e.  bed  roughness  and  suspended  sediment-induced  stratification) 
on  the  bottom  boundary  layer  flows. 

To  clarify  the  interactions  illustrated  in  the  diagram,  the  generation  and 
effects  of  surface  waves  will  be  traced  through  the  diagram.  A  similar  exercise 
for  currents  might  be  a  worthwhile  investment  of  time  for  the  reader:  these 
diagrams  contain  a  great  deal  of  information  on  the  complex  intereurtions  in 
the  shelf  system  and  studying  them  can  provide  some  insight. 


Starting  at  the  top  of  Figure  A.l,  note  that  large-scale  atmospheric  models 
drive  predictions  of  wind  fields,  which  then  can  be  applied  to  the  prediction 
of  wind-generated  waves.  The  elements  of  the  diagram  are  emphjisized  so 
that  it  is  easier  to  place  the  discussion  in  context  of  the  figure.  Depending 
on  the  wave  height  and  period  and  the  water  depth,  the  wave  may  break,  it 
may  generate  oscillatory  motion  to  some  depth  in  the  inviscid  core,  or  it  may 
generate  significant  velocities  (up  to  100  cm/sec  or  more)  at  the  seafloor.  In 
the  latter  case,  the  effects  of  the  wave  can  be  modelled  as  a  wave  boundary 
layer  in  the  bottom  few  centimeters  of  the  water  column.  If  the  bottom 
velocities  are  sufficiently  large,  bottom  sediments  will  be  moved,  perhaps 
forming  ripples.  These  effects  are  included  in  a  movable  bed  roughness  model, 
which  combines  with  other  elements  of  bottom  roughness  to  influence  the 
form  of  the  wave  bottom  boundary  layer.  This  is  the  first  element  of  feedback 
in  the  system.  If  currents  are  present,  the  nonlinear  interaction  of  waves  and 
currents  also  must  be  modelled,  along  with  its  effect  on  the  roughness  and 
the  wave  boundary  layer. (  The  roughness  model  is  in  pwt  a  derivative  of 
semi-empirical  models  of  bedform  development  and  near  bed  transport,  which 
explains  the  connection  to  the  bottom  boundzury  condition  models  of  the 
lower  right-hand  corner.) 

Having  determined  the  bottom  roughness  and  wave  bottom  velocities,  the 
total  boundary  shear  stress  can  be  calculated  using  a  friction  factor  or  drag 
coefficient  model.  Empirical  initiation  of  motion  models  use  the  skin  friction 
component  of  the  shear  stress  to  calculate  conditions  for  sediment  entrain¬ 
ment;  other  semi-empirical  models  are  used  to  calculate  the  bottom  boundary 
condition  for  suspended  load,  the  reference  concentration.  The  total  shear 
stress  is  also  used  in  an  eddy  viscosity  model  of  turbulent  production  auid  ver¬ 
tical  diffusion  of  mass  and  momentum.  The  vertical  diffusion  and  reference 
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concentration  are  combined  in  a  suspended  sediment  model.  Stratification  by 
sediments  can  influence  vertical  diffusion,  in  another  feedback  mechanism, 
and  they  jointly  influence  the  calculation  of  the  mean  velocity  profile.  The 
mean  velocity  profile  has  been  generated  using  models  of  geostrophic  and 
Ekman  layer  flow. 


Appendix  B 


Friction  Factor,  Shear  Stress  and  Shear 
Velocity  Solutions 


The  characteristic  boundary  shear  stresses  and  shear  velocities  are  calculated 
from  the  instantaneous  boundary  shear  stress.  The  instantaneous  boundary 
shear  stress  is  defined  in  GMG  using  a  quadratic  drag  law,  as  defined  in 
Equation  3.21: 


-  1  i.  /  1  2\t  «  V  , 

ro  =  +  t;  )(-  -t) 

‘  (u*4-V*)5  (u*+V*)a 

where  u,v  are  the  i,y  components  of  a  combined  wave  and  current  reference 
velocity  close  to  the  bottom  (though  we  are,  for  the  moment,  assuming  that 
the  wave  and  near-bottom  current  are  collinear  in  the  z-direction) .  /e*  is  the 
combined  wave  and  current  friction  factor.  The  characteristic  shear  stress  in 
the  wave  boundary  layer  (roe*,  =  p  is  defined  as  the  maximum  value  of 
Equation  3.21.  For  the  current  boundary  layer,  Toe  (=  is  calculated  by 
time  averaging  equation  3.21.  The  solutions  for  the  shear  velocities  are: 
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where  a  and  V  j  are  functions  of  the  maximum  and  time-averaged  velocities, 
respectively,  in  the  wave  boundary  layer,  u#  is  a  representation  of  the  velocity 
of  the  mean  flow  in  the  wave  boundary  layer,  so  ^  is  a  representation  of  the 
relative  strength  of  the  mean  versus  the  maximum  oscillatory  flow  in  the  wave 
boundary  layer.  The  value  of  /«„  is  determined  using  these  definitions  and 
the  wave  velocity  profile.  The  value  of  few,  is  calculated  implicitly  from  the 


equation: 


10.097(^)S4-l’  +  210.097(^)i4llAk“*-=  (B.l) 

Ah  f  A  Ah  f*  2a*  4  4a  a 

Jew  Jew 

where  <f>e  is  the  angle  between  the  wave  and  current  directions.  Ah  is  the 
bottom  excursion  amplitude  for  the  wave,  defined  in  Equation  3.26,  and  T 
is  derived  from  the  equation  for  the  wave  velocity  (Equation  B.2)  to  be: 

_ 1 

2^j  (ker’2?|-fkei*2?o^)^ 

The  solution  for  the  wave  momentum  equation  inside  the  wave  boundary  layer 
is  not  explicitly  of  interest  for  the  present  problem,  though  it  is  necessary  for 
the  calculation  of  the  boundary  shear  stress.  The  solution  immediately  above 


the  bottom  is: 


..  f,  ,  1  ^  +  1*154  +  »| 

Uw  —  tl6[l  +  n  I  T*  J 

^ker2^|  +*kei  2^| 


(B.2) 


where  ^  =  -f-  and  £o  =  Ker  and  kei  are  Bessel  Functions:  tabulated 
solutions  to  a  particular  form  of  differential  equation.  The  derivation  and 
background  for  the  wave  velocity  profile  and  friction  factor  equation  are  cov¬ 
ered  in  some  detail  in  Grant  and  Madsen,  1979. 
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Appendix  C 


Model  Results:  Five  More  Wave  Cases 


Plots  of  results  for  five  wave  conditions  discussed  in  the  Sensitivity  section 
(Section  5)  are  included  here.  The  list  of  wave  and  sediment  conditions 
represented  in  these  plots  is  presented  in  Table  5.2  on  page  88.  The  five  wave 
conditions  represented  are:  moderate  windsea,  laurge  storm  with  small  fetch, 
moderate  swell,  large  storm  with  large  fetch,  and  extreme  swell. 
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Figure  C.3:  Predicted  near*bottom  load  and  transport  (cm’/cm/sec)  values  for 
coarse  silt  in  moderate  windsea.  ranges  from  7-8  cm;  |  ranges  from  6  to  18  m. 
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Figure  C.6:  Predicted  near-bottom  load  (^)  and  transport  (cm’/cm/see)  values  for 
very  fine  sand  in  a  large  storm,  early  stages.  ranges  from  9-10  cm;  |  ranges  from  7  to 
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Figure  C.9:  Predicted  near-bottom  load  and  transport  (cm*/cm/sec)  values  for 

fine  sand  in  moderate  swell.  T=20  sec,  and  the  wave  height  in  SO  m  of  water  would  be 
2.2  m.  Ur  is  reference  current  velocity.  5*  ranges  from  29  to  34  cm;  g  ranges  from  12  to 
39  m. 
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Figure  C.18:  Predicted  near-bottom  load  (^)  and  transport  (cm’/cm/»cc)  values  for 
very  fine  sand  in  extreme  swell.  6^  ranges  from  29-40  cm;  |  ranges  from  10  to  54  m. 
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m  Abstract  (Limit:  200  words)  ~ 

^  This  thesis  presents  an  application  of  the  Grant-Madscn-Glcnn  bottom  boundary  layer  model  (Grant  and  Madsen,  1979;  Glenn  and 
Grant,  1 987)  to  predictions  of  sediment  transport  on  the  continental  shelf.  The  analysis  is  a  two-stage  process.  Via  numerical  experiment, 
we  explore  the  sensitivity  of  sediment  transport  to  variations  in  model  parameters  and  assumptions.  A  notable  result  is  the  enhancement 
of  suspended  sediment  stratincation  due  to  wave  boundary  layer  effects.  When  sediment  stratification  is  neglected  under  conditions  of 

large  wave  bottom  velocities  (i.e.  u^  >  40  cm/sec)>conceniration  predictions  can  be  more  than  an  order  of  magnitude  higher  than  any 
observed  during  storm  conditions  on  the  continental  shelf. 

A  number  limitations  to  application  emerged  from  the  analysis.  Solutions  to  the  stratified  model  are  not  uniquely  determined  under 
a  number  of  cases  of  interest,  potentially  leading  to  gross  inaccuracies  in  the  prediciton  of  sediment  load  and  transport.  Load  and  sediment 
tran.sport  in  the  outer  Ekman  Layer,  beyond  the  region  of  emphasis  for  the  model,  can  be  as  large  or  larger  than  the  near-bottom  estimates 
in  some  cases;  such  results  suggest  directions  for  improvements  in  the  theoretical  model. 

In  the  second  step  of  the  analysis,  we  test  the  ability  of  the  model  to  make  {redictions  of  net  sediment  transport  that  arc  consistent  with 
observed  sedi  ment  deposi  tional  pattems^ata  from  the  Mid- Atlantic  B  ighi  and  the  Northern  California  coast  are  used  to  dc  fine  reasonable 
model  input  to  represent  conditions  on  iwVifferent  types  of  shelves.  In  these  examples,  the  results  show  how  the  intensification  of  wave 
bottom  velocities  with  decreasing  depth  Can  introduce  net  transport  over  a  region.  The  patterns  of  erosion/dcposiiion  arc  shown  to  bo 
strongly  influenced  by  sediment  stratification  and  moveable  bed  roughness.  Also  predicted  by  the  applications  is  a  rapid  winnow  Ing  out 
of  fine  grain  si/c  components  when  there  is  even  a  small  variation  of  bed  grain  size  texture  in  the  along-flow  direciion 
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